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Abstract 



We investigate experimental consequences of the postulate that fundamentally photon 
propagation is governed by an SU(2) rather than a U(l) gauge principle. In the context of 
thermodynamics the SU(2) Yang-Mills theory is assumed to be in its deconfining phase. In 
addition to the results already available, we compute the dispersion law for the longitudinal 
mode, which does not thermalize with matter, at one-loop accuracy. We analyze radio- 
metric and bolometric methods to detect a gap in the black-body radiation spectra for low 
frequencies and at low temperatures. We simulate the expected experimental results for 
total-power and Dicke-switch radiometry. These simulations indicate that the predicted 
SU(2) signature is measurable with presently available detector technology appealing to 
both radiometric methods. We also investigate non-homogeneous thermodynamics (ther- 
momagnetic effect) in the adiabatic approximation. 
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Chapter 1 

Introduction 



The discovery of the Cosmic Microwave Background (CMB) by Penzias and Wilson [52] is 
easily one of the pivotal moments in the history of cosmology and the development of the 
Hot Big-Bang model. COBE further confirmed that the CMB power spectrum follows to 
a high accuracy ~ 10~ 2 a Planck distribution in the frequency region of z/=60-600GHz, 

gang EE]. 

Recently, however, measurements of CMB line temperatures T performed by the 
ARCADE 2 instrument |6T] at very low-frequencies (v = 3 GHz ... 90 GHz) show a statis- 
tically significant (5 a) unexplained power excess. Combining the ARCADE 2 data with 
data from the relevant literature in the frequency region from 22 MHz to 10 GHz results 
in a background power law spectrum of the form [1] 



T(u) = To + T R 1-1 , (1.1) 

with T = 2.725 K, v = 1GHz, T R = 1.19 ± 0.14 K, and a spectral index of f3 = 
—2.62 ± 0.04. It is further argued in |14] that this deviation from the perfect black-body 
relation T(v) = const is neither an artifact of galactic foreground subtraction, nor can it 
be attributed to distant point sources. 

Remarkably, a thermodynamical field theory developed in |25[ [28] could be able to 



explain the power excess in (1.1). Namely, it is argued that the gauge group U(1)y of 
the photon field in the Standard Model (SM) is not fundamental but merely emerges from 
a SU(2) Yang-Mills theory in a specific thermodynamical phase: The deconfining phase. 
This is called the SU(2) CMB hypothesis. The theory is determined by a single parameter, 
and can be fixed by specifying the critical temperature T c below which the thermal SU(2) 
gas undergoes a phase transition to the preconfining phase. 

Assuming the SU(2) CMB hypothesis, the excess power at low frequencies measured 
by the ARCADE2 radiometers can be explained as a result of evanescent photon fields 
thermalized near the critical temperature T c = T CMB of the deconfining phase |30j . 

That such an explanation for the excess power of the CMB is possible based solely 
on gauge symmetry considerations and finite temperature effects motivates the search 
for further experimental evidence in favor of the SU(2) CMB hypothesis. The ARCADE2 
observation serves to fix T c , which we use to predict further observable consequences that 
can then be tested by experiment. 

The purpose of this thesis is the investigation of possible realizations of such an ex- 
perimental test. We discuss the characteristics of the thermal radiation spectrum at low 
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temperature and frequency under the SU(2) CMB hypothesis. It will turn out that SU(2) 
gauge dynamics manifests itself through the appearance of a temperature-dependent fre- 
quency cutoff v*(T) below which photon propagation is forbidden. The cutoff frequency 
v* thus defines a screening region. At low frequencies and temperatures comparable to 
that of T c , this effect results in a black-body anomaly: The vanishing of the power out- 
put of a black-body cavity within the screening region. This will turn out to be the key 
experimental effect to search for when probing for SU(2) physics. 

This thesis is organized as follows. Chapter [2] reviews elementary properties of gauge 
theory and corresponding symmetry groups. A short review of Maxwell electrodynamics 
introduces the required mathematical formalism which is then generalized to nonabelian 
gauge symmetries. Basic concepts of topological configurations in SU(2) Yang-Mills theory 
are introduced. In the final section, thermal field theory is shortly summarized. Chapter [3] 
serves as an introduction to our main subject: The thermodynamics of SU(2) Yang-Mills 
theory, in the context of the nonperturbative approach developed in |25l I28j . After a re- 
view of the relevant theory, we proceed to calculate the dispersion law of the transverse 
and longitudinal polarization modes of the SU(2) photon. We also investigate the thermo- 
dynamical consistency of one-loop resummed energy density and pressure. In chapter [4] we 
investigate the possibility of generalizing an important assumption made in chapter [3] while 
deriving the thermal ground state, namely thermodynamical homogeneity (and isotropy). 
It will turn out that, within the restrictions imposed by an adiabatic approximation, cer- 
tain deviations from homogeneity are indeed allowed. Chapter [5] is the most important 
section of this work. After a quick introduction to the basics of radiation transfer, antenna 
theory and radiometry, we investigate the implications of the SU(2) CMB hypothesis on the 
radiometric measurement of thermal photon radiation. We show how the existence of a 
screening region can be established through purely radiometric means. In the final section, 
we comment on the possibility of a bolometric experimental setup. The appendices contain 
the details of certain derivations. Essential for chapter [4] is the result derived in appendix 
|A| where a particular solution to the SU(2) Yang- Mills equations in Minkowskian signature 
is derived. We have attached a CD with the Mathematica source code used in this work. 

Unless stated otherwise, we always work in natural units h = c = ks = !■ 



Chapter 2 

Gauge Fields and Symmetries 



In this chapter, we give a very summarized survey of gauge theory and related concepts, 
in order to provide the reader with the minimal background needed for chapter [3] We 
first review the basics of gauge theory and symmetry groups, then we briefly show how 
Yang-Mills theory emerges by extending the U(l) symmetry to the nonabelian gauge group 
SU(N). Next, we discuss topological field configurations in the context of Euclidean Yang- 
Mills theory. Finally, a brief summary of thermal field theory is given for completeness. 



2.1 Gauge theories and symmetry groups 
2.1.1 Maxwell electrodynamics 

In the Lagrangian formalism, the dynamics of the photon field in the presence of an 
external charged field of current density are described by the Maxwell Lagrangian density 

C M = --F^F^-ej^. (2.1) 

where the coupling constant e characterizes the strength of the electromagnetic interaction. 
It is identified with the electric charge of the current inducing external field. 

The field strength tensor F^ v can be defined by the action of the commutator [D^, D v ] 
on an arbitrary Lorentz scalar 

F^§ = i[£>„,D„]*, (2.2) 

where the covariant derivative is defined in terms of the gauge potential 

D„ = d„ - ieA^ . (2.3) 

The familiar electromagnetic vector fields E and B are recovered from the field strength 
tensor by the identities 

E i = -F 0i , (2.4) 
B i = -l<**F jk . (2.5) 
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It can be shown (Bohm-Aharonov effect) that the gauge potential A^, not the field strength 
tensor F^, is the fundamental physical field of the theory |57] . 

The Maxwell field strength tensor can in turn be written in terms of the gauge potential 



Afj, by inserting equation (2.3) in (2.2) 



F, 



d»A v 



duA L 



(2.6) 



Note that the electromagnetic coupling e does not appear in the expression (2.6), since the 
single gauge field commutes with itself and terms such as e 2 [A^,A u ] vanish. The photon 
cannot interact with itself, and does not carry charge. We will later generalize Maxwell 
electrodynamics to allow for a selfinter acting (also charged) gauge potential multiplet. 
The inhomogeneous Maxwell equations are given by the Euler-Lagrange equations of 



motion derived from (2.1) by the principle of least action 



d^=f. (2.7) 
The homogeneous Maxwell equations follow from the Bianchi identity [3j [15j ISZ] 

dpF^ = , (2.8) 



where the dual field strength tensor F^ u is defined by 
pixu _ Z € l J - l ' a Pp (Tp , 

Here e^ uap is the totally antisymmetric Levi-Civita tensor with e 



(2.9) 



0123 



1. 



The dynamics of the system described by equation. (2.7) and (2.8) are invariant under 
local field shifts 



A„ 



A' 



A„ + -d^, 
e 



(2.10) 



with £ = an arbitrary Lorentz scalar function. Transformations of the form (2.10) 



are called gauge transformations. They leave the field strength tensor (2.6), and thus the 



equations of motion (2.7) and integrability conditions (2.8), unmodified 
pv v p'pv = pv v , 



(2.11) 



This gauge invariance is a manifestation of redundancy in the degrees of freedom of the 
theory: The photon is a massless spin-1 particle with only two transverse polarizations, 
yet in (2.1) the dynamics are given by a 4-component gauge potential An. The excess 



polarization modes are nevertheless necessary for a covariant formalism. Local gauge free- 
dom also implies that the gauge potential transforms under Lorentz rotations as a 4-vector 
only up to gauge transformations. A Lorentz transformation between two frames is always 
associated with a gauge transformation |5T 



In the framework of differential geometry, the field A^ with transformation law (2.10) 
can be interpreted as the connection on a principal \J(l)-bundle over the space-time man- 
ifold M with Minkowskian signature (+,—,—,—), U(l) being the group of lxl unitary 
matrices [21 [15]. The Maxwell gauge potential A^ is therefore also called the U(l) gauge 
field. Its emergence can be understood as the result of demanding invariance under local 
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phase transformations of the quantum state vector characterizing the physical system 
[57] 

|V) -> W) = Sl(x)\rf>), (WO = <V#>- (2-12) 

Here, fl(x) = e l ^ x ^ is a space-time function taking values in U(l), which we write for 
simplicity as £l(x) G U(l). Theories which exhibit invariance under a continuous group 
of local symmetry transformations are called gauge theories. The group U(l) is the gauge 
group of the Maxwell theory. After a U(l) phase transformation, the Maxwell gauge field 



An transforms as in (2.10), which may be rewritten as 



Ap -»■ A'p = fi^fit + -fiS^fit , (2.13) 



for any Q(x) £ U(l). 

The geometric formalism just described can be generalized to larger groups G, while 
respecting Lorentz symmetries and covariance [3J. From ( |2.12 ) it is clear that G must be 



either unitary or anti-unitary. Following Yang and Mills [ 71] . we replace phase transforma- 
tions with non-abelian isotopic spin rotations of SU(2) valued fields. The theory of SU(N) 
groups is reviewed in the next section. 

2.1.2 Lie groups and Lie algebras 

The U(l) group underlying Maxwell electromagnetism is the abelian group of phase ro- 
tations (two arbitrary phase rotations commute), and a generic element of U(l) may be 
written as with real parameter £. In general, a continuous group G of order n may 
be parametrized by a set of n real parameters. A Lie group G of order n is a continuous 
group where the set of parameters live in a differentiable manifold M n . The operations of 
group multiplication and inversion are smooth maps in M n . A compact Lie group is a Lie 
group with compact manifold. 

In the neighbourhood near the identity element 1, M n looks locally like W 1 with flat 
Euclidean metric. This region can always be given a coordinate basis. The n coordinate 
unit vectors t a , a = l..n are called the generators of the group G. Every element g of G 
near the identity can thus be written as 

g = l + iu a t a , (2.14) 

with uj a infinitesimal parameters. Note that for unitary groups U(N), t a is hermitian, and 
for special unitary groups SU(N), t a is also traceless. Elements further away from the 
identity are generated by the exponential map e iw *°. In the case of a matrix group, the 
exponential map is given by the usual series 



°° „-71 



n=0 



Let gi and 32 be two elements close to the identity, from (2.14) we write them as 
gi = 1 + iojft a and g<± = 1 + ibi\t a . Their product is given, up to first order, by the 
expression g\g2 = 1 + i(ujf + w^)^ + 0(w 2 ), so up to first order they both commute, 
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[0i>52] = 0(uj 2 ). Up to second order, however, the expression Q\g2Q\ 1 ^ s n °t the 
identity, but 



9ig29i l g 2 l = 1 - uifoj b 2 [t a , t b ] . (2.16) 



Nevertheless, (2.16) is still close enough to the identity such that it can be written as 
in (2.14). This means that the commutator [t a ,i b ] must be a linear combination of the 



generators t 

[t a ,t b ] = if abc t c . (2.17) 



The factors f abc £ M are called structure constants, and can always be chosen to be 



completely antisymmetric. The generators t c together with the commutator relations (2.17) 
define the Lie algebra of the group G, which we denote as q. The rank of the Lie algebra is 
defined as the number of mutually commuting and linear independent generators. These 
generators form the Cartan sub-algebra of the group. A Lie group G containing no invariant 
subgroup other than the identity and G itself is called simple. A Lie group G is called 
semi-simple if it does not contain any abelian invariant subgroup other than the identity. 

Finally, we review some concepts from the theory of representations |19j . A linear 
representation R of a group G is a map which associates every element g of G with a linear 
transformation R(g) acting on a vector space V of dimension m. For a given basis of V, 
the transformations R(g) can be written as matrices, also of dimension m. The elements 
of a matrix group, for example, build a representation themselves, called the fundamental 
representation of the group. For a Lie group with n generators, the structure constants 
also define a matrix representation, the adjoint representation, made up of n x n matrices 

(Rc)ab = ~if abc ■ (2.18) 

Further discussion of Lie groups, Lie algebras, and their representations can be found in 
[T9] , In the following, we specialize for the group SU(2). 



2.2 SU(2) Yang-Mills Theory 

Following Yang and Mills |71j . we extend the gauge symmetry of the vacuum (j^ = 0) 
Maxwell equations to the non- abelian compact Lie group SU(2). The group SU(2) is defined 
as the group of 2 x 2 unitary matrices with unit determinant. It is a continuous Lie group 
of rank 1, and its the Lie algebra is defined by three traceless hermitian generators which 
we may choose as t a = \o a , a = 1..3, where a a are the Pauli matrices. Then tvt a t b = \$ ab , 
and the structure constant is the totally antisymmetric tensor f abc = e abc . The matrices 
t a define the 2-dimensional fundamental representation, while the 3-dimensional adjoint 



representation is given by the structure constant as in (2.18). 



Define now the Yang-Mills Lie-algebra valued connection A^ 

A,{x) = A-(x)t a , (2.19) 

where the A^(x) are (Minkowski) space-time-dependent real functions. The gauge poten- 
tials A a ^ now carry space-time index ji and internal symmetry index a. They are the basic 
dynamical variables of SU(2) Yang-Mills theory. 
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Generalizing (2.13), we demand that under SU(2) gauge transformations, the Lie- 
algebra valued potential transform as 

A„ -> A'p = SlA^tf + -fifyfit , (2.20) 

where the local transformation Q(x) = e iojata is now an element of SU(2), Q,(x) 6 SU(2), 
and g is the Yang-Mills gauge coupling. Note that this implies that the individual gauge 
potentials A^ = 2tr A^t a transform in the adjoint representation, up to an inhomogeneous 
term. A field (j) a transforming homogeneously cj) a — > <p la = fi^y (f> b under the adjoint 
representation = 2tr t a Q t b is also called an adjoint field. 

As in Maxwell electrodynamics, we introduce the SU(2) gauge covariant derivative D^. 
The operator is constructed so that a field and its covariant derivative transform in the 
same way. Acting on a field eft that transforms bi-homogeneously under the fundamental 
representation, (j) — > Qcf)tP , acts as a SU(2) matrix 

Dp = l 2X 20,i - igAp . (2.21) 
This results in a bi-homogeneous transformation law for D^(j) 

-> QD^rt . (2.22) 
Similarly, on an adjoint field (fr a , acts as 

(D^) a = [D^ <P] a = d^ a + gf abc A b ^ c , (2.23) 
so that the transformation for the adjoint field {D^(j)) a reads 

(D,^ KliiD^f. (2.24) 



A Lie- valued Yang-Mills field strength tensor F^ v can be defined as in (2.2) by its action 
on a field <f> living in the fundamental representation 

F^(f> = -[D^Djt. (2.25) 
9 



Calculating the commutator in (2.25) gives 

F^ v = d„A„ - d v A„ - ig [Ap, A v ] , (2.26) 
written in components 

F% = d,A a v - d v A% + gf ahc A\A% . (2.27) 



Observe that, compared with the strength tensor (2.6) of the U(l) theory, expression 
(2.27) has an extra term which is antisymmetric and bilinear in the gauge fields A^, and 
proportional to g. It describes gauge field self-interactions, and is what makes non-abelian 
gauge theories fundamentally different from their abelian Maxwell counterpart. 

By construction, the Yang-Mills field strength tensor F^ v transforms in the fundamental 
representation 



(2.28) 
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This allows us to define a gauge invariant pure Yang-Mills SU(2) Lagrangian 



C*Y M — 



^ Tpafiu rpa 

4 

--trF^F, 
2 



(2.29) 



The term pure refers to the absence in (2.29) of any field other than the SU(2) gauge 



field. From the behaviour of the field strength tensor under gauge transformations (2.28) 



it is easily seen that the Lagrangian defined in (2.29) is SU(2) gauge invariant. Further 



Lagrangian (2.29) respects Lorentz invariance, renormalizability and CP-Invariance. A 



term proportional to F^ u F apu , where the SU(2) dual field strength tensor is defined as in 



( ]2~9l ), F a>iu = \e pvpa F apa , may be added to the SU(2) Lagrangian ( |2.29[ ) if the demand for 



CP-Invariance is relaxed, but this will not concern us in this work. 
The SU(2) Yang-Mills action 

S Y m is defined as 



d x Cx 



(2.30) 



The Euler-Lagrange equations of motion can be derived from (2.30) by the principle of 
least action. We obtain the SU(2) Yang-Mills equations of motion 



DpF 1 ™ = . 



(2.31) 



As in the Maxwell case, the Bianchi identity of differential geometry results in the integra- 
bility conditions j3[ HI] 



DyF*" = . 



Note that (2.31) follows from (2.32) if the field strength tensor satisfies 



rpa rpa 

in which case F^v is called selfdual, or 



pa 



(2.32) 



(2.33) 



(2.34) 



in which case F pu is called anti- selfdual. Because of their importance, such field configu- 
rations are discussed in the next section. 



2.2.1 Selfdual field configurations 

We have seen that for (anti) selfdual field configurations 



pa 



±F a 

fJ.V > 



(2.35) 



the SU(2) Yang-Mills equations (2.31) are automatically satisfied, since (2.31) then follows 
from the geometric identity (2.32). This is significant, since equation (2.35) is first-order, 



yet (2.35) automatically implies the second-order Yang-Mills equations. 



2.2. SU(2) YANG-MILLS THEORY 



9 



There are no known solutions for the selfdual equation (2.35) in 4D Minkowski SU(2) 



Yang-Mills theory. Consider now Euclidean 4D Yang-Mills, which is arrived at by letting 
the metric become g^ u — > —5^ u and rotating the real time-like component of the 4-vectors 



x 

a a 
A 



IX 



-%A\. 



(2.36a) 
(2.36b) 



with real x 4 , A\. It can be shown that the Euclidean Yang-Mills action Sy M in R (that 



is, the action (2.30) after rotation (2.36)) can be decomposed as 



1 



d 4 x 



±F a^pa v + 1 



pa fpa 



(2.37) 



This is the celebrated Bogomol'nyi decomposition [6]. The action (2.37) is invariant un- 
der space-time translations. The canonical Euclidean energy-momentum tensor 9^ u then 
follows from Noether's theorem [20] as the Noether current resulting from translation in- 



variance, and can be calculated from the Lagrangian (2.29) after rotation (2.36) (up to a 
4-divergence term). In the case of pure Yang-Mills theory it can be written as |55| 



1 



F, 



F^x) (F uX + F uX ) +(ix+±v 



(2.38) 



Note that for (anti) selfdual field configurations the energy-momentum tensor (2.38) 
vanishes. This guarantees that (anti) selfdual field configurations have zero energy-momentum, 



and do not propagate. Further, the Euclidean action (2.35) becomes 



1 



(2.39) 



with +/— sign for selfdual/anti-selfdual field configurations. 

To further interpret this result, rescale the field strength by absorbing the coupling 

i pa 



constant in the gauge field: gA^ 
to 



A„. Then F« v 



9 V 



and the expression (2.39) goes 



1 



d 4 iF a %. 



(2.40) 



Consider now field configurations with finite action (2.40). Following [5], we assume finite 
action configurations define a field strength which vanishes at infinity. That is, F a ^ y {x) —> 
sufficiently fast as |x| — > oo. The gauge potential A^{x) of finite action configurations 
should therefore also tend to zero. But because of gauge freedom, this limit is zero only up 
to a gauge transformation ( 2.20[ ) 



lim A^(x) 
\x\— >oo 



n( x )d fl n j (x) 



(2.41) 



The gauge transformation Q,(x) need only be defined at large say over a sphere S 3 at 
the boundary of Euclidean 4-space-time. In this sense, £l(x) defines a continuous map 



SU(2) 



(2.42) 



But such a map has a well defined integer invariant: The Pontryagin index (topological 
charge) k which counts the number of times the map £1 wraps around the SU(2) manifold. 
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Thus, in Euclidean Yang-Mills theories, asymptotically flat potentials can be distinguished 
by its Pontryagin index. It is a topological invariant of the field configuration: Smooth 
deformations of the gauge field away from the Euclidean boundary S 3 cannot change the 
value of k. For the group SU(2), the index k is given by the gauge invariant expression [3] 



8Tr 2 k 



1 



d 4 xF a ^i> 



/IV 



(2.43) 



For (anti)selfdual field configurations, the topological invariant (2.43) coincides exactly 



with the Yang-Mills action (2.40). In general, the Euclidean Yang-Mills action (2.37) is 



given by the topological term (2.43) plus a sum of bilinears in the strength field. The 



minimal- action configuration within a given topological sector of topological charge k is 



thus given by (anti) selfdual field configurations with action S 



2.2.2 Topological classification 

Here we introduce some terminology we will make use of later while discussing SU(2) Yang- 
Mills thermodynamics in chapter [3] Following Gross, Pisarski and Yaffe |21J . classical gauge 
field configurations in Euclidean space-time with finite action can be classified according 
to topological charge, holonomy and magnetic charge. 



Topological charge 

As we saw in the preceding section, a well defined index k, topological in nature, can be 
defined as 



8vr 2 A; 



d 4 x F^F^ 



(2.43) 



This expression divides different field configurations into topological sectors classified by 
the integer k. 

Let us summarize here the results of the preceding section. In each topological sector 



k, the configuration with minimal Euclidean action (2.37) is given by the solutions to the 



selfdual equations (2.35). The Yang-Mills equations of motion are automatically satisfied 



for selfdual fields, thus reducing second order field equations of the type (2.31 ), (2.32 ) to the 



first order selfdual equation (2.35). Selfdual field configurations in Euclidean space-time 



cannot propagate, since they have a vanishing energy- momentum tensor (2.38) 



Holonomy 

The holonomy of a field configuration is given by the eigenvalues of thePolyakov loop 
defined as a time-like Wilson loop L 



L(x) = V expi / dr A^t, x) , 
Jo 



(2.44) 



where V denotes path ordering and A4 is the Euclidean-rotated gauge field (2.36). 



Magnetic charge 

The magnetic charge is defined as the coefficient m in the expression for the magnetic flux 



B 



m 



|if< 



evaluated at spatial infinity [8]. 
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2.3 Thermal field theory 



Recall from quantum statistical mechanics that the grand canonical partition function Z(f3) 
for a quantum system, which is defined by the Hamilton operator H and vanishing chemical 
potential /i = 0, can be written as an integral over a complete set of eigenfunctions \q) of 



H 



tr e 



-/3H 



dq(q\e-P H \q) 



(2.45) 



where as usual f3 = 1/T is the inverse temperature. It can be shown |35| [42J that by 



rotating to Euclidean space-time (2.36), the partition function (2.45) for a quantum field 
q(x) can be written as a transition amplitude by integrating over all field configurations 
that periodic in Euclidean time. In the case of pure SU(2) Yang-Mills theory, the role of 
q(x) is played by the gauge field A^(x), rotated to Euclidean time. The pure Yang-Mills 
partition function is then given by 



Z 



VA^e^YM 



(2.46) 



A M (0,x)=A M (/3,x) 



where the field configurations are Lie-Algebra valued A^ = A^t a with group generators 



t a normalized as usual, and the Euclidean Yang-Mills action Sy M was defined in (2.37) 
Integration (2.46) goes over all periodic field configurations A^t = 0, x) = A^(t = j3,x). 



In this formalism, the Euclidean space-time M goes over to the Euclidean torus 
5i x R 3 . Selfdual solutions on M 4 — > S\ x R 3 are called calorons. 



Expression (2.46) is used to derive other thermodynamical quantities such as the pres- 
sure, energy, etc. Perturbatively, the partition function InZ can be expanded in bubble 
diagrams. Further details can be found in the literature |35[ I42j . 



Chapter 3 



SU(2) Yang-Mills Thermodynamics 
in the Deconfining Phase 



In this chapter, we review a nonperturbative approach to the thermodynamics of the SU(2) 
Yang-Mills gas developed in |25l I29j , It is shown in |25| that such an approach results 
in three distinct thermodynamical phases: confining, preconfining and deconfining. Each 
phase is characterized by the emergence of a scalar field describing the thermal ground state 
dynamics. In this thesis we are interested in the physics of the deconfining phase, and the 
possibility of experimental detection of ground state effects in this regime. Experimental 
evidence for SU(2) dynamics hinges on the appearance of a frequency cutoff v* below 
which thermal propagation of the SU(2) photon is forbidden. This effect can be detected 
by measuring the antenna temperature of an antenna immersed into a thermalized photon 
gas. The radiometry of the SU(2) photon is discussed in chapter [5] 



This chapter is structured as follows. In section [3T] we briefly review the nonperturba- 
tive approach to the deconfining phase, outlining the basic assumptions and results. Section 
3. 2 1 introduces the resulting effective theory and some of its consequences, for example, the 



emergence of nonzero vacuum energy. In section 3.3 we list the basic thermal excitations 
of the theory and give the Feynman rules necessary to calculate effective radiative correc- 
tions. The purpose of this introductory review is not to give a complete description of the 
theoretical approach in |25|. [29] , but to introduce the concepts and mathematical results 
necessary for an understanding of the physics behind the v* cutoff. 



In section 3.4 we review the hypothesis that thermalized photon propagation is des- 
cribed by the tree-level massless mode of a Yang-Mills SU(2) theory in the deconfining 
phase. This hypothesis fixes the critical temperature of the deconfining phase to be slightly 
lower than the baseline CMB temperature T CMB = 2.726K. 



Our own results are presented in the last two sections. In section 3.5 we calculate radia- 
tive corrections to the propagation of massless modes, both transverse and longitudinal. A 
formula for the temperature-dependent frequency cutoff u*{T) of the transverse polariza- 
tion is derived, and the longitudinal polarization is shown to be nonpropagating. Finally, 



section 3.6 is dedicated to a discussion of thermal resummation and thermodynamical 



consistency. 
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3.1 A nonperturbative approach to Yang-Mills thermodynamics 

The sector of topologically nontrivial, selfdual, nonpropagating pure Yang-Mills configura- 
tions on S± xR 3 Euclidean space-time (calorons) is in general not accessible by perturbation 
methods. The quantum amplitude A of such configurations is non-analytic in the coupling 
constant g, A ~ exp (— ^V), and thus missed at all orders of perturbation theory. 

A thermodynamical approach to this problem was developed in |25| [29] that solves 
this difficulty by integrating out topologically nontrivial degrees of freedom through a 
spatial coarse-graining over such field configurations. This results in an effective theory 
containing topologically trivial excitations interacting with an effective macroscopic scalar 
field (f) a living in the adjoint representation of the gauge group SU(2). The field (fi a also 
introduces with itself a mass scale A which is the only parameter of the theory. 

The effective theory features a thermal ground state above which topologically trivial 
excitations fluctuate. In order to construct such an effective theory two basic constraints 
for infinite- volume thermodynamics should be obeyed |25j . First, in the absence of ex- 
ternal perturbations, any nonzero field resulting from the ensemble average and spatial 
coarse-graining over nonpropagating field configurations must be, (modulo gauge trans- 
formations), spatially homogeneous and isotropic. Furthermore, such a field can only be 
constant over all space, and rotationally invariant. Similarly, the gauge invariant field 
resulting from such averaging and coarse-graining operations cannot explicitly depend on 
time. The basic considerations leading us to the unique definition of a macroscopic scalar 
field are sketched in this section. 



3.1.1 Basic roadmap 

The nonperturbative approach to SU(2) Yang-Mills Thermodynamics developed in [25j 
129] uniquely derives an effective field theory starting from the fundamental Euclidean 
Yang-Mills action considered in the thermodynamical limit of infinite volume. Here, we 
summarize the main line of analysis underlying the emergence of a macroscopic scalar field 

BE- 



Consider the topologically nontrivial, selfdual sector of periodic Euclidean field con- 
figurations. Independently of their topological parameters, these field excitations 
have a vanishing (classical) energy-momentum tensor, and thus do not describe prop- 
agating degrees of freedom. 

Let a Lorentz invariant gauge multiplet cj) a be defined by coarse-graining over such 
configurations in [T] The multiplet (f) a is thus by construction nonpropagating (recall 
section 2.2.1). But this means (j) a acts as a static, spatially homogeneous background. 



Spatial isotropy would be broken by such a Lorentz tensor unless <j) a were a scalar 
field. Further, spatial homogeneity requires that the scalar field <p a have a constant 
modulus \(f)\. Since (fi a by construction cannot propagate and emerges by coarse- 
graining over periodic configurations in Euclidean time, its equation of motion pro- 
duces - in a given gauge - an Euclidean-time periodic phase t/> a (r + /3) = <^ a (r) with 
no space dependence. Also, since the classical action of selfdual configura tions with 
topological charge k is independent of temperature /3 _1 , S = 87r (see 2.2.1), no 



explicit (3 dependence is allowed in the definition of (j) a . Finally, if the scalar gauge 
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multiplet is thermalized with the sector of trivial topological charge k = 0, it must 
carry charge under the gauge group SU(2) (that is, transform nontrivially under 
gauge rotations). Then the index a is just a SU(2) valued color index. 

4. It can be shown [23] that there exists a unique construction for a family of phases 
{0 a } living in the adjoint representation satisfying the conditions described in [3] 
This construction is necessarily nonlocal, and only admits contributions from \k\ = 
1 (Harrington-Shepard) (anti)calorons of trivial holonomy. Non-trivial holonomy 
calorons can be apriori excluded from contributing to (f) a since such configurations are 
known to be unstable under trivial-topology fluctuations |9] (The average Polyakov 
line is zero, (trL) = [8j). Configurations of topological charge \k\ > 1 are also 
excluded, by Lorentz symmetry arguments [23j so that <p a emerges from spatial coarse 
graining over only trivial holonomy calorons of topological charge \k\ = 1. It turns 
out that the set {4> a } defines the kernel of - and completely determines a - second 
order differential operator D. The fact that the scalar field (j) a both has vanishing 
Euclidean energy density and belongs to the kernel of T> can be exploited to fix the 
gauge invariant effective potential V ((f)) uniquely for (j) a . This action also introduces 
a mass scale as an integration constant, the Yang- Mills scale A. 

5. Since no object other than the inert adjoint scalar field 4> a may emerge from the 
sector of topologically nontrivial fundamental configurations, the coupling of (fi a with 
the sector of zero topological charge must be minimal. That is, it should be realized 
via a covariant derivative = d^ + a^. Higher dimensional operators are forbidden 
by perturbative renormalizability [25]. 

Thus, we arrive at a scalar field <j> and an effective Lagrangian £ eff describing inter- 
actions with the sector of topologically trivial charge k = 0. The effective theory described 
by the scalar field <ft and the Lagrangian £ cS is discussed in some detail in the next sections. 



3.2 SU(2) thermal ground state in the deconfining phase 
3.2.1 Inert adjoint scalar field 

As described in |3.1.1| upon spatial coarse graining over the caloron sector of the funda- 
mental SU(2) theory a Lorentz invariant multiplet (j) a is generated which describes ground 
state physics up to a phase transition temperature given by the Yang-Mills scale of the 
theory. The modulus of 4> a does not depend on space or time, and is given by 



^tr02 (3.1) 



A 3 /3 
~2n' 



(3.2) 



with the inverse temperature (3 = A. The equations of motion ruling over the dynamics of 
(f) a are derived from an (effective) Euclidean Lagrangian of the form 



C<b = tr 



(d T <j>f + V [cj) 2 )] , (3.3) 
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ffl (in 



Figure 3.1: A spatially homogeneous and isotropic static field 4> a is generated by spatial 
coarse-graining over interacting calorons of trivial holonomy and topological charge \k\ = 1, 
shown in (i) as fluctuations inside a spherical integration volume. The inert field <p a and the 
pure-gauge configuration which solves the effective Yang-Mills equations together define 
a thermal ground state, (ii) depicts domainization occurring from interactions between 
calorons. Domain walls intersections are identified as magnetic monopoles. Domainization 
effects can be taken into account by loop-expanding the quantity under investigation about 
the thermal ground state in the effective theory. Figure taken from |32) . 



the potential V is given by the expression 

V(<f,)=tr^, (3.4) 



with the inverse of the scalar field defined as <p 1 = j^m. The mass scale A emerges as 
an integration constant and serves as the Yang-Mills scale of the theory, defining a phase 
transition temperature T c = X c ^ which bounds the deconfining phase from below. The 
critical dimensionless temperature A c is defined by the evolution of the gauge coupling e(A) 
in the effective theory (see section 3.2.2[ ). At A c , the gauge coupling exhibits a logarithmic 



singularity so that e(A) = -4.591og(A - A c ) + 18.42. Numerically, A c = 13.87 [3T]. The 
transition temperature T c is the only parameter of the theory. To make contact with exper- 
iment, this parameter must be fixed so that quantitative calculations may be performed. 
This will be done in section 13.41 

Up to global gauge transformations, a solution to the Euler-Lagrange equations of 



motion determined by (3.3) and (3.4) is given by 



_^ tl exp(±— i 3 r). (3.5) 

The field <j) represents a spatially homogeneous background for the dynamics of coarse- 
grained, propagating, trivial-topology gauge fields. It will also break the gauge symmetry 
SU(2) down to its Abelian subgroup U(l). The effect of this symmetry breaking will be to 
confer a temperature-dependent mass to two of the three available vector modes. This is 
discussed in section [3.3.11 
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3.2.2 Effective action and bare thermal ground state 



In section 3.1.1, it was shown how to construct a scalar field resulting from spatial coarse- 



graining over Harrington-Shepard calorons. The procedure integrates out the topologically 
nontrivial sector of the fundamental theory. The dynamics of the topologically trivial 
remainder is then determined by an effective theory valid up to a scale \<p\, characterized 
by an effective Lagrangian £ cB describing the interaction of the scalar field (f) with the 
topologically trivial sector (k = 0). It has been shown [25] that £ cff can be uniquely 
defined - up to a maximal energy resolution \<f>\ - and is given by 



C,, H [a. 



tr 



+ 



A 6 



(3.6) 



where = a" t a is a propagating, trivial-topology gauge field in the adjoint representation 
remaining after coarse-graining. The effective field strength is defined as usual 



G 



ze[a M , a v \ 



(3.7) 



with e the effective gauge coupling. The effective coupling constant e is not the same as 
the coupling constant of the fundamental theory. The evolution of this coupling constant 
as function of dimensionless temperature e(A) follows from demanding thermodynamical 
consistency |25j . As already mentioned, the coupling behaves logarithmically at the critical 
temperature A c , e(A)|^ c ~ — log (A — A c ), and the theory undergoes a second order phase 
transition to the preconfining phase. We remark here that, for temperatures higher than 
A c , local gauge transformations can always be found such that the Polyakov loop Pol[a^ g ] 
is transformed from Pol = — 12 to Pol = 12- This reflects the electric Z2 degeneracy of 
the thermal ground state estimate, and the theory for A > A c is thus deconfining. Further 
details can be found in the literature [25, 29j. The covariant derivative imposes minimal 
coupling between the effective k = gauge field and the \k\ = 1 coarse-grained inert scalar 
field, and is written as 



D„ 



te a 



(3.8) 



Note how the effective action in (3.6) at resolution \(f>\ splits into a classical Yang-Mills 
term describing the propagation and interactions among topologically trivial fluctuations, 
an interaction term between k = and \k\ = 1 configurations given by the square of the 
covariant derivative acting on tfi, (D^cp) 2 (Higgs mechanism for two out of three propagating 
directions in the SU(2) algebra), and a potential V for the field (f> [25J. 

The Euler-Lagrange equations governing the effective theory can be read from the 



Lagrangian (3.6) 



ie[4>, D„(, 



(3.9) 



The inert scalar ground state configuration (3.5) and the pure gauge configuration 



. 2ir 

ep 



(3.10) 



are ground state solutions of (3.9) since G^la 1 ^] = D u [a^](f> = 0. Apart from (small) 



radiative corrections and up to global gauge transformations, (3.5) and (3.10) constitute a 
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good first estimate of the thermal ground state of the effective theory in the deconfining 
phase. Accordingly, a good first estimate of the energy density of the thermal ground state 
is given by 

A 6 

p gs = C cft [a s ;\ = tr = 4vrA 3 T , (3.11) 

which is interpreted as a temperature dependent cosmological constant. Together, <j) and 
a® 5 represent an a priori estimate of the bare thermal ground state which turns out to be 
quite accurate. 



For later use, we write (3.11) in SI units. The temperature T is rescaled as T — > ksT 1 



the mass scale A as A — > ^, and the energy density of the thermal ground in the deconfining 
phase depends on the temperature through the expression 

pSS=47r (^) kBT ' (3 ' 12) 

This energy density will later serve as a criterion for thermodynamical consistency in the 
presence of external fields. 



3.3 Thermal quasiparticles and Feynman rules 



3.3.1 Mass spectrum 

In the effective theory of deconfining Yang-Mills thermodynamics, each broken generator 
t a of the original gauge symmetry, causes the original massless mode associated with t a to 
acquire a mass m a |29j 

(3.13) 



-2e 2 tr 



To interpret (3.13) physically, we impose unitary gauge for the thermal ground state 

= 2 |0| t 3 , (3.14a) 
a es = 0. (3.14b) 



In this gauge we have 

m = mi = m,2 = 2e(A)|0| = 2A 



e(A) 



2Ae 



2vr. 



-T 



A 



y/X \ A ' J V -2kT 

7773 = . 

In the fundamental theory, all gauge modes are massless, while in the effective theory two 



(3.15a) 
(3.15b) 



is 



of the three gauge modes acquire mass at tree-level. Note that the mass m in (3.15a 
temperature-dependent. Recall that, at the critical temperature A c = jT Cl the effective 
coupling e(A) shows a logarithmic singularity e(A)|^ c ~ — log(A — A c ) for A \ A c . This 
signals a phase transition to the preconfining regime. Near this critical temperature, the 



mass m in (3.15a) diverges, and interactions between massless and massive modes decouple. 



That is, radiative corrections to ( 3.15b[ ) becomes exactly zero. 

The unbroken gauge subgroup U(l) is responsible for the remaining massless mode after 
coarse-graining. It is this remaining symmetry, and the vanishing of radiative corrections 



to (3.15b) at T c , that will lead us to speculate about the fundamental origin of the Standard 



Model photon (see section 3.4). 



But first, we review the Feynman rules in the effective theory. 
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7A b 




Figure 3.2: 3- and 4- vertices. Taken from |59j 



3.3.2 Finite temperature Feynman rules in unitary gauge 

In unitary-Coulomb gauge, which is a physical, completely fixed gauge, the real-time prop- 
agators of the fields are obtained by a Wick rotation to real timd^of the usual Matsuraba 
imaginary time propagators |35[ |4"T| |4"2] . In the real-time formulation of finite-temperature 
field theory [ 4~T] it possible to discern quantum from thermal fluctuations. To tell the for- 
mer from the latter is necessary since quantum fluctuations away from the mass-shell will 

1 2 

be limited by the physical resolution \<j>\ in the unitary-Coulomb gauge. Also, the fields 
may only propagate thermally [25]. The Feynman rules for the SU(2) vertices are taken 
from the literature as well, by substituting the fundamental gauge constant g — > e. 

To conform with the nomenclature of |59j we refer to the fluctuative SU(2) directions 
1 and 2 as tree-level heavy (TLH) and to direction 3 as tree- level massless (TLM). In 



unitary-Coulomb gauge (3.14) the real-time bare thermal propagators in momentum space 
are then given by [59j 



TLH,0, x 



D 

n TLM,0/ \ 
D ab,»l (P) 



-8 ab D^ u [2tt5( P 2 - m 2 )n B (\ PQ \/T)] 



'hib 



P, 



p* 



+ 2ir5(p 2 )n B (\p \/T) 



P 2 



where ub(x) = l/(e x — 1) denotes the Bose-Einstein distribution function and 



(3.16) 
(3.17) 

(3.18) 



The four-vector u„ = (1, 0, 0, 0) specifies the restframe of the thermal bath and the 3-space 
transverse projection operator Ptt u is defined as 



P°° = p§i = if = o . 

p« = ^'-py7p 2 . 



(3.19) 
(3.20) 



TLM modes carry a color index 3 while TLH modes have a color index 1 and 2. The three- 

lr This method becomes problematic for self interactions and loop calculations of higher order than two 
It has been shown that such difficulties do not apply to our approach to SU(2) thermodynamics 
|26| . Note also that a change of signature from Euclidean to Minkowskian is trivial as far as the physics 
of the ground state estimate of section |3.2.2| and the associated dynamical gauge-symmetry breaking are 
concerned: neither ground state pressure, energy density or quasiparticle mass depend on Euclidean time. 
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and four-gauge-boson vertices are given as (see figure 3.2) 

T mabM 1, k) = e(27r) 4 5(p + q + k)e abc 

x \sT{q ~ PY + 9 vp ik - q)» + g">(p - k) u ] (3.21) 
r WaT>^, S ,r) = -ie\2n)*5(p + q + s + r)l 

eabee C de(g^g va -g^ a g up ) 

+ e ace e bde WgP° -g^g vp ) 

+ e ade e bce (9 flV 9 pa -9 fip 9 Va )}- (3-22) 

Finally, recall that, before calculation, each loop-diagram has to be divided by i and the 
number of its vertices (symmetry factor) |41] , 

The maximal off-shellness of effective gauge modes as well as the momentum transfer 
in four- vertices are constrained by coarse-graining down to a resolution \<f>\ [25[[28]. TLM 



modes (3.17) may never be off-shell beyond 

\p 2 \<\4>\ 2 - (3.23) 
TLH modes on the other hand may only fluctuate thermally, that is, on-shell. Condition 



(3.23) also fixes the momentum transfer in a three- vertex by momentum conservation. 
Momentum transfers in an effective four vertex cannot be larger than \(j)\ 2 |25| . 

We now have all the information on hand that we need to calculate in the effective 
thermal theory. The only step left is to fix the scale A so that we may get quantitative 
results. This is achieved by identifying the Standard Model U(l) propagating photon with 



the massless SU(2) mode in the effective theory. This is discussed in section 3.4 

Afterwards, we introduce the polarization tensor for the TLM mode, and proceed to 
calculate its influence on photon propagation. This is the subject of section [3~5| 



3.4 The SU(2) CMB hypothesis and some implications 

The nonperturbative approach to the thermodynamics of the SU(2) Yang-Mills gas re- 
viewed in the preceding sections can now be summarized as follows: An adjoint scalar field 
(j) a is defined by coarse-graining over trivial-holonomy field configurations with topological 
charge \k\ = 1. A thermal ground state composed of (f) a and the pure-gauge field configura- 
tion a^ s , considered together with topologically trivial (|A:| = 0) gauge field fluctuations, is 
seen to effectively describe the fundamental SU(2) Yang-Mills theory at finite-temperature. 
In the effective theory, two of the three fundamental gauge modes of the local SU(2) sym- 
metry become massive due to <j) a by the Higgs mechanismus. The remaining U(l) symmetry 
is responsible for the masslessnes of the third gauge mode. 

At this point, it is worthwhile to consider whether this symmetry is realized in nature. 
An intriguing possibility is the identification with the U(l)y gauge group of the Standard 
Model. As argued in |17| 124"! 127] [30l [58], this would immediatly fix the critical temperature 
T c with the cosmic microwave background temperature T CMB = 2.725 K. It follows for the 
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SU(2) Yang-Mills scale 

A = A CMB = 2vr *f*£ = 2vr (3.24) 

= 1.70433 • 10~ 23 J 
= 1.06376 • 10" 4 eV, 

with A c = 13.87. Once we have fixed the Yang-Mills scale A = A CMB we can proceed to 
calculate the polarization tensor of the TLM mode. This is done in the next section. 

3.5 Polarization tensor and photon propagation 

Having identified the massless mode of the coarse-grained effective theory with the familiar 
U(l) photon, we now proceed to investigate what effect the existence of a thermal ground 
state may have in the physics of photon propagation. We expect the thermal ground state 
to interact with the massless mode via the massive vectors. This interaction is encoded 
in the polarization tensor U^ u (pQ,p), which we compute perturbatively at one-loop level. 
The rest of the chapter is dedicated to the main results of this computation. 

3.5.1 Prerequisites 

The dispersion relations for the massless gauge modes are given by the poles of the full 
thermal propagator D^fj^ik). In general, any dressed thermal propagator T>^ v can be 
written as the geometric series 

V^ = V% + V ^U^V% + ... , (3.25) 

where V^ v is the bare thermal propagator and the polarization tensor IP^ is defined as 
the sum of one-particle irreducible (1PI) diagrams with 2 external lines [53]. The series 



(3.25) can be written as 

V^ = V% + Vl a WPV pu . (3.26) 

tensor 



Defining the inverse propagator as T>^T> al ' = ', we get an equation for the polarization 



Tl^ = V-l-V\l. (3.27) 



Further, is transverse 

= . (3.28) 

Here, p^ is the 4-momentum vector of the thermal excitation. By separating the space-time 
transverse projector operator 

P^H^-pV/p 2 , (3.29) 

in the space transverse projector Pt and space longitudinal projector Pl 

= pM = pfi = , 

P^ = 5 ij - pY /p 2 , (3.30) 
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we can expand the space-time transverse polarization tensor as 



(3.31) 



with scalar functions G and F. Before analytic continuation from imaginary time, the 
full interacting propagator D^^f,(k) for the massless mode of the effective theory can 
be written in the Coulomb gauge in the Euclidean formulation as the linear combination 
[351 E9] 



D 



TLM 
ab,fj,u 



(P) 



P 



1 



+ 



P 



1 



(3.32) 



^G-p 1 p 2 F-p 2 

where = (1,0,0,0) specifies the restframe of the many-body system. Note that, as 
expected, while the first term in the right hand side of expression (3.32 ) describes transverse 
propagation, P£ v ~ (G — p 2 ) -1 , the second term describes longitudinal propagation 
P£ U D^ M ~ (F — p 2 )~ l - The corresponding dispersion relations are then determined for 
all modes by the poles of (3.32). For the transverse modes, this results in 



Po = P 2 + G(p ,p) . 



(3.33) 



Separate pq = oo — ry, and assume weak damping 7/w <C 1. The transverse dispersion 
relation is then given by 



oo 2 (p) = Rep 2 = P 2 + ReG(o;(p),p) 



71PJ 



1 



2w(p 

Similarly, for the longitudinal mode 
w£(p) = p 2 +ReFK(p),p), 



ImG( W (p),p). 



7l(p) 



1 



-ImF(wi(p),p) , 



(3.34) 
(3.35) 

(3.36) 
(3.37) 



2w L (p) 

where we added the subscript L to differentiate from the transverse relations. Taking p 
parallel to the z-axis, we can extract from (3.31 ) expressions for the scalar functions G, F 
as components of H^ u 



G(w(p),p) = nn(w(p),p) = n 22 (w(p),p), 

„2 



%(p),p) 



p 



n 00 . 



(3.38) 
(3.39) 



3.5.2 Polarization tensor at the one-loop level 

In order to calculate the polarization tensor one has to sum all 1PI diagrams to the required 



precision. Figure 3.3 shows the diagrams contributing to the photon polarization tensor 
IP' - ' at the one-loop level as a function of incoming momentum p. The total polarization 
tensor is just the sum of the two diagrams 

ip» = n7( P ) + ri£>). (3.40) 
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fi, a = 3 




v, b = 3 



fi, a = 3 




u.b = 3 



Figure 3.3: The photon is represented by single- line propagators, double- line propagators 
correspond to massive modes. Taken from |59j . 



Diagram A can be written by following the Feynman rules described in section 3.3.2 which 
results in 



i 



d 4 k 



2i J (2vr) 4 



e' e ace [g^(-p - k) x + g^ x (k -p + k)" + g^(p - k + pf) 



X (Scd) 9 pa 



kpk a 



2it5{k 2 -m 2 )n B (\kQ\/T) 



k 2 — m 2 

x e dbf [g™(-k - pY + g™(p + p - kf + g K °(-p + k + kf] 
(p - k) x (p - k) h 



x (Sef) 9\k 



(p — k) 2 — m/ 



(p — k) 2 

+ 2tt5(( P - kf - m 2 ) n B (\p - k \/T) 

(3.41 



Here m = 2e\(f)\ = 2ey / A 3 /27rT is the induced mass of the two other gauge modes in the 
effective theory (see ( 3.15[ )), and n B {x) is the Bose- Einstein function. 
Diagram B is associated with the following expression 



(5^(-<W (9P* ~ ^f) ~ m 2 )n B (\k \/T)] 

x (-ie 2 )[e abe e cde (g^g^-g^g^) 



+ e ace e bde (g^gr° - g^g^) + e ade e bce (g^ g^ - g^g™)] . 

(3.42) 

As shown in [46], G and F may be calculated selfconsistently by numerical methods from 
this expressions for the polarization tensor. This is described in the following sections. 



3.5.3 Transversal modes 

Transversal propagation is determined by the scalar function G(lj, p), which can be cal- 
culated from the polarization tensor by use of equation ( 3.38[ ). In both (3.41) and (3.42), 
the off-shellness of the k integration is restricted in the effective theory by the maximal 



resolution \4>\ after coarse-graining, see (3.23). This results in the support condition 



{p + kf 



\G + 2p- fc + 4e 2 |^| 2 | < |</>|' 



A 3 

2irT ' 



(3.43) 
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where we have used p 2 = G and k 2 = Ae 2 \(j)\ 2 . The k integration in (3.41) and (3.42) is to 



be taken over all values of k such that (3.43) is fullfilled. It can be shown that diagram 



A, which is given by equation (3.41), is identically zero over the support (3.43) [46J. Since 



(3.41 ) is the only imaginary contribution at the one-loop level, there is no finite transverse 



photon width (3.35). 



The scalar function G(u, p) can then be calculated up to one- loop from diagram B, that 



is equation (3.42), by integrating over the support determined by the inequality (3.43). To 



evaluate (3.42), transform the integral to dimensionless variables 



£ = 



h 
\<f>\ 



l 



p = 



k 2 + k\ , 



arctan 



k 2 



(3.44) 



Equation (3.38) gives G when p is taken parallel to the z-axis. Recall the definition of 



dimensionless temperature 
2vrT 

A 



A 



(3.45) 



and obtain an integral expression for the radiatively induced screening function G, nor- 
malized to G -> G/T 2 



£ _ nil _ 

rp2 ~ ^2 — 2~'2 



J d£ J dp e(A) 2 A- 3 (-A + 



d£ J dp h G (£,p,A) 



4e(A)^ 



n B ( 2tt\- 3 / 2 y/p 2 + ( 2 + 4e(A)< 



+ e + 4e(A) 2 



(3.46) 
(3.47) 



where we have noted the A-dependence of the coupling e by writing it explicitly. The 



integration must be done over the support (3.43), which we also write in dimensionless 

(3.48) 



variables. Define the dimensionless momentum X 

IpI N 



X = 

T T ' 
and the support function s± 



s± (€,pA, x ,g) 



0^ + ^ (± v / ^+5Vp 2 + e 2 + 4e2 - ex) + 4e 5 



Then the support condition (3.43) can be written in dimensionless form as [46 

s± (£,p,\x,g) < i , 



(3.49) 



(3.50) 



with g = G/T 2 . The doubling of inequality (3.43) comes from writing the delta function 
in (3.42) in terms of A: |59j 



5{k 2 



m 



5 (k° + ^\k\ 2 + m 2 ) +5(k°- y/\k\ 2 + m 2 



(3.51) 
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Define the function Hq{X,X,q) as the integral of Kq (£, p, A) summed over the support 
defined by s± (£, p, A, X, < 1 

/oo /*oo 
d£ / dp0(l-s a ({;,p,\X,g)) h G (£,p,\) . (3.52) 
-00 Jo 



where 9{x) is the Heaviside step function. Equation (3.46) then becomes 



g- H G (X,X,g) = 0. 



(3.53) 



It turns out that (3.53) is enough to uniquely determine G (46]. For a given temperature T 
and momentum p = |p| z, the integral (3.52) can be calculated numerically as a function 
of A, X and g. The value of g consistent with its defining equation (3.47) can then 



be determined by searching for the root of expression (3.53) using standard numerical 



methods. This root g of (3.53) happens to be unique for each momentum and temperature. 



This is all done using the Mathematica numerical software. Our algorithmic procedure is 
then as follows: For a given fixed temperature T and various values of X within an interval 
[0, X max ], the value of g consistent with (3.47) is calulated by finding the root of (3.53). 



Each point is then used to interpolate g(X,T) in X over [0,X max ]. 

In figure 3.4a we have plotted the resulting function G(X,T)/T 2 for different tem- 
peratures. Note the existence of both an screening region where G(X, T) > and an 
antiscreening region where G(X, T) < 0. Asymptotically, G goes to zero, its maxima and 
minima also tend to zero with increasing T. 



G 



With G nonzero we expect the dispersion law (3.34) to deviate from the massless case 
= 0. Defining the dimensionless frequency 



Y 



T 



we write the transverse dispersion law (3.34) as 



Y 



X 2 + 



G(X,T) 
T 2 



(3.54) 



(3.55) 



This dispersion relation for the transverse modes is plotted in figure 3.4b for different 
temperatures. The cutoff frequency uj* is defined as 



uj*(T) = ^fG{X = 0,T) = Y(X = 0,T)T = Y*(T)T. 



(3.56) 



As seen in figure |3~4b"{ the cutoff frequency is the minimal frequency available for transverse 
propagation at a given temperature. From the same figure is also clear that Y* decreases 
with increasing temperature. The functional relation Y*(T) is very important for our 
experimental considerations, since it defines the upper boundary of the screening region 
Y < Y* where photon propagation is forbidden. The experimental consequences of such 
screening effects will be investigated in chapter [5] In the next section, an expression for 
Y* (T) will be derived from curve fitting. 



3.5.4 Some properties of the photon screening function G 

In this section we give numerical fits for three characteristic points of the screening function. 
We work in SI units, so as to build an intuition on experimental orders of magnitude. 
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(a) The dimensionless screening function G/T 2 calculated for a given temperature T and 
momentum X — pa/T by solving for the unique root of equation ( |3.53[ ). To obtain a 
complete plot in a given interval of X, the values of G for a given T were successively 
calculated for different values of X partitioning the interval. The curves shown correspond 
to the temperatures T = 5K (red), T = 8K (blue) and T = 10K (green). 
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(b) The dimensionless dispersion law Y(X) = 
the value of the cutoff frequency Y* = Y(X 
plotted temperatures are T = 5K (red), T = 
dispersion law Y — X is shown as a black line 



^X 2 + G/T 2 calculated from G. Note that 
= 0) for each temperature falls with T. The 
8K (blue) and T = 10K (green). The U(l) 



Figure 3.4: Transverse screening mass G/T 2 (a) and dispersion relations (b). 
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In section 3.5.3 we denned the cutoff frequency as the minimal available frequency for 
any given temperature, hu*(T) = c\J G(|p| = 0, T). We calculated v* for temperatures 
between 2 and 165 Kelvin. The results are depicted in figure 3.5 Figure 3.5b shows how 
v* grows with decreasing temperature, reaching a maximum at about 5.4 Kelvin of about 
v* = 17GHz, then decreases and suddenly drops to zero below 4.5K. To calculate the fit, 
we have assumed a power law functional relation of the form 



v*(T) = c + bT* . 



(3.57) 



We excluded from the fit data points below 7K, as seen in figure |3T5b| The Mathematica 
ire was 

v*{T) 



software was used to fit the data points to (3.57), resulting in the expression 

-0.53 



GHz 



42.71 



+ 0.22. 



(3.58) 



How well this fit matches the data is shown in figure 3.5 



Another interesting characteristic frequency is the crossover frequency v c , defined as 

G(u e ,T) = 0. (3.59) 

This is the frequency at which the screening function G goes from screening to antiscreen- 
ing. Just as v*, v c drops to zero below 4.5 Kelvin. A fit of the data points to a power law 
of the form (3.57) gives as result 

1.12 



GHz 



1.83 



+ 13.48. 



(3.60) 



This is plotted in figure 3.6 



In chapter [5] we introduce the radiance function of a black emitter. It will be shown 
that the ratio between an SU(2) and an U(l) emitter is given by the term 



c 2 G 

(H 2 



1 (*-"•) • 



(3.61) 



The maximal difference between U(l) and SU(2) radiance in the antiscreening regime, 
where G is negative, is then given at the frequency vm where G/v 2 reaches a minimum 



G(v M ,T) G(u,T) 



(3.62) 



Thus vm is the frequency of maximal antiscreening. A fit of the data to a power law of 
gives 



the form 3.57 



vm{T) 
GHz 



3.45 



1.08 



+ 12.90. 



(3.63) 



Figure 3.6 shows the calculated points overlaid with the fitted curve. 
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(a) The set of cutoff frequencies {v*} calculated from a discrete set of temperatures {Ti} between 2 and 
165 Kelvin, shown here as black dots, overlapped with the curve v *(T) = 0.22 + 42.71 • T~ ' 53 fitted from 
the dataset {v*,Ti\. 



„ 10- 
X 

o 



T[K] 

(b) Closeup of the dataset {u* , Ti} near the tem- 
perature T = 4.5K for which the cutoff frequency 
v* rapidly falls to zero. The temperature region 
between 2K and 7K was not included in the cal- 
culation of best-fit parameters. 




T[K] 

(c) Closeup of the dataset {v* ,Ti} near T = 7K 
which was used as the lower bound of the curve 
fit. 



Figure 3.5: The cutoff frequency u* was fitted to a power function in T, using a dataset 
{v*,Ti} calculated from discrete set of temperatures {Tj} between 2 and 165 Kelvin. The 
original set had i = 1, .., 306 data-points. The first 50 of those, corresponding to tempera- 
tures between 2K and 7K where ignored in the calculation of best-fit parameters. Taking 
v* in GHz and T in Kelvin, the best fit for the remaining 256 points in the dataset is given 
by the power law u*(T) = 0.22 + 42. 7T" - 53 . 




3.5.5 Longitudinal mode 

Consider now the propagation of longitudinal modes, characterized by the poles of the 
term oc UnU u in (|3~32j). On shell, p 2 = F, and (|3~39j) simplifies to 



X 2 EE P 



J"2 



noo 

J>2 ' 



(3.64) 



To calculate this expression numerically, we first assume that diagram B of figure 3.42 



is 



the only contribution to the polarization tensor. Afterwards, we check that diagram A 
vanishes on-shell F = p 2 . The expression for II 00 



3.3 



in figure 
equation (j3.42h 



Il'g can be read from 



ttOO 

J>2 



dp 2e 2 \- 3 



3 + 



4e 2 



Vp 2 + e 2 + 4e 2 



d£ /" dp (£,p,A) 



(3.65) 
(3.66) 



The integrations are subject, as in the transverse case, to the following constraint 



s±(€,p,\,X,f) < 1, 



(3.67) 
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with the support function s± defined as 
s ± fop,\,XJ) = 



f\3 \3/2 



(2tt)S 



7T 



(3.68) 



with / = F/T . As we did in section 3.5.3 for the transverse case, we express the inequality 



(3.50) by using the Heaviside step function 6(x), and define a function Hp (X,X, f) as the 



integral of hp (£, p, A) in (|3.66 ) summed over the support determined by (3.68) 



/oo roo 
d£ dpe(l-s a (£,p,\,XJ)) hp(£,p,\) . (3.69) 
,. _._ _ -co J 



Then the consistency condition for the integral equation (3.65) is just the expression 
X 2 - H F (X,X,f) = 0. 



(3.70) 



Using the same method as in section 3.5.3 which is to search for the root of equation |3.70| 
in /, given fixed values of X and A, results in numerical difficulties, due to the spectral 
function F(X) being double valued. The longitudinal screening function F, in contrast to 
G, separates in different branches. Therefore, the searching algorithm has to be modified 



to look for the roots of (3.70) in X, with fixed values of / and T. This root turns out to be 
unique. The search results are displayed in figure 3.7a, where we have plotted F(X)/T 2 
for different temperatures. Note the existence of two branches of F. Furthermore, X is 
confined to a finite domain between X = and X = 0.4, as shown in the plot. This strange 



behaviour is reflected in the longitudinal dispersion law of figure 3.7b 



Using similar numerical methods as those described above, we have checked that dia- 

a2 



gram A, given by expression (3.41 ), vanishes on shell F = p 



In dimensionless quantities the longitudinal dispersion law reads 



T 



X 2 + 



F(X,T) 

rp2 



(3.71) 



This relation is plotted in figure 3.7b for three different temperatures. Notice the ex- 
istence of various branches, consequence of the double- valuedness of F(X). The lower 
branch shows a dispersion law with negative slope, reaching its maximum available spatial 
momentum at the point where Y{X) = X. Further, the upper branch of Yl becomes su- 
perluminal at very low X, or large wavelength X -1 , and diverges to infinity very quickly. 
This suggests that the longitudinal modes are unphysical. 

Consider also the magnetic-electric dual interpretation for the thermal SU(2) Yang- 
Mills theory described in |25], where the propagating electric field in the SU(2) theory 
becomes a magnetic field in the thermal photon plasma. Recall that longitudinal oscil- 
lations describe electric charge density waves |35j . In the dual interpretation, the long 
wavelength branching modes in figure |3.7b become magnetic density waves, and cannot 
transfer energy in any interaction with electric charges. We gather that no thermalization 
can occur for these modes. 

From the above arguments we can only conclude that the longitudinal modes are, if not 
unphysical, at least unable to propagate nor transport energy. Thus, we may ignore the 
longitudinal modes in our further discussion. The only excitations of the SU(2) thermal 
photon gas we will consider are the transverse modes of the tree-level massless gauge field. 
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(a) The dimensionless longitudinal screening function F/T 2 . As seen in the plot, F(X) is double 
valued, so the searching algorithm used for the transverse modes had to be modified to look for 
the roots of ( |3.70| l in X, with fixed values of F/T 2 and T. This root turns out to be unique. The 
curves shown correspond to the temperatures T = 5K (red), T = 8K (blue) and T = 10K (green). 
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(b) The longitudinal dispersion law Yl(X) 
branches, consequence of the double- valuedness of F 
T — 8K (blue) and T = 10K (green). The U(l) dispersion law Y 



2 + F/T 2 calculated from F. Note the two 
The plotted temperatures are T = 5K (red), 



X is shown as a black line. 



Figure 3.7: Longitudinal screening mass F/T 2 (a) and dispersion relations (b). 
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3.6 Resummation of one-loop radiative corrections 

Apart from the effect on photon propagation, the transverse screening function G may 
be used to calculate thermodynamical properties of the SU(2) gas. Trough the procedure 
of resummation, corrections to the energy density and pressure induced by interaction 
of SU(2) photons with the thermal ground state can be extracted. We will find such 
correc-tions to be linear in the temperature. The thermodynamical consistency of this 
resummation approach is also discussed in this section. 



3.6.1 Prerequisites 



Recall from section 3.2 that the ground state of the thermal SU(2) Yang-Mills gas is 
composed of a scalar field <fi (3.5) and a pure-gauge configuration (3.10). From the effective 



action (3.6) the energy-momentum tensor 6^ v of this configuration may be calculated 



' [ii/ 



-4-irATg 



[IV • 



(3.72) 



Take the SU(2) gas to be a perfect fluid in thermal equilibrium, 9& = (P gs + p gs )u fl u l/ + 



P es 9[tu, with the restframe specified by the 4- velocity vector u^. Then we read from ( |3.72 ) 
that p gs = — P gs = 47rAT. As given in (3.11), the (classical) ground state configuration has 
a nontrivial energy density which grows linearly with temperature. For later comparison 
with higher-order corrections, we write this result in the form 



rp4 



2(2tt) 4 A" 3 ~ 3117.09 A" 



(3.73) 



where as before the dimensionless temperature is defined as A 



2ttT 
A ■ 



Note also that the ground state pressure P gs in (3.73) is negative. This may be in- 



terpreted microscopically as magnetic monopole-antimonopole |37[ l38l H31 |5T] attractive 
interactions due to the dominance in thermal fluctuations of small-holonomy caloron con- 
figurations [9]. 

Equation (3.73) gives the tree-level contribution to the gas thermodynamics. What 
about higher order corrections? As shown in |18J . linear contributions AP gs,1_loop and 
^^gs,i-ioop ^ Q p gs an( j ^ gs ^ reS p ec tively, arise from the free fluctuations of tree-level massive 
modes at high temperature. This correction is given by 



Ap g 



,1 — loop 



AP g 



, 1 — loop 



rp4 



-2(2ir) 4 \- 



-3117.09 4 A" 



(3.74) 



2e|, 



8 ^% i s the temperature scaled mass 



where a — T — A3/2 

Note that, since both P gs and AP gs ' 1_loop are linear in temperature, a Legendre trans- 
formation to the energy density p = — P must vanish. But p ss + A/9 gM ~ loop = also 
vanishes, therefore both corrections are thermodynamically consistent. 

In the following, we show how a correction to the base energy and pressure of the 
ground state can arise by assuming selfconsistent propagation at the 1-loop level of tree- 
level massless modes. This correction will depend linearly on temperature, again signaling 
ground state effects on photon propagation, and vice versa. 
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Figure 3.8: The Dyson series for the dressed propagator of the massless mode, containing 
only tadpole diagrams, since it is the only nonvanishing diagramm on the induced mass 
shell, see section 3.5 A thin line refers to the propagation of the tree-level massless, a 



thick line to the propagation of tree-level massive modes. Taken from [12] . 



3.6.2 Resummation 



A resummed perturbation expansion is given by replacing the bare propagator with the 
dressed one. In figure 3.8 the diagrammatic expansion of the dressed propagator up to 

as given by |3.41[ vanishes in the effective 



one-loop is shown. Since diagram A in figure 3.3 



theory (see section 3.5 ), we only need to consider the one-loop tadpole (diagram B in figure 



3.3). We can then calculate resummed thermodynamical quantities by the usual method 



of expanding In Z in bubble diagrams, but replacing the bare propagator with its resummed 
(dressed) version. 

But we already calculated the dressed one-loop propagator (determined by the one-loop 
polarization tensor). As we reviewed in 3.5.1 the dressed propagator in figure 3.8 as a 



result of photon interactions with the ground state, displays a modified dispersion law for 



the propagating, transverse massless modes. From section 3.5.3 we recall the dimensionless 
dispersion relation 



Y 



X 2 + 



G(X,T) 

J^2 



(3.55) 



and the definitions Y = X = -r^. Thus, we can insert the modified dispersion relation 



T ' 



(3.55) in the thermodynamical integrals for energy density and pressure 
P 



1 



7T^ 



dpp 



oJ3u _ I 



P = J d PP 2 m(l-e~ 



(3.75) 
(3.76) 



derived from the partition function for the noninteracting (zero chemical potential) massless 
boson with two degrees of freedom |35j 



-2V 



d 3 p 

(27p 



1 



: /3uj + lnl - e 



(3.77) 



Then, we get the 1-loop resummed energy density and pressure 

•jl — loop, res 



P 



^4 

pi — loop, res 



VX 2 + G/T 2 



1 f 

vr 2 Jo p ^xz+G/T2 



1 



1 f°° 

—s / dXX 2 \n(l -e 
n Jo 



), 



(3.78) 
(3.79) 
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where we have rescaled by T 4 to make both expressions dimensionless. Define the dimen- 
sionless corrections AP and Ap to the contributions of free, massless particles as 

q! — loop, res q\ — loop, res ^-2 

A P^ r4 =^~"15' (3 ' 80) 

^^^Dl — loop, res pi — loop, res ^-2 

AP= =^ = — ^ _ (3.81) 



rp4 



rp4 



Corrections (3.80) and (3.81) can be calculated, by first calculating G as shown in 
section 3.5.3 as a function of temperature. Figure 3.9 depicts both AP and Ap as functions 



of the dimensionless temperature A in the high-temperature region 100 < A < 500 (recall 
that A c = 13.87). The calculated data points can be fitted to a curve assuming power law 
dependency 



Ap 
AP 



c p A 7 
cp\ 5 



The best fit is given for the parameter values 



c p = 3.9577, 
c P = 8.4963, 



7 
6 



-3.02436, 
-3.00904. 



(3.82) 
(3.83) 

(3.84) 
(3.85) 



First note that the resummed corrections are linear, similar to (3.73), but three orders 



of magnitude smaller. Equations (3.84), (3.85) can thus be interpreted as corrections to 



the tree-level ground state estimates. The linear dependency is then not surprising, since 
selfconsistent propagation is affected by the quasiparticle mass of tree-level heavy modes in 



the tadpoles of figure 3.8 which in turn reflect the a priori estimate for the energy density 
of the thermal ground state. 

In calculating the contribution of higher irreducible loop orders to the polarization 
tensor we observe a hierarchic decrease and expect a termination at a finite order. For 
the effective theory of SU(2) Yang-Mills thermodynamics described in this chapter, it can 
be shown that the contribution of higher irreducible loop orders to the polarization tensor 
quickly decreases at higher loop orders |26[ 136] . Further contributions of these loop orders 
to the one-loop resummation calculated here should therefore be small enough to be safely 
ignored. 



3.6.3 Thermodynamical consistency 

We have seen that the a priori pressure and energy estimates for the ground state plus 
corrections, which incorporate fluctuations of the tree-level massive modes, obey the usual 
Legendre relations of the form p = T^ — P. What about the resummed corrections 
defined in fl3~80"l ), ( |3~8lj )? 

-ioop,rc S ag the Legenc i re transformation of AP 1 " 100 ^ 8 , 



Directly defining ApV 



Ap 



1 — loop. res 

L 



JAP 



1 — loop, res 



dT 



AP 



1 — loop, res 



(3.86) 



we can compare both the energy density A / 9^~ loop ' rcs resulting from the Legendre relation 

loop, res ca l cu lating the ratio of both func- 
Api/Ap. Such a comparison is plotted in figure 



(3.86) and the resummed energy density Ap 1 
tions A^ loop ' rcs / Ap 1 



-loop, res 

clear that the ratio is far from unity, and both densities show opposite signs. 



3.10 



It is 
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(a) The resummed one-loop dimensionless correction to the energy density, as denned in 
equation ( |3.80| l . Fitting of the data to a power law of the form ( |3.82[ | results in the expression 
Ap = 3.9577- A" 3 ' 02436 . 
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(b) The resummed one-loop dimensionless correction to the pressure, as defined in equation 
(3.811. Fitting of the data to a power law of the form (3.831 results in the expression 



Ap = 8.4963 ■ A~ 3 ' uuyu4 . 

Figure 3.9: The calculated high-temperature A dependence of (a) Ap and (b) AP. Dots 
correspond to computed values, the line is a numerical fit to the power functions defined 



in equations (3.82) and (3.83). 
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Figure 3.10: The A dependence of the ratio where Ap^ lo " i> ,!> 



«riioro A /i 1 lo °p. res = rpi jg defined in 
equation (j3.S6h , and A is the dimensionless temperature. The plot shows how the correc- 



tions to pressure and energy density due to one-loop resummation as calculated in the text 
are not thermodynamically self consistent, since radiative corrections shifting the dispersion 
law of tree-level massive modes were not considered in the calculation. 



How should we understand this apparent inconsistency? Remember that by using the 
bare (free) propagators of the TLM vector modes in our calculations, we ignored one-loop 
corrections to the propagation of these modes. Without taking into account radiative 
corrections for the massive modes, such inconsistencies were to be expected. One may stay 
consistent only up to a given loop order, and then every resummed propagator must be 
considered. 



Nevertheles, the linear correction given by expressions (3.82) and (3.83) with fitted 



parameters (3.84) and (3.85) can still be interpreted as the contribution to the thermal gas 



energy and pressure by photon propagation among stable monopole-antimonopole pairs, 
which in turn influence the ground state physics. 



Chapter 4 

Inhomogeneous Thermodynamics 



In the nonperturbative approach to the thermal SU(2) Yang-Mills gas discussed in chapter 
[3j the assumption of homogeneous and isotropic thermodynamics was of central impor- 
tance. In particular, it led to the inference that the modulus \(j)\ of the adjoint scalar field 
(p a should be constant (see section 3.1), and that the ground state field configuration 
must be pure gauge (see equation (3.10 )). The thermal ground state composed by the inert 
scalar field (j) and the pure gauge solution then serve as a homogeneous and isotropic 
background to the dynamics of topologically trivial fluctuations. 

Thus, the question arises as to whether this assumption can be relaxed. Namely, if it 
is possible to expand the nonperturbative treatment of chapter [3] to consider the case of 
inhomogeneous thermodynamics, and whether the physical description of a thermal ground 
state made up by a scalar field together with a background gauge field configuration remains 
valid. The second question that arises is, under what circumstances is such a description 
not valid? Is it always the case that a ground state (0, a^ s ) exists over which trivial 
fluctuations propagate? 

Both problems are not only of theoretical interest, but important for experimental con- 
siderations as well. If we hope to realize an SU(2) thermal photon gas in the laboratory, 
we need to specify under what experimental circumstances the theoretical formalism re- 
viewed in chapter [3] is valid. Further, and in the same vein, we need to know if there are 
any possible experimental conditions where the ground state description of the SU(2) gas 
breaks down. 

This chapter is an attempt to answer both questions and gain insight into the possibility 
of experimentally realizing the thermal SU(2) ground state. 

First we tackle the problem of generalizing the approach of chapter <^ to inhomoge- 
neous thermodynamics. A particular solution to the SU(2) equations of motion with scalar 
field has been found, which in the context of nonperturbative Yang- Mills thermodynamics 
and under the adiabatic approximation allows for temperature profiles without disturbing 
the ground state. The derivation of this solution is given in appendix |A| In this chapter, 



we use a practical ansatz for the spatial a vectors, then calculate with the help of (A. 17 



and (A. 24) the rest of the gauge and scalar fields. With a solution at hand, we investigate 
what kind of thermal gradients are produced, and under which boundary conditions. 

For the second problem, we invoke a recent resonant cavity experiment at low temper- 
atures by Tada et. al. |66j as an example of how the thermal SU(2) ground state may be 
deformed to be indistinguishable from the dynamics of a U(l) gas. 
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4.1 The SU(2) thermomagnetic effect 

In this section, we present a solution to the SU(2) Yang-Mills equations of motion that does 
not assume space homogeneity. Instead, the temperature profile is allowed to vary along a 
coordinate axis. As we will see, under certain assumptions, this temperature profile implies 
the existence of non-trivial gauge field configurations, giving rise to a magnetic field, which 
may be interpreted as emerging from the temperature gradient. This has some parallels 
with the thermoelectric effect in solids. Indeed, the SU(2) thermomagnetic effect can be 
understood as a rearrangement of topological magnetic charges in the ground state because 
of thermal diffusion. 



4.1.1 Basic approach 

We will abandon the requirement of constant temperature for the thermal system, and 
instead assume a space-dependent temperature distribution in one dimension, T(x) = T(z). 
Such a temperature distribution will in turn affect the thermal ground state through the 

modulus of the scalar field = |$|(^) = \J 2t!t(z) • Here we use the notation of 

appendix [A] and label the inhomogeneous scalar field as capital <£, to differentiate from 
the homogeneous inert field (j> of chapter |3| We are interested in the effects of a non-zero 
temperature gradient on the gauge field configurations which solve the Yang-Mills equation 
of motion. 

In the SU(2) gas, position-dependence of the modulus of <I> implies the position- 
dependence of the mass for the two massive modes. We expect that the pure-gauge ground 



state configurations of the form (3.10) will no longer be solutions of the equations of motion 
with space-dependent ground state, but will be replaced by a non-homogeneous configu- 
ration reacting to the deformation of <3?. The non-homogeneities in the gauge field lead to 
the emergence of non-scalar (isotropy breaking) quantities such as a magnetic field V x a 3 . 

The deformation of will also influence radiative corrections for the propagation of 
the transverse massless mode. Since radiative corrections are small in the effective theory 
compared to the tree-level (bare) ensemble, we may ignore such corrections. Therefore, 



we only search for classical solutions to the Yang-Mills equations of motion (3.9) with a 
space-dependent scalar field |$|(-z) taken as a source term. This is carried out in appendix 

El 



4.1.2 Adiabatic approximation 

An important consideration we must bear in mind is that the field energy density of such 



inhomogeneous field configurations described in section 4.1.1 must never be greater than the 

energy density of the SU(2) ground state over the entire domain of the fields. Such an energy 

hierarchy is condition sine qua non to ensure the existence of the coarse-grained scalar 

field. Otherwise, the ensemble dynamics could no longer be described by topologically 

trivial fluctuations over a static background, and the effective theory described in chapter 

[3] would not be applicable. We therefore demand that the energy scale of the scalar field 

$ be larger than any other scale of the system. This is the adiabatic approximation. 

In our particular solution to the equations of motion with z-dependent scalar field <£, 

l 
_■ 



a nontrivial magnetic field B with energy density p-Q = i | !B | 2 emerges. The adiabatic 
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approximation then demands 
p B (z)<p^(z). 



(4.1) 



Condition (4.1) guarantees that the effective theory summarized in chapter [3] remains a 
good description of the SU(2) ground state. 



4.1.3 Restrictions on possible solutions 

To use the results of appendix |A| we first need a good ansatz for the vector a 3 . To have a 
physical interpretation the solutions derived from the ansatz should be such that 

1. The temperature field T(z) has a parametrization such that a baseline temperature 
T m i n = lim^oo T{z) can always be specified, with T m ; n greater than the critical 
temperature T c of the deconfining phase, T m ; n > T c ~ 2.73 K. 

2. The conditions needed for the validity of the adiabatic approximation, described in 
section 



4.1.2 particularly equation (4.1), must hold. 



Note that condition [l] rules out using simple power laws such as a 3 ~ z 7 as ansatz, 
since after inserting in equations (A. 17') and (A.24), the temperature field always shows 
~ z -2 " 7 which is either or infinity, depending on the sign of 7. One 



the limit lirn^oo T(z) 
could try numerically searching for solutions a 3 (z) of the equation 



1 



3\2 



(a 3 ) 




A 3 
2vr 



T-\z) 



(4.2) 



Equation (4.2) is arrived at by inserting (A. 17') in (A. 25) 



for a specific distribution T(z) 
and taking a 3 = a 3 , = a 3 . 

But this was too much of an overhead for our purposes. In the next section, we will 
offer a simple ansatz which fulfills both requirements [T] and [2] using just two parameters. 
We then discuss how it describes deformations from homogeneous thermodynamics. 



4.1.4 A particular ansatz 

In this section, we will introduce a simple rational function as an ansatz to derive in- 
homogeneous field and temperature configurations. We will check that the requirements 



outlined in section 4.1.3 be fulfilled, then discuss the derived field solutions, in particular 
the resulting temperature distribution. This exercise will help us build a basic intuition for 
the possible ways that an homogeneous temperature distribution may be deformed without 
disturbing the ground state physics of the SU(2) gas. 

In the following, unless otherwise stated, we work in natural units h = ks = c = 1, 
so energy and temperature have dimension eV, length dimension (eV) -1 , etc. We use 
dimensionless coordinate variables, which are made dimensionless by rescaling with the 
factor k^- (the factor 2ir is for numerical convenience only). Recall that A is the SU(2) Yang- 
Mills scale and has dimension of energy, [A] = eV, so lengths are rescaled by multiplication 
with j-. Coordinate variables are then rescaled as z 



2tt' 



Define now the z-dependent function 



9 (z; 90, 9i) 



9o 



2tt) z 2 + 1 



(4.3) 
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with dimensionless parameters go and g\. By definition, the function 4.3 has dimension of 
energy squared, [g] = (eV) 2 . 

From appendix |a| we learn that from a reasonable ansatz for a 3 (x) = a 3 (z) we can 
calculate a complete solution to the equations of motion (A.l). Using (4.3) as an ansatz, 
we write 



a\(z) 



4{z) 



a 3 (z) 



9(z) 



(4.4) 



Note that the energy dimensions on the right and left hand side of (4.4 ) coincide. Definition 



(4.4) requires the parameter g$ to be positive. We also take g\ > so as to tentatively 



allow all z > 0. From (A. 17'), (A. 24) and (3.2), we can calculate 



a (z) 



\t>\ 2 (z) 



T(z) 



dh 3 



A 



2tt V 2e 2 a 3 ' 



1 



■i\2 



A\ 2 dja 



2vr J 4e 2 a ' 



A 3 



27r|$l 2 ' 



(4.5) 
(4.6) 
(4.7) 



The magnetic field for this field configuration can be calculated from (4.4) 



\B\(z) = Iroto 3 ! = \/2— dza 3 , 
I I 2vr 



PB (z) 



1 



IBI 2 . 



(4.8) 
(4.9) 



Finally, recall the energy density of the thermal ground state, which is given by equation 



(3.11) 



pg* (z) = 4vrA 3 T (z) . 



(3.11) 



A z-l 



Notice that lim^^oo |$| 2 = lim^^oo ^(a 3 ) 2 = (^) 2 go, since from (4.5) we read linig 

dja 



z and lim; 



0. Then from equation (4.7) it is seen that 



lim T{z) 

z— >oo 



(4.10) 



So the baseline temperature at infinity can be fixed by adjusting g$. This fulfills our first 
condition. In order to have T m i n > T c we need 



%Tikg„ l > T c 



A 

2?r' 



go < 



16vr 2 
^7 



(4.11) 



For A c = 13.87 [31J, this means < 11-39. We also demand that the temperature profile 
decrease asymptotically to T c , but this means that a 3 must be monotonically increasing 



for large z. From (4.4) and the definition of g (4.3) we see that this can be achieved by 
having g\ < 1. 
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The second condition, that the magnetic field energy density always be smaller than 
that of the ground state, is difficult to check analytically for all z. But we can at least 
make sure that it holds asymptotically. Recall that in the limit of large z, |$| 2 behaves as 
^(a 3 ) 2 . Then the ground state energy density at large z is given as 



lim p s 

z— >oo 



32tt 2 AV 



The magnetic field energy density at large z can also be calculated 



lim p B (I) 

2— »OC 



A V 



— I hm(9 5 a 3 ) 2 

Z7T / z^co 



9o 



2vr J 



Condition (|4.1l) then reads 



go < (\/27r)167r 2 



(4.12) 



(4.13) 



(4.14) 



But y/2iv > X c 1 , so the restriction T m j n > T c already ensures the energy hierarchy required 



by the adiabatic condition in the large z limit. For small z, we will need to test (4.14) 
explicitly. 



Summarizing, we see that for < g\ < 167r 2 A c 1 and < g\ < 1, the ansatz (4.14) 



inserted in (4.4 ) results in field configurations that do indeed fulfill the conditions in section 



4.1.3 at least asymptotically. We test these configurations in the next section. 



4.1.5 Some examples 

Now we try our solu tion using ansatz parameters go 
are plotted in figure 



4.1 



gi = 0.1. The resulting fields 
in units of the SU(2) Yang- Mills scale ^~ an d m dimensionless 

coordinates. 

Note how a 3 (z) (figure 4.1a) tends to a constant value at infinity. Plotting ao{z) 

(figure 4.1b), we encounter our first problem, clq(z) becomes imaginary at Zq = 0.39, 
which is also the point where the real part of ao(z) becomes zero. The temperature field 
(figure 4.1c) shows non-monotonic behaviour in the region where Reao(i) = 0. It assumes 
0, which contradicts the requirement that T{z) 



This is not 



a minima at T 

surprising. In section 4.1.4 we showed that the function g, defined in (4.3) and used as 



ansatz for the field a 3 fulfills the requirements outlined in section 4.1.3 asymptotically. We 
need to restrict the domain of our solutions from z £ [0, oo] to z £ [Zrmm °°] with z m i n large 
enough that the problems outlined above disappear. 



Figure 4.1c shows two maxima for the temperature field T(z). After the second maxima, 
T{z) appears to decrease asymptotically to a constant value, which is the behaviour we 
are looking for. It seems therefore appropriate to define z m i n to be the value of z at which 
the second maxima of T(z) occurs, since this would give us a monotonically decreasing 
temperature distribution. In this case then, z ma x = 0.836. Note also that z max > zq. The 
domain [z m i n ,oo] should therefore be free of imaginary fields. 

For simplicity, we define a coordinate translation z — > z — z max , and consider only 
z > 0. This is the region we are interested in from now on. 



In figure 4.2 we have plotted the translated fields ao and T. We interpret figure 4.2b 



physically as a temperature gradient. In figure 4.2c we have plotted this gradient in SI units. 
We see that for parameters go = 1 and g\ = 0.1, the resulting temperature distribution 
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"a 0.4 



< 0.2- 



(a) The field a 3 (z) plotted in units of the SU(2) 
Yang-Mills scale =p, calculated from 14.41) 



(b) The field a (z) plotted in units of the SU(2) 
Yang-Mills scale ~, calculated from ( |4.5[ ). The 
real part of ao is plotted in black, the imaginary 
part in blue. 



1200 



< 




(c) The temperature field, plotted in units of the SU(2) Yang-Mills scale calculated 



from (4.7 1. The maxima of T(z) in the domain where Im ao(z) = is shown by the 



blue line and lies at z max — 0.836. 



Figure 4.1: Various fields calculated from ansatz (4.3) with parameters go = 1, g% = 0.1. 
z is dimensionless. 
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(a) The field ao(z) plotted in units of the SU(2) 
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(b) The temperature field T(z) plotted in units of 
the SU(2) Yang-Mills scale after translation 
z — > z — Zmax so that it has a maxima 5 = 0. 
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(c) The temperature field T(z) from figure 
maximum value of T max = 59.26K and asymptotically decreases to a value T m 



plotted in SI units. It starts from a 



31.04K. This gradient is a result of the ansatz in (4.3 1 with parameters go = 1 and 
<7i=0.1. 

Figure 4.2: The fields (a) ao(z) and (b) T{z) are translated in z so that the maximum of 
the temperature now lies at z = 0. The ansatz parameters are the same as in figure 4.1 
go = 1, gi = 0.1. z is dimensionless. In (c) we have plotted the temperature field in SI 
units. It decreases from a, maxima Tmax — 59.26K to a baseline value Xmin — 

31.04K. 
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decreases from a maxima T max — 59.26K to a baseline value T m [ n — 31.04K. The length 
scale of this gradient is of the order of centimeters. 

This temperature gradient is a scalar perturbation of the equations of motion for the 
SU(2) fields, and causes the field configurations (4.4), (4.5), which are static nonzero solu- 
tions of the perturbed equations of motion. Then we can interpret the resulting magnetic 
field B(z), plotted in figure 4.3a as a thermomagnetic field generated by perturbing the 
coarse-grained SU(2) Yang-Mills thermal ground state with a temperature gradient of the 
form (4.7). The interpretation that the thermal gradient causes the magnetic field and not 
the other way around is forced upon us by the requirement that the ground state energy 
scale, and therefore the thermal gradient, is the dominant energy scale of the system. Note 
from figure 4.3a that it goes to zero very quickly. 

The thermomagnetic field in figure 4.3a contributes to the energy density of the (per- 
turbed) vacuum with ps as given in equation (4.9 ). This energy density is plotted in figure 
in units of the SU(2) Yang-Mills scale Comparing ps with the energy density of 



4.3b 



the ground state (figure 4.3c), we see ps lies about 8 orders of magnitude below p 9S over the 
entire domain. The condition (4.1) imposed by the adiabatic approximation is satisfied. 
Solutions calculated from go ~ 1 and g\ ~ 0.1 are therefore acceptable deformations of the 
thermodynamically homogeneous SU(2) Yang-Mills ground state. 

To show the effect of parameter g\, we increase g\ while leaving go constant, and 
calculate another set of solutions using go = 1 and g\ = 0.999. We have plotted the 
resulting energy density and temperature profile in figure 4.4 Note that the baseline 
temperature T m i n = 31.04 K remains unchanged, since it only depends on go as was shown 
before. The new maximum, however, lies at T max = 31.67K. Increasing g\ has resulted 
in a smaller value T max . This parameter is therefore responsible for setting the maximal 
temperature of the profile. It can also be checked numerically that in the limit g\ — )■ 1, 
d z T(z) — > and the temperature is constant over the entire domain, with T(z) = T max = 
Tmin = 87rAgQ . Ansatz (4.3) shows a continuous limit (numerically) to the homogeneous 
case. 

Finally, we remark that in the context of nonperturbative SU(2) thermodynamics, a 
dual interpretation is necessary [25J, in which the magnetic fields in the effective theory 
are understood as electric fields in nature. The thermomagnetic effect described here is 
therefore a thermoelectric effect in nature. 



4.2 Decoupling from the ground state 

We have seen how it is possible to disturb the homogeneity of the thermal ground state 
in such a way that it still remains an effective description of the thermal SU(2) gas. The 
deciding criterion is the energy scale hierarchy. As long as the ground state energy density 
is the dominant scale, deformations from the totally homogeneous case are still consistent 
with the existence of a coarse-grained ground state. In our particular solution, such distor- 
tions were also accompanied by non-trivial gauge field configurations responsible for the 
emergence of a magnetic field. In order to remain consistent with the adiabatic approx- 
imation, the field energy of such configurations must always be considerable lower than 
that of the thermal ground state. 

In our experimental endeavours, for reasons discussed in chapter [5] it will not only 
prove useful to know under what conditions the temperature or field energy deformations 
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(a) The modulus of the magnetic field B(z) plot- 
ted in units of the SU(2) Yang-Mills scale cal- 



culated from (4.8 1. The coordinate z is plotted in 
meters. 
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(b) The magnetic energy density ps(i) plotted in 
units of the SU(2) Yang-Mills scale calculated 



from (4.91. The coordinate z is plotted in meters. 
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(c) The magnetic field energy density Pb(z) (in red) calculated from the generated 
magnetic field B(z), and the energy density of the thermal ground state p gs (T(z)) (in 
black), plotted in units of the SU(2) Yang-Mills scale p 9 " lies roughly 8 orders of 
magnitude above ps- 

Figure 4.3: The magnetic field B(z), its energy density ps(-S), and a comparison with the 



ground state energy density The ansatz parameters are the same as in figure 4.1 

go = 1, gi = 0.1. The coordinate z is plotted in meters. 
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(a) The temperature field T(z) plotted in SI units. It starts from a maximum value of 
Tmax = 31.67K and asymptotically decreases to a value T m i n = 31.04K. This gradient 
is a result of the ansatz in (4.31 with parameters go = 1 and g\ — 0.999. 




0.20 



(b) The magnetic field energy density Pb(z) (in red) calculated from the generated 
magnetic field B(z), and the energy density of the thermal ground state p gs (T(z)) 
plotted in natural units. 



Figure 4.4: (a) Temperature field and (b) energy densities calculated from ansatz (4.3) 
with parameters go = 1, g\ = 0.999. The coordinate z is plotted in meters. 
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conserve the ground state, but also under what conditions they cause a decoupling from 
it, effectively creating U(l) conditions. In this section, we investigate the possibilities of 
creating an effective U(l) system. 



4.2.1 Electromagnetic field decoupling 

As explained in the preceeding section, the energy scale of the ground state energy must be 
the largest scale of the system if the effective theory is to remain valid. But what happens 
when this hierarchy is destroyed, and the SU(2) ground state energy density becomes small 
relative to another energy scale? 

An example of such an experimental situation is given by the single-photon counting 
experiment described in [66| . 85 Rb atoms in the Rydberg state flying trough a 

resonant cavity are excited upwards to the lll P3/2 state by absorption of thermal photons 
of frequency 2527 MHz, and then are counted by a method of selective ionization. Thus 
the experiment is sensible to the mean photon number n(T) inside the cavity at a given 
temperature and cavity resonant frequency. The measurement was performed for temper- 
atures T ranging from 67 mK up to IK. During the measurements, a static electric stray 
field of magnitude \E\ ~ 25 was present in the cavity. No deviation of n(T) from the 
U(l) expected Bose-Einstein distribution was observed. Yet under the SU(2) hypothesis, 
absolutely no photon propagation should be possible. 

From what we have learnt in this chapter, we can guess that this could have to do with 
the energy hierarchy inside the cavity. The energy density pe of the external electric field 
in SI units is given by 

PE = ^E\ (4.15) 

where e = 8.854187817 10 -12 J/(V 2 m). The energy density of the SU(2) thermal ground 
state p EB in SI units is 

pT = 4* A* -^T , (4.16) 



where A = 2tt k ^ Tc , T c = 2.725 K and A c = 13.87, see section 3.4 To maintain energy 
scale hierarchy, the electric field energy density should be smaller than the energy density 
of the ground state, pe < p ss - For an static electric field magnitude of \E\ ~ 25mV/cm, 
this would imply a ground state temperature of the order of at least 10 3 K. That is, 
the effective radiation temperature implied by the stray electric field \E\ is at least three 
orders of magnitude above the thermalized radiation at 1 K (and below) inside the cavity. 
With such a large energy disparity, a decoupling of the ground state physics from the 
propagation properties of photons must have occured. Ground state induced effects on 
the photon disappear and the dispersion law becomes indistinguishable from that of the 
conventional U(l) theory. 

As we will see in chapter [5j this decoupling effect can be used to compare and test for 
the presence of an SU(2) thermal ground state. 



4.3 Summary 

In this section we have investigated the different possibilities of deforming the thermal 
ground state. We have learnt that energy scale hierarchy determines the validity of the 
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effective theory and thus the presence of SU(2) effects. 

To realize SU(2) black cavities experimentally, we should therefore be wary of stray 
electric or magnetic fields with energy densities much larger than that of the effective 
ground state such as those present in the cavity experiments of Tada et. al. |66j . For 
example, at T = 10 K, the energy density of the ground state becomes comparable to that 
of an electric field of magnitude \E\ ~ 2.47mV/cm. To protect scale dominance of the 
thermal ground state, the cavity should be shielded from stray electromagnetic fields much 
smaller than 2.47mV/cm, of the order magnitude of at most \E\ ~ 10~ 3 mV/cm. 

Similarly, to mirror the conditions of the Tada et. al. experiment and produce a 
black U(l) emitter at T = 10 K, the cavity should be transcended with electric fields of 
magnitude of at least \E\ ~ 10 3 mV/cm. 

This knowledge will be useful later when we design an experiment to compare between 
both regimes. 



Chapter 5 

Radiometry and Bolometry of SU(2) 
Photons 



In this chapter we first review some basic definitions regarding energy transfer, with the 
goal of arriving at an expression for the spectral black-body radiance without assuming a 
U(l) dispersion law. Next we examine the basic concepts of antenna theory and re-derive a 
formula for the antenna temperature. Again, this is done without making any dispersion law 
explicit. Both sections are included in this work since, apart from reviewing basic concepts, 
we need to check that the usual formulae given in the literature also apply in the case of a 
non-trivial dispersion law. This will, indeed, turn out to be the case (up to the appearance 
of a screening region, as described in chapter |3j), but it is not a priori clear and needs to 
be shown explicitly. The third section describes briefly what a radiometer system is, and, 
most importantly, the possible sources of measurement error and the maximum possible 



radiometer sensitivity. Finally, sections 5.4 and 5.5 discuss the possibility of detecting the 
SU(2) thermal ground state by radiometric means and possible experimental realizations, 
relying on concepts introduced throughout the chapter. 

Unless stated otherwise, we always use natural units, h = c = kg = 1, until sec- 



tion 5.4 in which we describe calculated predictions of experimental results using SI units 



exclusively. 



5.1 Thermal radiation and energy propagation 

Following |22[ 170] . in this section we review basic concepts on radiation propagation. The 
results derived in this chapter will be important when characterising thermal energy trans- 
fer. Specifically the expression for black-body radiance is obtained without making use of 
any explicit dispersion law. The SI units of important radiometric quantities introduced 
in this section will be written beside their equation number. 

Recall that the fundamental self-interacting SU(2) gauge modes become weakly inter- 
acting thermal quasiparticles in the effective theory \27\ |4"T] . This allows the local spectral 
decomposition of any quantity X characterising some aspect of energy transfer by photon 
propagation. The spectral components or X^ of X can always be defined by only 
considering contributions from monochromatic radiation of mode k or frequency w. Then 
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Figure 5.1: The flux through an area element da propagating along the subtended angle 
dQ, see text. 

the total quantity X on a system is given by the integration over its available modes 

X = JdkX k = JdojX w . (5.1) 

By a change of variable we can relate both fc-mode and w-mode spectral decompositions 
X u = \d k u k \- 1 X k . (5.2) 

5.1.1 Radiant power 

In the following we will deal with the characterisation and measurement of the electromag- 
netic field in vacuum. The labels energy or radiant energy will always refer to the energy 
stored in the propagating photon field as reviewed in Appendix [B] Radiant power $ is then 
defined as the change in radiant energy Q per unit time 

* = • [W] (5 - 3) 

We will also often refer to "energy rays" or "energy beams". The ray concept, borrowed 
from geometrical optics, provides an intuitive characterisation for the propagation of energy 
from source to detector. This intuitive description is made precise by identifying "ray" with 
the Poynting vector of the propagating electromagnetic field. When expanding in single 
Fourier modes, the Poynting Vector is parallel to the wave vector k of the plane wave, see 
section IB .21 

5.1.2 Radiance 



Consider now an elementary beam of incoherent radiant power $ in a homogeneous and 
isotropic medium passing through a surface element da, which propagates further within 
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the solid angle d£l with axis vector n, see Fig. 5.1 Then the radiance (or brightness) at a 
point in da, defined as the radiant power per unit solid angle (in the direction of the ray) 
and per unit projected area (perpendicular to the ray), can be written as 

L " dnW^) = dn( d dalo S 9)- [Wm^sr- 1 ] (5.4) 

Consider now two area elements da\,dai, separated by a distance R, through which an 
incoherent beam of radiant energy propagates along a ray vector n. Let 9±, 02 be the angles 
between n and each of the surface's normals. Then the solid angle subtended by each area 
element at the opposite end is 

ddi ■ h dai cos 0j . 

= -hs- = — hs— > ' = 1 ' 2 - 

It follows that dVt\da2 cos 9% = d^dai cos#i. 

Since it was assumed that the beam passing through do\ and da2 is contained entirely 
within each area element, we conclude from energy conservation that the power d& passing 
through both areas is the same. Then d 2 Q± = d 2< 3?2 and one obtains for the radiance Lj at 
each element 

_ d 2 <5> _ d 2 $ _ 
1 dQ,2 {da\ cos 9\) d£li (da2 cos 02) 2 

L is thus conserved along a ray within an homogeneous and isotropic medium. This is why 
radiance is a useful concept for describing energy transfer. 

5.1.3 Transfer of radiant energy 

Let a radiating source be characterised by its radiance L. Then the power flux S emitted by 
the source of radiance L and passing through an area element da is given by an integration 
over the solid angle Q s subtended by the source 



S = dQ cos9L(9,(f)) , (5.5) 
Jn s 

where 9 is the angle between the normal da to the area element and the axis h of f2 s , 



compare with Fig. 5.1 Accordingly, the total power is found by integrating S a over a total 



area A of a detecting device 

$= f daS a . (5.6) 



J A 

Note that the linear superposition of power is only possible by assuming no interference 
effects, that is, by assuming incoherent radiation. Since we will be dealing with thermalized 
radiation, this assumption will be met ideally. 

5.1.4 Radiance formula for thermalized black-body emitters 

We will now derive the known Planck radiance formula for black-body emitters in thermal 
equilibrium with radiation |54j . Special attention will be paid to the assumptions under 
which the derivation holds, and no special dispersion law will be singled out. We believe it 
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Figure 5.2: Spherical volumes Vi and V 2 

necessary to do this calculation in detail since most presentations in the literature assume 
U(l) dispersion relations. The dispersion law for the U(l) photon is just the linear function 
= c |k|, but in the SU(2) case, Wk = \/k 2 + G and it is not clear that the classical U(l) 
result for the black-body radiance remains valid, since the functional dependence of the 
energy density for each Fourier mode of its wave vector is nontrivial (see section 3.5). In 



the following, in order to write thermalized radiation as a sum over plane waves, we will 
also need the results from Appendix [B] 

Consider a spherical black-body surface volume element dVi of radius T2 positioned 
at the centre of an isothermal spherical black enclosure V\ at the same temperature, and 
let the surfaces dV\ , d V2 be thermally coupled by exchange of radiant energy. Thermal 
coupling must hold for every available mode of the cavity, so thermal equilibrium must hold 
for available mode individually. Assuming energy and charge conservation, the continuity 
equation 



(B.58) 



dt d 6 xp+ d> d a • S = , 

Jv JdV 

must hold, and the radiant power emitted by 8V2 in the mode corresponding to the wave 
vector k and collected by 8V\ must be 



k 2->l 



-d t Q 



k 2 



-dt 



d 3 x/9 k 



d a ■ Sk . 



(5.7) 



dV 2 



Here the integration volume V 2 C is the volume bounded by dV<i and infinity. From the 
symmetry of the problem and thermodynamical equilibrium, the Poynting vector S of the 
thermalized electromagnetic field at dV2 is normal to the surface, and constant everywhere. 



This allows us to write the spectral expression (5.7) as 



$k 



2->-l 



|Sk| 



dV 2 



da = I Sk| 47T72 , 



(5.8) 



where |Sk| is the magnitude of the Poynting vector for the mode k at the boundary 8V2- 
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Consider now the spectral power ^k i->2 emitted by dV\ and absorbed at dVi- Using 



(5.6), we may write it as 



$ki^2= da tin cos n L k (#,</>) , (5.9) 
Jav2 Jn ± 

where 6 n is the angle between the normals of each area element at dV2 and the solid angle 
Q± = 411 subtended by dV\. From spherical symmetry, 9 n = 0. Further, it is easy to 
show |22j that for the black-body emitter, its spectral radiance Z/kftft must, by definition, 
be independent of the polar angles 9, (f) 

Ac {0, <t>) = L k bb (e,<f>) = L k bb (0,0) = L k bb ■ (5.10) 
We may then write 

$ki-+2 = / da dQ L k bb = 47r L k bb / da = lQit 2 L k bb r\ . (5.11) 
JdV 2 J fix JdV 2 

Since we have assumed thermalization for each available mode k in the cavity, it follows 
that $k 2-s-i = ^ki-s-21 arid we obtain a general expression for the spectral black-body 
radiance Lk bb m a thermalized cavity 

Lkbb = -;— |Sk| • (5.12) 

47T 

Our conclusion is also valid for any other geometrical configuration. This can be argued 
as follows. 

In general, cavity geometry influences the homogeneity and isotropy of thermal radia- 
tion inside the enclosure. For a rectangular cavity of side lengths Lj, i = 1,2, 3, and volume 
V = L1L2L3, for example, the correction 5D{v) caused by geometric discretization effects 
to the ideal spectral density Dq(v) = (^) of electromagnetic modes in a lossless closed 
cavity of same volume V is given by the expression (in SI units) |U 

5D(v) /c\2 1 

- ' ' ;(Li + L 2 + L 3 ). (5.13) 



Dq(v) \vJ %W 

But by choosing the geometric dimensions of the cavity to be large enough, these dis- 
cretization effects can be made arbitrary small in the spectral region of interest. With this 
in mind, here and in the following we will always presume the cavity to be large enough 
so that homogeneity and isotropy may be assumed up to any order. 

Take now a black-body cavity with arbitrary geometry, and choose an spherical in- 
finitesimal volume element dVo positioned at a point ro inside the enclosure. In thermal 
equilibrium, the radiation field is homogeneous and isotropic, so that the propagation 
of radiation from the cavity walls through dVo may be modelled by assuming an outer 
spherical surface So centered in ro emitting radially at oIVq. By the same argument, since 
dVo encloses homogeneous isotropically propagating thermal radiation, it may be seen as 
a spherical black-body emitting at So. But this results in exactly the same geometrical 
configuration as before. We conclude that for a cavity with arbitrary geometry and suffi- 
ciently large dimensions^ that discretization effects can be ignored, with walls in thermal 



equilibrium with the radiation field equation (5.12) holds. 



1 What '"sufficiently large'" means depends of course on the cavity geometry. In this work, we restrict 
ourselves to rectangular cavities, where corrections are given by | pU3] l, 
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We now use the results of Appendix [B] to express L^b as an integral over Fourier modes. 
Recall the expression for the energy density p of thermalized U(l) photon radiation at 
temperature T = first derived by Planck |54j (here written in natural units h = c = 
k B = l) 



P(P) 



d 3 k 



(2vr) 3 e^* - 1 



(5.14) 



To describe the both the thermalized U(l) and SU(2) photon field using the expressions 
derived in Appendix [Bj we demand that equation 

P=% = l\- r \ ^3 (2A^ k ) u k , ( |R47} ) 



V 



(2tt) £ 



which gives the energy density of a given U(l) photon configuration in vacuum (defined by 
Nfe), to coincide with (5.14). The factor 2 in (B.47 ) is due to summing over two transverse 



polarization states. The energy density of the SU(2) photon gas is given by (B.59), which 
sums only over transverse polarization states, since the longitudinal polarization, which has 



a very narrow wave vector support, cannot transport energy (see section 3.5.5). Equation 



(B.59) has the same form of (B.47 J, but with w k given by the SU(2) transverse dispersion 
law (3.34) (written here in SI units) 



(/ia; p ) 2 = (cp) 2 + G>,p;r,A) . 



[J 2 ] (5-15) 



That is, we can use the same spectral decomposition (B.47 ) when describing the total field 



energy in both theories, with the only replacement being the non-trivial dispersion law for 
the photon of SU(2) Yang-Mills thermodynamics (see chapter pj. This is important, so 
we will repeat it here: We use the same expression (B.47 ) for both the U(l) and the non- 



abelian SU(2) Yang- Mills photon, the vast dynamical difference between both is already 
taken into account by the different dispersion laws. 



Comparing both (5.14) and (B.47 ) gives us the correct thermal normalization factors 
for the spectral decomposition (B.47 J 



N k 



IV n B (/So*) 



2w k 



(5.16) 



where ns(/3o;k) = (e^ k — l) 1 is the Bose-Einstein distribution function. 
The spectral energy density is then given as 

q r r d 3 k 

V " 



dkpb 



1 



dk k 2 ns(/3wk) 



(2-iry TT- 
=>• Pk = -4 n B (/3wk)k 2 w k - 
The energy propagation velocity Vk of a mode k can be defined [7] so that 
Pk v k = S k • 



(5.17) 



(5.18) 



For each mode, it can be shown that this propagation velocity, by definition, is just the 
fc-gradient of the dispersion law [7] 



v k = V k Wk , 



(5.19) 
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whatever the functional expression for the dispersion law is. 



Inserting (5.17), (5.18), and (5.19) into (5.12), we obtain the spectral radiance of a 



black-body emitter in thermalized radiation 



L\t bb 



1 

4tt 



I Pk V k o; k | 



1 



471-3 e /3^ k _ I 



k 2 |<9kW k | 



(5.20) 



Using (5.2 ), we can also write the same expression as a function of frequency. The radiance 



per frequency interval of a black-body emitter in thermalized radiation is then, written in 
natural units, given as 



bb 



w k 



1 



4^3 e /3w k _ 1 



[eV 4 ] (5.21) 



To write (5.21) in SI units, remember that the (angular) frequency and wave vector are 
rescaled as w k — > 7iw k , k — > cftk and the inverse temperature as /3 — > jr^f ■ Moreover, 



radiance has dimension eV 4 in natural units, so we need to divide expression (5.21) with 
(he) 2 to return to SI units. We obtain 



hv 



bb 



271-2 e hu/k B T _ l 



[Wm~ 



sr 



^Hz" 1 ] (5.22) 



The additional factor of 2tt comes from the change of variable doj —> 2irdv. The formulae 



(5.21), (5.22) characterise radiation in thermal equilibrium with a black-body, indepen- 



dently of any functional relation between k and w k . 



Later, (5.22) will be used to calculate the power output of any measuring device ob- 



serving thermal radiation. 



5.1.5 Radiance of the SU(2) thermal gas 



In this section we compare radiance functions of both the U(l) and SU(2) theories. The 
radiance function L^)(u,T) for the U(l) photon gas in SI units can be read off from Eq. 
(5.22). Inserting the U(l) dispersion law hw^ = ch\\t\ = hv results in the expression 



2h 

c 2 e 



hu/k B T _ i 



[Wm^sr'Hz" 1 ] (5.23) 



To write the radiance function 



v, T) for the SU(2) case, we recall the SU(2) disper- 
y/(chk) 2 + G(oj, Kk; T, A). Recall that the SU(2) photon gas 
shows a screening region so that propagation is forbidden below a cutoff frequency v* , see 



sion relation (3.34) faj k 



section 3.5 That is, the only modes allowed in the cavity are those with frequencies v > v* 



That only those modes can contribute to the spectral radiance is ensured by multiplying 
with a ^-function. In the end, we obtain 



2h v 6 

e hv/k B T _ i 

c 2 G 



L u bb 



1 



[huf 



c 2 G 
(hv) 2 

0(v-v 



[Wm^sr^Hz -1 ! 



(5.24) 
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Note that outside the screening region, that is, for v > v* , the radiance (5.24) for the SU(2) 

times the U(l) radiance (5.23). As discussed in [461 



c 2 G 



gas factorizes into the factor (1 - 

[59], the function G goes to zero power-like for increasing temperature and exponentially for 
increasing frequency. So we expect a return to the U(l) case for high enough temperatures 



and energies, since in that limit the factor (1 



1, V* 



0, and both radiance 



functions L\l(y,T) and L u A (v,T) become indistinguishable. 

In the case of low energy and temperature, however, we do expect the SU(2) photon 
gas to show considerable deviation from the U(l) case. This suggests the possibility of 
experimental detection of the SU(2) gas by radiation measurements at low frequency and 



temperature. This possibility will be explored further in section 5.4 



Examples of the calculated radiance for different temperatures are plotted in figures 
5.3 to 5.6 To plot the SU(2) radiance, the cutoff frequency v* was calculated as shown in 



section 3.5.4 Note the screening region at low frequencies where the radiance of the SU(2) 
gas vanishes. The crossover frequency u c , discussed in section 3.5.4 is also easily visible in 
all plots. 



5.2 Basics of the theory of antennas 

In order to test experimentally the most important prediction of the SU(2) photon theory, 
and thus ascertain the presence of a thermodynamical ground state composed by topolog- 
ical field configurations as described in chapter |3j we need to be able to detect thermal 
radiation at low frequencies and temperatures. In the preceding section we derived an 
expression for the radiance of black-body emitters which characterises energy propagation. 
Now we turn to the problem of how to measure such an altered propagation. 

The experimental science and technology of the measurement of electromagnetic radi- 
ation is called radiometry |67|. A typical radiometric set up always consists of a radiation 
source, an aperture or optical system, and a detector. The aperture collects the radiation 
from the source and redirects it to the detector. The detector converts radiation into a well 
defined signal that can then be measured within a certain precision and further processed. 

In the microwave region of the electromagnetic spectrum - which, as we will later see, 
is where the predicted differences between U(l) and SU(2) propagation occur - an antenna 
is often used as the aperture. An antenna is defined |40|, I6T] as the region of transition 
between an electromagnetic wave propagating in free space and a guided wave propagating 
in a transmission line 

In this section, we review some concepts relating to energy collection through antennas. 
The most important result will be an expression for the antenna temperature Ta, which 
quantifies the power output of an antenna in relation to its incident thermal radiation, 
derived here without assuming any explicit dispersion law. The detection of the collected 
thermal radiation is the subject of section [5~3] 

Further properties that characterise antennas are geometric in nature, such as the 
antenna radiation pattern and the antenna effective area. Geometric factors also determine 
such coefficients as beam efficiency, which characterises undesired signal contribution by 
side-lobe radiation. These characteristics will play a minor role for us since we will be 

2 A transmission line is a device allowing for propagation of radio-frequency energy from one point to 
another. Antennas radiate or receive energy, while transmission lines transport energy. See also section 
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Figure 5.3: Comparison of the spectral radiance of an SU(2) photon gas (red) and a U(l) 
gas (grey) at a temperature T=5.4K. Note the spectral gap in the low-frequency region. 
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Figure 5.4: Zoom-in of figure 
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T=8.0K 




10 20 30 40 50 



v[GHz] 

Figure 5.5: Comparison of the spectral radiance of an SU(2) photon gas (red) and a U(l) 
gas (grey) at a temperature T=8.0K. Note the spectral gap in the low frequency region. 
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Figure 5.6: Comparison of the spectral radiance of an SU(2) photon gas (red) and a U(l) 
gas (grey) at a temperature T=11.0K. Note the spectral gap in the low frequency region. 
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dealing with isotropic and homogeneous radiation. Still, these geometric factors need to 
be understood before we arrive at a formula for T4. We briefly review these subjects before 
arriving at our main expression for the antenna temperature. 



5.2.1 Antenna parameters 

This section serves as a quick introduction to all antenna parameters we will require in the 
following sections to describe and calculate through our experiment. A thorough review of 
the subject can be found in [40j. 

Consider a transmitting antenna placed at the origin. For sufficiently large r, we may 
consider the radiated modes as plane waves. This zone is called the far field, and we always 
work in this limit unless stated otherwise. 

With the antenna positioned at the origin, we expect the far field electric and magnetic 
vectors to be mutually perpendicular and normal to the unit vector f. The Poynting vector 
S must be radial. In the case of isotropic antenna radiation, the total radiated power given 
by integrating over the surface of a spherical volume V centered at the origin becomes in 
this limit 

$ = / da S = 4nr 2 S r , (5.25) 

JdV 

since the radial component S r is constant over all directions. In general, antenna radiation 
patterns are not isotropic, and almost always distributed around a narrow solid angle called 
the main lobe. S r is no longer isotropic and the polar dependency of antenna radiation is 
characterised by the antenna radiation pattern 

F(r,8,cj ) ) = r 2 S r (9,ci ) ) , (5.26) 



so that the total radiated power in the far field is given by 



$ = J anF(e, 



(5.27) 



where we have omitted the r-dependence, which in the far-field reduces to equation (5.26) 
F{9,4>) = T 2 S r . For an isotropic antenna, (5.27) reduces, of course, to (5.25). 



Imagine now a radiating antenna which we supply with input power Pi. From the 
second law of thermodynamics we may assume that a fraction of P, is dissipated as heat 
from the antenna. Let P Q be the actual radiated power from the antenna aperture, and 
define the antenna radiation efficiency r\ as the ratio P Q /Pi. We can define a polar version 
of this factor as well, the antenna gain function G(9,(p), as the ratio between radiated 
power per solid angle d&/d£l = F (9, (j>) and total input power per solid angle Pi /An 



4vrF( 



Pi 



(5.28) 



Defining the normalized antenna radiation pattern 



(5.29) 
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the pattern solid angle dp, which characterises the effective width of the radiation pattern, 
Op = f dnF n (6,(j)) , (5.30) 



in 



the antenna directivity D (9, (ft) 

D (9, 0) ee 1 Fn }°f (f) - = D F n (9, </>) , (5.31) 
^ J 4w dn F n {9, <j>) 

and the maximum directivity Dq 

D = 4^ = 4tt = AnS r (9,4>) \ max 

~ n P j A7T dnF n (9,<p) f 47r dns r (e,<f>) ' l ' J 



we can write (5.28) as [6 



G (9, 4) = r,D Q F n (9, cj>) = r, An . (5.33) 

lip 

Thus an equation for the radiated power per solid angle emerges 



Note that (5.34) depends linearly on the input power Pi. Moreover 77, fip and the polar 
angle dependency F n (9, (j)) are all intrinsic characteristics of the radiating antenna. That 
is, once they have been calibrated, the power output of an antenna is then fully determined 
by the input power Pi feeding the antenna. 

So far we have considered transmitting antennas, and defined the antenna radiation 
pattern and directivity as the characterising angular distribution of transmitted power. 
What about the power gain of an antenna when collecting radiation? The reciprocity 
theorem [60, 70j shows that, under general assumptions (e.g. isotropy of the medium), the 
transmitting and receiving antenna parameters are the same. The directivity and radiation 
pattern characterise the angular cross section of the antenna, and it is the same for both 
transmission and reception of power. In this sense, problems of radiation emission and 
radiation detection by an antenna are interchangeable. 



5.2.2 Antenna effective area 

Power transmission and reception by an antenna may be described, in the far field at least, 
by borrowing from optics the idea of wave diffraction by a geometric aperture A. This 
results in the Fresnel-Kirchhoff diffraction formula [60j, which reduces the computation 
of the antenna far field radiation, given the field at the physical antenna aperture Ap, 
to a scalar diffraction problem with effective geometric antenna aperture Ae- It is found 
that this effective aperture, or effective area, depends on the wavelength of the diffracted 
radiation. Since we would like to apply the same methods to the radiometry of SU(2) 
photons, we must be careful that these results should not depend explicitly on the dispersion 
law. This is a priori not clear and needs to be checked. That the Fresnel-Kirchhoff 
diffraction formula does not depend explicitly on a U(l) linear dispersion law and thus 
holds also for SU(2) photons by implicitly allowing for functional dependencies such as 
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(3.34) is checked in appendix [C| Here we use this result to arrive at an expression for the 
antenna effective area which we need to relate electromagnetic fields in the far field to 
electromagnetic fields in the aperture, to calculate the incoming power at the aperture. 
This is done by first calculating the transmitted power through the aperture, and then 
invoking the reciprocity theorem. We follow somewhat the presentation in [67J but without 
assuming any explicit dispersion law. 

Consider a transmitting antenna, with its physical aperture area Ap, and a monochro- 
matic field distribution E^ of constant phase defined over the surface Ap. This is realised 
by having the Poynting vector of the monochromatic planar wave be normal to the aperture 
surface. We say that the aperture is illuminated by the monochromatic field distribution 
Ek of constant phase. Evaluating the total power <£> crossing through the aperture gives 



$ = / dx a dy a h ■ S k (x a ,y a ) = / dx a dy a Sk (x a , y a ) 
J A J A 

/ dx a dy a |£ k ( Xa j ya ) | ; 
JAp 



1 |k| 

2wk 



(5.35) 



where n is the normal vector to the aperture, and we have expressed the modulus of the 



Poynting vector S k by that of the electric field (compare with (B.64)). The integration 
variables x a and y a parametrize the aperture surface. 



Because of energy conservation, (5.35) must be equal to the power radiated in the far 



field by the antenna, where radial symmetry holds and the monochromatic field propagates 
as an spherical wave after diffraction through the physical aperture 



dfi r S r 



■in 



Hence 



dQ S r 



Itt 



— / dx a dy a \E^(xa-,ya) 
2r 2 w k J A 



(5.36) 



(5.37) 



Take the maximum value Sq of the power density radiated by the antenna. Assume a 
symmetric radiation pattern distribution along the z-axis (for example that of a horn 



antenna), then the maximum can be calculated by using equations (C.18), (C.17) and 
setting 6 = 4> = M 



5 = 5 k (0,0) 



1 Ikrlkl 



8^ 2 -2 



T" W k 



|/»k (o,i 



i |krik| 

Itt 2 r 2 w k 



dx a dy a E k (x 



Inserting both (5.37) and (5.38) in (5.32) gives 

2 

( dXndVn Eh- (Xr, . Vr,) 



D, 



o 



1 iki= 



Ia p dx a dy a Ek ( 

Xa j Da 



vr J dx a dy a \E k (x 

aiVa, 



|2 ' 



(5.38) 



(5.39) 



This is the maximum directivity Dq for a monochromatic field E^ diffracting through the 
physical antenna aperture Ap. For constant illumination, £" k = const, the above equation 
reduces to 



Dn 



7T 



(5.40) 
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Otherwise, from Schwartz's Inequality for complex valued functions follows 
/ dx a dy a E k (x a ,y a ) <A P dx a dy a \E k (x 

J A P J Ap 



This motivates the definition of an effective antenna aperture Ae [S 

2 

J Ap dx a dy a E k ( %a i ya ) 

A E ~ 



j Ap dx a dy a \E k (x 

so that in general Ae < Ap and 
1 



2 ' 



(5.41) 



Do, 



7T 



\KA E . 



(5.42) 



Each monochromatic mode k interacts with the physical antenna as if its aperture had the 
effective area A E = Ae^ and directivity D = Z?o.k-Fn,k (9,<ft), see Eqs. (5.41) and (5.42). 



In the effective Fresnel-Kirchhoff theory, the antenna pattern solid angle also becomes a 
function of wave number k 



n 



(5.43) 



Recall from the reciprocity theorem that this also holds for an antenna receiving power. 

Now let an antenna sit at the origin while collecting monochromatic power of wave 
number k from a source emitting with radiance (with respect to the origin). What is 



the power output Pa of the antenna? From (5.5) and (5.6) we can write the expression for 
spectral power as 



da / dft cos#L k (0, ( 



Before applying this formula to our antenna problem, we must take the area integrated 
over as the effective antenna area A = Ae^- Also, we must replace cos# with r]F n k, since 
the portion of energy passing through the effective area is no longer given by projecting 
the incoming ray to the surface normal, but by the normalized radiation pattern multiplied 



with the radiation efficiency factor. It is r}F nk that projects energy collection, see (5.34). 



Finally we arrive at a formula for the monochromatic power collected by the antenna 



Jn s 



(5.44) 



with £l s the solid angle subtended by the source. Then the total power output Pa received 
by an antenna of radiation efficiency r\ and radiation pattern F n (6, (ft) from a source with 
radiance L k (8, (ft) subtending the solid angle Vl s results from integrating over all available 
modes 



P A = v d|k| P k 



dkA k I dnF n ^{6,(ft) L k (0,, 



(5.45) 
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5.2.3 Antenna temperature 



The concept of antenna temperature was first introduced by R. Dicke in a famous 1946 
article on microwave radiometry [10J, where the power output Pa of an antenna enclosed 
by a perfect black-body cavity was related to the temperature T of the radiation inside the 
enclosure, by equating Pa with the thermal (Nyquist) noise of a matched resistor. Here 
we re-derive the relation between radiation temperature and antenna power output, which 
will be needed to characterise our experiment. 

To calculate the total power output Pa of an antenna completely immersed in thermal 
radiation inside a black cavity, we integrate (5.45) over the full solid angle Vt s = Air and 
each available mode k. The thermal radiance is given by L k bb- Since thermal radiation 
is unpolarized by definition, and antennas are polarized receivers |40| [67] we must divide 



expression (5.45 ) by two, since the antenna can only pick one of the two available transverse 



polarizations. This gives 



P A = ~V J dkA *J 4 dOF„ ik (fl,0)L k 
= X -r) JdkL kbb A k J^ dnF n<k (6, 



bb 



which by (5.30) and (5.43) simplifies as 

p A = ~V [ dkLkbbAkttpk = \rj I dkL k 



bb 



4tt^ 
W 



Changing variables to cancel the squared wavelength factor A 2 = 



(2-) 2 



we obtain finally 



Pa 



1 [a t ^ 

-r? J du k L u bb — ^ 

L0 k 



1 

2tt' 



V / dw k 



(5.46) 



Note that since equation (5.46) is expressed only in terms of the frequency, it is independent 
of any specific dispersion law for photons in the thermal gas. The antenna power output 
Pa is only sensible to the available propagation modes of the thermal gas, which constrain 



the frequency spectrum available for integration in (5.46). 



To write (5.46) in SI units recall again that converting from natural units one has 
uj k —7- hv, (3 — > y^t- ^ e ra ti° V i s unit- less. To have appropriate dimensions, expression 



(5.46 ) must be divided by h. Then the total power Pa collected by an antenna of radiation 



efficiency r\ fully immersed in thermal radiation is written in SI units as 



Pa 



du 



hu 



,hvlk B T _ l 



(5.47) 



In the Ray leigh-, Jeans regime hv ^ fc^T, and we can approximate the exponential 



function in the denominator of (5.46) as e hu / k B T = \ _|_ JkkL. -(-.... In this regime, 



Pa = k B r]T I du . 



(5.48) 
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Let the antenna be connected to a receiver of bandwidth Av = V2 — v\, with v<i > v\ and 
both frequencies inside the Rayleigh- Jeans zone. If Av lies entirely within the available 
cavity frequencies, the power delivered from the antenna then results from integrating 



(5.48) over the whole bandwidth Av 

P A = k B r]T Av . (5.49) 

This simple but important equation confirms what we discussed before: The output power 
at the terminals of an antenna fully immersed in thermal radiation is solely determined 
by the radiation cavity temperature T and accessible bandwidth Av. Let the radiation 
temperature in the cavity be T = Tr. We can then define the temperature Tj\ of an 
antenna placed inside the cavity, collecting radiation over a finite bandwidth Av, as 

T A = r,T R = -^j- . (5.50) 
kbAv 

If Av does not lie within the region of allowed frequencies, there is no power output 
possible, Pa = and the antenna temperature is zero. 



The thermal U(l) gas 

For the U(l) thermal photon gas, propagation over the entire spectrum v > is possible 



(up to geometrical discretization effects which we may ignore, as discussed in section 5.1.4 ) 



Then (5.49) for the U(l) case is just 
ks r]T Av . 



A 



(5.51) 



The thermal SU(2) gas 



As we saw in section |3T5| the frequency spectrum available to the SU(2) photon is bounded 
from below by a cutoff frequency v*(T). Below v*, no propagation is possible. Let the 
detectable antenna bandwidth be Av = V2 — v\ > 0. An antenna completely bathed in 
thermal SU(2) radiation will provide the power output 



p SU(2) 

A 



k B rjTAv ifv 1 >v* 

k B rjT (z/2 — v*) if v\ < v* < v 2 

otherwise . 



(5.52) 



If the bandwidth Av lies entirely below the cutoff frequency v* , the antenna power output 
of a thermal SU(2) photon gas is zer cQ This is in contrast to the U(l) case ( |5.51[ ) which al- 
ways produces a non-zero antenna temperature proportional to the radiation temperature 
in the cavity and bandwidth. The SU(2) gas dynamics manifests themselves as an antenna 
temperature gap below v* This is the key experimental prediction from SU(2) thermo- 
dynamics, as developed in chapter [3j that is accessible by radiometric means. To measure 
the antenna temperature some type of radiometer must be connected to the output of the 
antenna. Such experimental configurations are described in the next section. 



Up to antenna and receiver noise, of course. This is discussed in section 5.3 
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5.3 Radiometer systems 

Before tackling the problem of experimental detection of the antenna temperature gap de- 



scribed in section 5.2.3 we need to familiarise ourselves with basic radiometric techniques. 
The purpose of this section is to give a basic review of the subject, a deeper treatment can 
be found in the literature |67| 163] . First we introduce basic concepts, then we describe in 
brief the two radiometer configurations one could employ to measure SU(2) effects. In this 
section we always take the radiometer bandwidth Au to lie outside the screening region of 
the SU(2) gas and within the Rayleigh- Jeans zone. This is just to simplify our introduc- 
tory review. Discussions about screening effects and their radiometric detection are left to 
section 15.41 



5.3.1 Basic definitions 

A radiometer system is a radiation measurement device composed of an antenna coupled 
to a receiver. The receiver selects a bandwidth Av for detection, and produces a signal 
characterizing the collected radiation. 

In the ideal case, power radiometers produce a voltage signal which is directly propor- 



tional to the power output of the connected antenna, and by way of (5.49) also proportional 
to the antenna temperature 

V = cGP A = cGAuk B T A . (5.53) 
Here G is the amplification factor of the radiometer, and c a calibration factor. A correct 



calibration is of absolute necessity if T4 in (5.53) is to be interpreted as a temperature at 



all. Different calibration methods are discussed in |63l 167] . In the following, we assume 
that the radiometer has been correctly calibrated, and absorb the factor c into G. 

In general, random noise (mostly thermal in nature) intrinsic to the radiometer appa- 
ratus will distort the output voltage ( |5.53 ) , so that it is not a strictly linear function of the 



antenna output power Pa- It can be shown however |63[ 167] that radiometer noise is well 
modelled by considering it as an additive contribution to the measured signal (the antenna 
power output). We thus define the radiometer noise power P/v as the power that would be 
measured by the radiometer viewing a source at zero absolute temperature. 

We can also write Pjy in terms of a temperature, the radiometer noise temperature Tjy = 
Pn/{G k B Av). Thus, radiometer noise may be interpreted as an additional temperature 
T/v collected along with the antenna temperature Ta by an ideal noiseless radiometer. The 
noise power (noise temperature) is an intrinsic characteristic of the radiometer. Allowing 
for thermal noise, the total radiometer output power then reads 

V = G{P A + P N ) = G Av k B (T A + T N ) , 

= GAuk B T a , (5.54) 

where we have defined the radiometer output temperature T Q = V a /(G Av kg). The ra- 
diometer output temperature characterizes the average power output of a radiometer. 
Thermal noise contributes to the stochastic fluctuations in the power output. These are in 
turn characterized by the radiometer sensitivity AT , which gives the minimum detectable 
change in the radiometer temperature. In the case of thermal radiation, both radiome- 
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ter temperature T Q and sensitivity AT a fully characterize the radiometer output]^] In the 
following, we describe two useful radiometer configurations, and give the corresponding 
formulas for radiometer temperature and sensitivity. 

5.3.2 Total power radiometer 

The total power radiometer is the most basic type of radiometer designed to measure 



antenna temperature. Its design is sketched in figure 5.7 An antenna of total temperature 
Ta is coupled via a transmission line to a receiver consisting of four basic stages. The 
first is signal amplification, which is carried out by an amplifier of gain G. This is followed 
by bandwidth selection by a filter of narrow bandwidth Aza In ihedetection stage, the 
input signal is converted to a measurable voltage. For a total power radiometer this is 
achieved by a semiconductor diode working as a square law detector: The diode output 
voltage is proportional to the input signal power Pa + Pn- in the fourth and final stage, 
the voltage signal coming out of the detector is smoothed by integration over a time period 
r to reduce random fluctuations. This integrator effectively acts as a low pass filter, 
suppressing frequencies higher than 1/r, where r is the radiometer integration period. Note 
that the integrator causes the radiometer output to be Gauss- distributed, independently 
of the stochastic distribution of the signal. This can be understood as follows: A single 
integration over a period r of a signal with arbitrary stochastic distribution is equivalent 
to taking the average of many single measurements of said signal. But from the central 
limit theorem, this average will tend to be noramlly distributed. The radiometer output 
after integrations is thus fully characterized by the average output signal and its variance, 
independently of the stochastic distribution of the input signal. 

The noise performance of the radiometer system (antenna plus receiver) can be char- 
acterized the radiometer noise temperature T/v. A useful model for Tjv is given below. 



The total power radiometer as depicted in figure 5.7 is therefore fully characterized by 
specifying a few parameters. While the stages of bandwidth selection and time integration 
are self-explanatory, we believe the subjects of signal amplification and radiometer noise 
temperature merit some further discussion. 



Signal amplification 

Signal amplification in modern radiometers is realized by semiconductor devices. In the 
microwave region (y ~ 1 — 100GHz), High Electron Mobility Transistors (HEMT) are the 
ideal solid state amplifier because of high signal gain and low noise. HEMTs are a form of 
Field Effect Transistors (FETs) so constructed that its charge carriers are confined to a very 
narrow 2-dimensional layer, forming a two-dimensional electron gas (2-DEG). Transport 
properties of this 2-DEG are superior to that of a usual FET since, due to the absence 
of donor impurities, scattering is greatly reduced, thus higher electron mobility and low 
intrinsic noise is achieved [TJ |6D]. This advantages comes at a cost, however. HEMTs at 
cryogenic temperatures suffer from relatively large fluctuations in the amplification factor 
(see e.g. |16j). Let the amplification factor variate in time as G(t) = G + SG(t). Then 
the standard deviation AG = \J (5G(t) 2 ), called gain variation, is the dominant factor 
influencing radiometer sensitivity AT. A clever technique conceived by Dicke |10j avoids 



4 The signal from a radiometer looking at thermal radiation is normal-distributed [67] ■ See also section 
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Figure 5.7: Diagram of a basic Total Power Radiometer: The input signal to the radiometer 
is characterized by the antenna temperature T4, which is determined by the radiation 
collected in the antenna, and the intrinsic radiometer noise temperature T/v, which fully 
characterizes thermal noise in the radiometer. The signal is amplified with gain G, then 
filtered to select a narrow bandwidth Aza The (microwave) signal power is detected by a 
square-law diode. Finally, the signal is integrated over a period r to smooth out random 
fluctuations. 



this difficulty by alternately measuring the signal and some other known reference at a 
higher frequency than the characteristic gain variation frequency of G(t), called the knee 
frequency / knee . 

The defining parameters of an amplifier are thus its gain G, its mean fluctuation or 
gain variation AG and its characteristic frequency /knee- The intrinsic amplifier noise 
contribution to total signal noise will be included in the receiver temperature noise T mc to 
be defined below. 



Noise and attenuation 

Thermal noise due to random motion of electrons in a conductor is the basic source of 
random noise in a radiometer, but not the only one. Other types of noise include quan- 
tum noise, shot noise and flicker noise [67J. As the name suggests, quantum noise is a 
consequence of quantum-mechanical effects in conducting devices and manifests itself as 
discrete random fluctuations in the transported signal. In the Ray leigh- Jeans limit such 
quantum effects can be ignored. Shot noise arises from current flow fluctuations in con- 
ducting devices due to the discrete nature of charge carriers. Flicker noise is caused by 
surface irregularities in semiconductors. All this fluctuations contribute to the noise output 
power, which can be characterised by a noise temperature T/v, as explained above. 

A perfect knowledge of each radiometer component would make it possible to calcu- 
late the random noise contribution of each waveguide, switch, resistor, etc. to the total 
radiometer noise. Instead, we use a very simple model that nevertheless is good enough 
for most purposes [67], where noise is characterized by a few measurable parameters. 



Let an antenna with radiation efficiency rj (see Sec. 5.2.1) be coupled through a 
transmission line of loss factor L to a receiver with characteristic receiver temperature T rec . 
The antenna and transmission line are assumed to be in thermo dynamical equilibrium, at 
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a physical temperature T p . The radiometer noise temperature T^ is then given by [B7] 

T N = (L- n) T pA + LT iec . (5.55) 

The antenna radiation efficiency r/ and loss factor L of the transmission line can always 
be measured. The temperature T PjA is the physical temperature of the antenna and trans- 
mission line coupling at the time of measurement. This is different from the antenna 
temperature Tj\ which characterizes the power collected by the physical antenna aperture. 
The receiver temperature T rec may be defined as the temperature that is measured by the 
receiver when disconnected from the antenna. It characterizes all noise sources internal 
to the receiver. It is strongly dependent on the physical temperature T p ^ rec of the receiver 
at the time of measurement. The value of T rec and its dependenc on Tp^ rec is specific for 
every receiver and must always be measured before any experiment. Here, however, we can 
only estimate T rec based on the values given in the literature for radiometers designed to 
measure radiation in the range of temperatures and frequencies we are interested in. We 
use a simple ad-hoc model 

'^roc( ^p.roc) — ^rec^p,rec ~i~ Crec • (5.56) 

where the receiver noise temperature T rec is a linear function of its physical temperature 
T piec . Parameters a rec and c rec are chosen so that T rec (T Ptiec ) roughly matches the same 
order of magnitude as given in the literature for a receiver working at temperature T pvec . 



This will be explained in more detail in section 5.4 



Total power radiometer output temperature and sensitivity 

As discussed before, the total power radiometer outputs a voltage signal V Q which may be 



interpreted as a radiometer output temperature T , defined in equation (5.54). Here and 
in the following we assume that the radiometer has been properly calibrated so that we 
may always calculate T Q from V a . 

The measured radiometer output temperature T Q of the total power radiometer is then 

T = Ta + T N 

= rjT A + (L-r])T p + LT rcc = r]T A + (L-r])T p + L (a rec T p , rec + c rec ) • (5.57) 

The radiometer sensitivity is determined by the average fluctuation of the output tem- 
perature. We anticipate that this average fluctuation is affected by the high gain variation 
of HEMT-type amplifiers. On the other hand, we expect the integrator to function as a 
low-pass filter and suppress fluctuation frequencies much higher than the inverse of the 
integration period. Both conjectures turn out to be true. Namely, it can be shown that 
the sensitivity of the total power radiometer with bandwidth Ai>, integration period r, 
gain G and gain variation AG is given by the expression [67J 



2 



with T a given as in (5.57). 

Note that both the integration period and gain variation appear as independent terms 
under the square root. This means that large uncertainties in the output temperature due 
to high gain variation AG cannot be compensated just by choosing a higher integration 
time. It is therefore very important to minimise gain variations as much as possible. Both 



expressions (5.57) and (5.58) will later be used to simulate experimental results. 
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Ta T re f 



Figure 5.8: Diagram of a basic Dicke Radiometer: The receiver input is rapidly switched 
between the antenna signal characterized by Ta and a reference signal with temperature 
T re f. The switching occurs at a frequency fg chosen so that it is greater than the charac- 
teristic fluctuation frequency of the amplifier, fg > /knee- Then for each switching period 
the amplifier gain may be taken to be constant, and a signal difference Ta — T re f can 
consistently be defined. See text for details. 



5.3.3 Dicke radiometer 



A basic Dicke Radiometer configuration is shown in figure 5.8 It is based on a clever 



switching technique introduced in 1946 by Dicke |10| that suppresses noise power out- 
put produced due to amplifier instabilities in total power radiometers. Dicke's method 
effectively prevents amplifier gain fluctuations from degrading the sensitivity of the ra- 
diometer |10| . This is realized by switching the receiver input back and forth between the 
antenna signal Ta and a reference signal T re f, at a switching frequency f$ greater than 
the characteristic amplifier frequency /knee- After amplification and detection, the output 
of the square-law diode is multiplied before integration by +1 or —1, depending on the 
position of the Dicke switch. The signal integrator then adds (or subtracts) the power 
coming out of the square-law diode into the total output signal V , which will be given 
by the difference between antenna and reference signals.. An expression for V a can be 
derived as follows: Let the switching frequency /g be large enough so that the amplifier 
gain G(t) = G + AG(t) may be regarded as constant over one switch period 1/fs, further 
let the integration period t be much larger than 1/fs- Then for each switch period, the 
power difference G(t)(PA — -Frcf) = G(t) Avks (Ta — T re f) at an arbitrary time t may be 
consistently defined, since within each period the amplifier gain can be assumed constant, 
G(t)»G(t + l/(2/ 5 )). 

The power difference is integrated over a period r, and the average output voltage after 
integration is given directly by |63] 

V = G(P A + P N ) ~ G(P rci + P N ) = G Av k B (T A - T ref ) = GAuk B T . (5.59) 



The noise temperature T/v cancels in expression (5.59), and the radiometer output tempera- 
ture T Q is just the difference between antenna and reference temperature. The average gain 
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G now multiplies a temperature difference. By choosing T re f sufficiently close to Tjy, we 
can thus suppress the effects of gain fluctuations in the measured signal. Ideally, T4 = T re f 
so that the output signal becomes precisely zero and the dependence on G vanishes. 

Note however that in Eq. (5.59) we have assumed the noise power output of the 
radiometer to be just an additive contribution Pjy to the total power output P Q . It was 
further assumed that the generated noise power is constant over both switch measurements. 
Nonlinear effects in the radiometer manifest themselves as a slight shift from the ideal case 



(5.59). This we discuss below. 



Dicke radiometer output temperature and sensitivity 

Here we give the output temperature and sensitivity formulas for a properly calibrated 
Dicke radiometer. We use (5.55) and (5.56) to model radiometer noise. A refinement of 
Eq. (5.59) must be made to account for nonlinear effects |33| . The radiometer output 
temperature T Q of a Dicke radiometer is then given by the expression 



To = Ta~ T rc f + T,ff . 



(5.60) 



The switching noise temperature T ori results from non-idealities in radiometer components, 
even when comparing between identical signal temperatures [33J. It cannot be included 
within the receiver noise temperature Tv, since T/v cancels in our simple linear noise model. 
As given in |61j . it is of the order of 1CP 1 K. 

The minimum detectable temperature difference of a Dicke radiometer is given by a 



modification of (5.58) to take into account the diminishing influence of gain fluctuations 



AT, 




(T A + T N ) 2 + (T iei + T N y 



+ {-g) {Ta 



-refJ 



(5.61) 



where T/v is defined as in equation (5.55). 
made negligible, by choosing T4 



As was expected, the effect of AG can be 



T re f to be as close to zero as possible. This is the great 
advantage of a Dicke radiometer setup. 



5.4 Radiometric detection of the SU(2) ground state 

We now arrive at the main section of this chapter, where we describe possible ways to 
detect the presence of a thermal ground state. 



5.4.1 Detection principle 

There are two ways in which the problem of experimentally testing the SU(2) photon 
hypothesis can be approached. 

1. Measure the total output power of the photon gas directly at the radiation temper- 



ature Tr and frequency v. From (5.52), a very sharp decline of output power below 



the cutoff frequency v*(Tr) would corroborate the predicted spectral gap, and the 
SU(2) hypothesis. 

2. Measure the output power difference between a prepared U(l) and an SU(2) photon 



gas cavity, both at radiation temperature Tr. From (5.51) and (5.52), this difference 



should be zero above the cutoff frequency u* (Tr) and non-zero below. 
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We saw in chapter [4] that the thermal ground state can be deformed by temperature 
gradients and still remain a good a priori description of the SU(2) gas in the adiabatic 
approximation. However, the thermal radiation formulas derived in this chapter always 
assume a homogeneous temperature distribution. To ensure the validity of these formulas, 
temperature gradients within the cavity should be carefully controlled to be as small as 
possible. 

From section 4.2 we also know that it is possible to deform the SU(2) ground state with 
static electric fields so that the thermodynamical assumptions from which it was derived are 
no longer valid. The emitted radiation then becomes physically indistinguishable from that 
of the U(l) gas. When probing for SU(2) thermodynamics one should be very careful with 
the experimental conditions so as to avoid this. The cavity containing the SU(2) thermal 
gas should be shielded from any stray electromagnetic radiation. For the low temperature 
region of experimental interest, this shielding should guarantee an electromagnetic field 
energy density of at most \E\ = 10 -3 mV/cm. 

On the other hand, the destruction of the SU(2) ground state gives us the possibility 
of generating U(l) behaviour to use as a reference. This will be exploited in the design 
of the experiment. The U(l) gas is realized by applying a sufficiently large static electric 
field. This is explained in the next section. 



5.4.2 Experimental setup 



The basic setup of our experiment is sketched in figure 5.9 
capable of operating in a specified temperature range T, 



Take a cryogenic chamber 
T mM , This range will be 



- mm • • * max 

determined later when the radiometer bandwidth is chosen. The dimensions of the cryo- 
genic chamber should be large enough to accommodate two identical rectangular cavities. 
Let each cavity be of side length Lj , i = 1, 2, 3, and take each Lj to be large enough so 
that modifications due to cavity geometry to the Planck expression for black-body spectral 
energy density (5.14) may be safely ignored. These corrections are given in ( |5.13 ). The 
spectral region of experimental interest is 5 GHz < v < 15 GHz (see below). In this fre- 
quency range, taking each length to be of the order of Lj ~ lm ensures that discretization 
corrections (5.13) are at most of the order 10 -4 |3] and may safely be ignored. 



Label each cavity A and B, as shown in figure 5.9 The walls of A and B are covered with 



a thick enough layer of microwave absorber of the largest possible emissivity e > 0.9999, 
so that radiation within each cavity is completely thermalized. As discussed in chapter 
[4j in order to avoid ground state decoupling effects, A and B must further be shielded 
from stray electromagnetic radiation. We will be working with temperatures of at least 
8K. The ground state energy density at this temperature is equivalent to the electric field 
density of an electric field of magnitude \E\ = 2.216mV/cm (see Eqs. (4.15) and ( 4,16[ )). 
We therefore require the cavity to be shielded from stray electromagnetic radiation down 
to \E\ = 10 _3 mV/cm, to be certain that no decoupling effects occur at all. 

Cavity B is further equipped with two capacitor plates positioned at opposite sides 
within the cavity, connected to an adjustable voltage source V c , so that a constant electric 
field of energy density p c can be generated inside B. 

A Dicke Radiometer is also brought inside the cryogenic chamber. The radiometer 
parameters are discussed in the next sections. Dicke switching occurs between the terminals 
of two identical horn antennas, labelled A and B. Each antenna is placed inside their 
corresponding cavity, as shown in the figure. 
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Figure 5.9: Basic Radiometer Setup for the detection of thermal SU(2) effects. Inside a 
cryogenic chamber at temperature Tr we place two identical cavities A and B, of side 
dimensions L ~ lm. The interior cavity walls of both cavities are fully coated with a thick 
enough layer of microwave absorber with e > 0.99 emissivity. (This is shown in the diagram 
only for one wall). Cavity B is further equipped with capacitor plates opposite to each 
other. They connect to a voltage source V c to generate a static electric field of adjustable 
energy density inside the cavity. Each cavity is also shielded from stray electromagnetic 
fields. Identical horn antennas are placed inside each cavity, both of which connect to the 
radiometer through a Dicke switch. The radiometer is also placed inside the cryogenic 
chamber, so that the whole system: cavities, antennas, couplers and radiometer is at the 
same thermodynamic temperature Tr. 



In order to specify the necessary experimental parameters, we must first choose the 
spectral region and working temperature the radiometer is going to probe. This is done in 
the next section. 



5.4.3 Working frequency and temperatures 

Here we determine the optimal frequency bandwidth for the receiver, and the corresponding 
radiation temperature range. 

In order to test for screening effects, we need to measure the propagation of thermal 



radiation inside and outside the photon screening region (see section 3.5.4). This region is 
determined by a frequency v* = v*(T), above and below of which photon propagation is 
allowed and forbidden, respectively. 

We first define the middle frequency v m and bandwidth Av which will be investigated 
in this experiment. Following [61J . we choose the 10 GHz Low channel which is used there 



with v m = 9.72 GHz and = 860 MHz. Now use the results of section 3.5.4 to plot v* 
as a function of T and select an adequate temperature region. 

Since we would like to avoid boundary effects such as having the frequency v* lie within 
the bandwidth region of the above 10 GHz channel, we avoid temperatures for which the 
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Figure 5.10: The red line plots the cutoff frequency u*, calculated as in section 3.5.4 



The region 500 MHz above and below v* is plotted in light red. The 860 MHz bandwidth 
available to the radiometer is shown as light blue, centred around the frequency v m = 
9.72 GHz. The region lying within available bandwidth and working temperature is shown 
in blue. The temperature regions were so chosen that the measured bandwidth lie at least 
500 MHz above or below the temperature-dependent cutoff frequency v* 



bandwidth boundary would lie within 500 MHz from the cutoff. The temperature ranges 



8K < T R < 14 K and 22 K < T R < 28 K match this condition. In figure [57TO| we have 
plotted the cutoff frequency v* and the region 500 MHz above and below it. The working 
bandwidth of 860 MHz around the middle frequency v m = 9.72 GHz is visible, along with 
the chosen temperature ranges. We therefore need a cryogenic chamber capable of working 
between temperatures T m j n = 8 K and T max = 28 K. 

Note also that with our choice of frequencies and temperatures, we still lie within the 
Rayleigh- Jeans region, since for the lowest temperature, (his m ) / (kBT m i n ) = 0.0583 <C 1. 

The rest of the experimental parameters are taken to be of the order of magnitude as 
those in the radiometer configuration used in [61 J. This is discussed in the next section. 



5.4.4 Radiometer parameters 

In this section, we will simulate the results of radiometer measurements predicted from 



the SU(2) CMB hypothesis described in section 3.4 The most important radiometric con- 



sequence of this hypothesis is the existence of a frequency cutoff v* below which photon 
propagation is impossible. Using the formulas for radiometer temperature and standard 



5.4. RADIOMETRIC DETECTION OF THE SU(2) GROUND STATE 



73 



Radiometer Parameters 




I 


II 


Ilia 


nib 


Antenna efficiency 


V 


0.99 


0.99 


0.99 


0.99 


Middle frequency 




9.72 GHz 


9.72 GHz 


9.72 GHz 


9.72 GHz 


Bandwidth 


tS.u 


0.86 GHz 


0.86 GHz 


0.86 GHz 


0.86 GHz 


Radiometer noise 


T 

-^rec 


0.1T+10K 


0.2T + 15K 


T + 15K 


T + 15K 


Radiometer offset 


^off 


200 mK 


500mK 


IK 


IK 


Gain variation 


(AG/G) 


-30dB=10~ 3 


-20dB=10~ 2 


-10dB=10" 1 


-20dB=10~ 2 


Loss factor 


L 


0.5 dB =10 05 


ldB=10 01 


5dB=10 a5 


5dB=10 05 


Integration time 


T 


0.1s 


Is 


10s 


10s 


Switching frequency 


fs 


75 Hz 


75 Hz 


75 Hz 


75 Hz 



Table 5.1: Radiometer parameters, see text. Note that Njb = 10 log(iV). 



deviation, (5.57), (5.60), (5.58), and (5.61), we can generate expected measurement plots, 
so as to have an intuitive idea of the order of magnitude of both the expected signal average 
and stochastic variation. But before using the radiometer formulas given above, we must 
specify certain characteristic parameters of the radiometer. 

The parameters used to generate simulation plots are given in table |5.1| These are 
taken to be of the same order of magnitude as those described in the literature |2j [33l |331 
l56l l6T| l62| 16^1 169] for radiometers used in experimental situations closely matching our 
own. That is, low temperatures and frequencies in the microwave region. We have divided 
them in three sets. Set I is the parameter set which corresponds to an ideal radiometer 
with above-average sensitive. The parameter set II corresponds to a radiometer with 
average performance. Set III results in a very noisy radiometer. This parameter set is 
further divided in Ilia and Illb, the only difference being the gain variation. The choice 
of parameters is explained in detail below. 

Corrugated horn antennas are a popular choice for microwave applications because 
of high beam symmetry, low side-lobes, low crosspolarization and low reflection over a 
wide broadband. Horn antennas may be designed to satisfy requirements such as working 
frequency, antenna gain and directivity |45| . This in turn specifies the expected antenna 
radiation efficiency, which we take to be 0.99. This sole parameter characterises the antenna 
response in our simplified example. Of course, a full characterisation including radiation 
beam pattern is needed to precisely describe an actual experiment, see for example [62J. 
Antenna efficiency is taken to be the same in all cases. 

Middle frequency and bandwidth of the radiometer were chosen following [61J. Radio- 
metric equipment in this frequency region is commercially readily available, further, the 
temperature region it imposes is such that radiometer and amplifier noise remains low 
enough for our purposes. Higher frequencies show too narrow temperature ranges to test 
for screening effects, compare with figure 5.10 and we would like to measure screening 
effects over a relatively wide range of temperatures. On the other hand, a lower frequency 
would imply increasingly large temperatures ranges to test for non-screening (note the 
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decreasing slope of u*{T) in figure 5.10). Noise and random measurement error increase 



with temperature, so we would like to avoid relatively large temperatures. We find that 
v m =9.72 GHz is a nice compromise. The temperature ranges for testing above and below 
v* are then given by 8K < T R < 14 K and 22 K < T R < 28 K. 

Receiver noise is determined by the thermodynamical temperature of the physical ra- 
diometer system. The HEMT amplifier noise is also dependent on the working frequency 
band [69]. Since we plan to access the SU(2) ground state by probing different temperatures 
Ti, receiver noise will in general change for each temperature. In the actual experiment 
this will become a source of systematic error, since for each run T rec will vary, however 
slightly, as a function of the working temperature. This variation should be calibrated for 
each Ti, so that for every measurement run we may extract the output power temperature 
T Q , whatever the radiometer working temperature. Further, the receiver noise tempera- 
ture T rec crucially limits the experimental precision of the radiometer. Here, we assume 



the simple linear model given in (5.56) 



2~rec(^p,rec) — QTec2p,rec ~\~ C r 



(5.62) 



where T PjTec is the physical receiver temperature. The coefficients a rec and c rec are chosen 
such that in the temperature range we work with, T m i n = 8K and T max = 28 K, they 
are of the same order of magnitude as the receiver noise found in the literature, e.g. 
|2j |3U [56] EU [6H [69], which we have estimated into best, average and worst performance 
(sets I, II and III). 

Radiometer offset is the baseline noise in Dicke switching measurements [33] . We have 
taken it to be of the order of 10 _1 K [61J. 

The power gain variation AG/G is the most important determinant of radiometer 
precision, as can be seen from the equations (5.58) and (5.61). (And as we will later 



confirm in the simulations). Commercially available HEMT amplifiers come with gain 
variation of the order of 10 -3 |5]. We take this value for the best case scenario I. We 
assume a poorer mean fluctuation of the order of 10 -2 for the average case II and "better" 
worst case Illb. For the worst case scenario Ilia we take fluctuations of the order 1CP 1 

Loss factor gives the signal power loss between antenna and receiver. Here we have 
taken to be of the order of 1 dB= 10 ' 1 for the average scenario, based on the |61] setup. 
The values for the I and III cases were taken to be some orders of magnitude above an 
below 1 dB. 

The integration time can be calculated so as to minimise radiometer uncertainty |70] , 
Here we take it to be r = 1 s in the average case and vary the value accordingly for best 
and worst case. Finally, the switching frequency fs when working in Dicke mode should be 
chosen such as to be greater than the characteristic fluctuation frequency of the amplifier, 



fs > /knee- Otherwise the Dicke switching formulas (5.60) and (5.61) do not apply. Here 
we use the value given in [61J, fs = 75 Hz. 



5.4.5 Experimental signature 



Having defined the most important experimental parameters in table 5.1 we may now use 



those values to estimate the radiometer output in the experiment. As already described, 
two experimental procedures are possible, total power mode and Dicke switching mode. 
We simulate an experimental run of 600 measurements. The integration time of the 



radiometer is given in table 5.1 for each set of parameters, so each measurement run has 
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a duration of 1 minute, 10 minutes and 100 minutes for each case I, II and III. Output 
temperatures for each single measurement were simulated by using a Gauss-distributed 
random generator with mean T Q and standard deviation AT Q , both of which were calculated 
using the appropriate formulae (5.57), (5.58) for a total power run and (5.60), (5.61) for a 
Dicke switch run. The radiometer noise temperature T/v given in (5.55) was used in both 
cases to calculate the noise contribution to the output signal. 

The temperature of the cryogenic chamber defines the radiation temperature Tr, since 
the black-boy cavities then have thermodynamical temperature Tr. The physical temper- 
ature of the receiver and of the antenna and transmission-line coupling, were also taken to 
be Tr, since the radiometer is placed inside the cryogenic chamber. That is, T p and T Pirec 
are taken equal to the radiation temperature T p = T Ptiec = Tr. 

Recall that the antenna temperature Ta is by definition the temperature of the collected 
radiation. Inside the black-body cavity, radiation is at thermal equilibrium with the cavity 
walls and Ta = Tr. This means that the antenna efficiency r\ cancels in the expression 
(5.57), since in our case, Ta = T p . 

The results of the simulations are discussed below. 



Total power mode 

In total power mode, we set the cryogenic chamber to a temperature Tr then let the 
radiometer observe the output A. For each Tr we run 600 measurements. The results are 
plotted below. We first describe the experimental procedure, then comment on each case. 

Start the first measurement run at Tr = 28K, then decrease the temperature to Tr = 
26K and measure again, and so on until Tr = 22K. If the radiometer parameters are such 
that the relative signal error is small enough for a distinct linear relationship between T Q 
and Tr to be identified, we may fit the four pairs (Tr, T q ) to a line 

T = A ac T R + C ac . (5.63) 

The subscript ac stands for "'above cutoff "\ We call this measurement series the above- 
cutoff series. 

Decrease now the temperature further, to Tr = 14K, and start the measurement series 
again going down to Tr = 8K. As above, if the measurement errors are small enough, we 
can fit the four measured data points to a line 

T = A bc T R + C bc , (5.64) 

where be stands for "below cutoff". Similarly, we call this measurement series the below- 
cutoff series. 

As an example, take the average radiometer (II) first. It has integration period r =ls, 
so each measurement run takes 10 minutes. The results are plotted in figure [5~lla| for all 
temperatures Tr. The red lines represent simulated experimental runs assuming the SU(2) 
hypothesis, the grey lines would be the result of U(l) dynamics. 

The standard error AT Q is also visible in figure [5~lla| as the width of the measurement. 
For parameter set II, it is of the order of 10 -2 , the dominant contribution coming from the 



average gain fluctuation. This is small enough to fit Eqs. (5.63), (5.64) 



Now consider the above-cutoff series. In this temperature range, the cutoff tempera- 



ture lies below the radiometer frequency bandwidth, compare with figure 5.10 We should 
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therefore measure a radiometer output temperature T Q consistent with a finite Tr temper- 



ature and equation (5.57), independently of whether the photon obeys SU(2) or pure U(l) 
dynamics. Figure 5.11a shows how both predictions overlap in this temperature region. 

Take now the below- cutoff series. Assume SU(2) dynamics. Then, in this temperature 
range, the radiometer is now directly looking at the screening region, compare with figure 



5.10 There can be no propagation inside the cavity within the frequency bandwidth of the 



radiometer. The effective radiation temperature Tr the radiometer sees is therefore zero. 



The expected output temperature inside this region is then given by (5.57) with Tr set to 



zero, Tr = OK. From this same formula, we see that the only output measured is radiometer 



noise. This is now the key experimental fact. If we compare the fitted parameters in (5.64) 



(5.63), for the SU(2) case, the should differ. The parameters above-cutoff characterize a 
finite antenna signal plus noise, while those below-cutojf pick up only the noise signal. For 
small standard variation in T Q both sets of parameters should differ beyond curve fitting 
errors. 

In the case of U(l) dynamics, on the contrary, both linear fits should be similar enough 
so as to have each fitted parameter lie within each others error bars. 



All of this is visible in figure 5.11a The red lines show a slope in the non- screening 
region of about A ac = 1.6, compared to the slope in the screening domain of Aj, c = 0.3. 
The gray lines, on the other hand, both show the same slope A^ c = A ac = 1.6 above and 
below the screening region. 

This is therefore the experimental signature of SU(2) Yang-Mills thermal ground state 
when probed by a total power radiometer. 

We now discuss the other sets of radiometer parameters. The results of the worst case 



radiometer Ilia measuring in total power mode are given in figure 5.12a| We see that in 
this case, the signal error is so great that no definitive linear relationship may be attributed 
between measured radiometer temperature T Q and the cavity radiation temperature Tr. 
It is still of note that an absolute difference is visible in the screening region between 
the U(l) case (grey) and SU(2) case (red), they do not overlap. But this depends on the 
calculation of predicted U(l) signals, which in turns depends on the full characterization of 
the radiometer noise temperature function. This would induce enough error dependencies 
that such a signal would lose any experimental significance. 

Amusingly, this worst case radiometer may be immediately made useful by just re- 
placing the 10 _1 gain variation amplifier with a quieter one of gain variation 10 -2 . Since 
gain fluctuations AG are the dominant error source in the receiver, this change is enough 
to make the Ilia radiometer a good enough Illb. This is seen in figure 5.13a| The high 



receiver noise still influences the absolute value of the output temperature (compare the 
T D ~ 225K at Tr = 28K with T Q ~ 61K at the same radiation temperature for the II 
radiometer), but the different slopes are clearly visible: A ac ~ 4.75 outside the screening 
region and A^ ~ 2.35 inside it. This shows that the most important radiometer parameter 
in our case is is the average gain fluctuation AG. 



Figure 5.14a shows the results from simulating a best case radiometer measurement in 
total power mode. The SU(2) signal is clearly visible, with A ac ~ 0.93 above-cutoff and 
-Abe ~ 0.08 below it. It is remarkable that such a strong signal may be recovered from a 
very simple radiometric technique. This of course assumes the availability of a best case 
radiometer as listed in table 15.11 
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Dicke switching mode 

We now discuss the radiometer working in Dicke switching mode. As explained before, 
in this mode the radiometer switches with frequency fg between measuring the antenna 
temperature T4 and a reference temperature T re f In this experiment we will use a U(l) 



cavity as reference. Such a cavity is prepared as described in section 4.2 

A voltage is applied to the condensator plates inside the shielded cavity B, so that the 
ground state is destroyed, and the thermal radiation behaves identically as in a U(l) gas. 



The field should be of the order of \E\ ~ 1 V/cm or greater, compare equations (4.15) and 



(4.16). This ensures a U(l) radiating cavity. 

The antenna B may now be used as a reference emitter for the Dicke switching ex- 
periment. Each experimental run of 600 individual measurements is done, as before, for 
temperatures between 8K and 14K (the below-cutoff series) and temperatures between 
22K and 28K (the above-cutoff series). 

Since both cavities find themselves inside the cryogenic chamber, both radiate at the 
same temperature Tr, outside the screening region. 



That is, for the above-cutoff series, we expect from formula (5.60) that the output 
temperature of the radiometer be just T g, which is of the order of 0.5 K. This should be 
the case wether the SU(2) CMB hypothesis is true or not, since above the cutoff frequency 
both U(l) and SU(2) cavities have the same power output. Since both signal and reference 
temperature are the same, we further expect the influence of gain fluctuation errors AG 
on radiometer sensitivity to vanish in this temperature region. 

It is inside the screening region, that is for the below-cutoff measurement series, where 
we expect large differences between the U(l) and SU(2) behaviour. In this temperature 
range the cavity B still emits U(l) radiation at a temperature T re f = Tr. If the SU(2) CMB 
hypothesis is correct, however, no photon propagation can occur in cavity A. The effective 



antenna temperature vanishes, Ta = 0, and formula (5.60) now results in T a = T g — Tr 
(in our plots we have redefined T Q to be positive, for convenience). The influence of gain 
fluctuation errors AG cannot vanish anymore, on the contrary, it should return to be 
the dominant fluctuation source. If the SU(2) CMB hypothesis is false, no change from the 
above-cutoff series should be noticeable, since both cavities would still have non vanishing 
output power. 

As an example, take radiometer II with average parameters. The output of a Dicke 



switching measurement run is given in figure 5.11b Outside the screening region, the 
generated values for T Q hover above the radiometer offset noise of T g = 0.5K, as expected. 
Inside the screening region, however, there is a big jump in the signal in the case of SU(2) 
dynamics. This is the experimental signature for the Dicke switching experiment. If photon 
dynamics are indeed described by an effective thermal theory with nontrivial ground state, 
then the signature event would be a nontrivial signal, of the order of 10K below-cutoff. 
For the average radiometer II such a signal is easily detectable. Figure 5.12b| shows 



the generated temperature output for the radiometer with parameter set III. And here is 
where the Dicke switching technique really shines. As long as the switching frequency has 
been chosen correctly (that is, smaller than the characteristic gain fluctuation frequency 
of the amplifier) the SU(2) signal is also clearly visible! The measured signal shows a large 
error inside the screening region but that is not a problem, since we can automatically 
detect the presence of SU(2) dynamics just by having a non trivial signal at all! Compare 
the difference in orders of magnitude from T Q ~ 40K below the cutoff, to IK above it. 
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Such a huge disparity in the Dicke signal output would immediately signal the presence 
of the SU(2) ground state. Since this signal is visible even for the worst of parameters, it 
is the ideal configuration for a radiometric experiment probing for the existence in nature 
of SU(2) physics. On the other hand, if the radiometer output below 14K stayed close to 
zero, the SU(2) CMB hypothesis would be automatically falsified. 

For comparison, diagrams 5.13b and 5.14b show the simulated outputs for a radiometer 
with parameters from the set Illb and I respectively. 



5.5 Bolometry of SU(2) photons 



Aside from the radiometric experiment described in section 5.4 a type of bolometric ex- 
periment can be conceived. Unlike radiometric experiments where a narrow measurement 
frequency band is used, bolometric experiments are by design generally broadband experi- 
ments, probing a wide region in the frequency spectrum. 

Here we review a type of bolometric experiment first discussed in [13j that detects 
the presence of SU(2) ground state dynamics by thermal cooling. The down side of this 
experimental setup is the need for special filters, which should ideally absorb all radiation 
above (or below) the cutoff frequency v* . The experimental setup is as follows. 

Let the apertures of an isolated low-temperature U(l) black-body at temperature T\ 
and that of an SU(2) black-body of identical characteristics at temperature T2 face each 
other, and exchange radiant energy. The U(l) black-body is prepared as before, by means 
of a static electric field. At low temperatures, of the order of 5K-10K, the magnitude of 
the electric field should be of the order of \E\ ~ 1 V/cm to guarantee a decoupling from the 



SU(2) ground state (compare with section 4.2). Couple the SU(2) black-body to a large 



heat reservoir, so as to keep its thermodynamical wall temperature T2 constant, and insert 
a low-pass filter within the common aperture so that photons above the cutoff frequency of 
the SU(2) cavity v*(T2) are completely absorbed. The filter works bi-directionally, that is, 
no matter from which direction the photons come from, they are absorbed if they propagate 
above v*{T2)- 

Our bolometric experiment is then so configured that, if the SU(2) CMB hypothesis is 
correct, power is transmitted unidirectionally from the U(l) cavity to the SU(2) cavity. 
Because of the spectral gap below is*, the SU(2) cavity cannot radiate in this region of 
the spectrum. Energy is transfered continuously from the U(l) cavity to the SU(2) cavity 
only, and thermal equilibrium can never be reached. This effect can then be detected as 
the lowering of the U(l) cavity temperature T\. 

Consider now the case T\ = T2 = T. Replace the low-pass filter from before with 
a high-pass filter which absorbs all radiation below the cutoff frequency u*(T). Let a 
monochromatic radiant beam of frequency v > u*(T) propagate from the SU(2) cavity 
into the U(l) region trough the aperture, along an axis h with solid angle df2 u(1 ) defined 



on a sphere centered at a point within the common aperture (recall Fig. 5.1). Similarly, 
have another radiant beam of same frequency v propagate from the U(l) cavity into the 
SU(2) region, along an axis — h with solid angle dS7 SU ( 2 ). 

First note that by Snell's law |39| . we can relate both solid angles as 



^^SU(2) 

ff$7 u(1) 



(5.65) 



5.5. BOLOMETRY OF SU(2) PHOTONS 



79 



Total Power Radiometry 
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(a) Output temperature T D measured by a total power radiometer. The graphic was generated 
by using a Gauss-distributed random generator of mean T a and standard deviation AT a , which 
were calculated using the corresponding formulas ( 5.57[ ) and ( |5.58[ ) 



Dicke Switching Radiometry 
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(b) Output temperature T measured by a Dicke switching radiometer. The graphic was generated 
by using a Gauss-distributed random generator of mean T and standard deviation AT , which 



were calculated using the corresponding formulas (5.601 and (5.611 



Figure 5.11: Simulated results for the set of radiometer parameters II given in table 
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(a) Output temperature T a measured by a total power radiometer. The graphic was generated 
by using a Gauss-distributed random generator of mean T and standard deviation AT D , which 
were calculated using the corresponding formulas ( 5.57[ ) and ( |5.58[ ) 
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(b) Output temperature T measured by a Dicke switching radiometer. The graphic was generated 
by using a Gauss-distributed random generator of mean T a and standard deviation AT D , which 



were calculated using the corresponding formulas (5.601 and (5.611 



;ure 5.12: Simulated results for the set of radiometer parameters Ilia given in table 
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Total Power Radiometry 
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(a) Output temperature T measured by a total power radiometer. The graphic was generated 
by using a Gauss-distributed random generator of mean T and standard deviation AT , which 
were calculated using the corresponding formulas ( 5.57[ ) and ( |5.58[ ) 
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(b) Output temperature T measured by a Dicke switching radiometer. The graphic was generated 
by using a Gauss-distributed random generator of mean T and standard deviation AT , which 
were calculated using the corresponding formulas (5.601 and (5.61 1 



Figure 5.13: Simulated results for the set of radiometer parameters Illb given in table 
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(a) Output temperature T measured by a total power radiometer. The graphic was generated 
by using a Gauss-distributed random generator of mean T a and standard deviation AT D , which 
were calculated using the corresponding formulas ( 5.57[ ) and ( |5.58[ ) 
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(b) Output temperature T a measured by a Dicke switching radiometer. The graphic was generated 
by using a Gauss-distributed random generator of mean T Q and standard deviation AT , which 



were calculated using the corresponding formulas (5.601 and (5.611 



Figure 5.14: Simulated results for the set of radiometer parameters I given in table 
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where the phase velocity v p h of energy propagation inside the SU(2) black-body is given 
by v p h = ^ = c 1 — (f]7p) ■ But this results in the identity 

Lf fe < 2 V,T)dft su(2) = L u v §(v,T)dn m) . (5.66) 

we recognize the left and right hand side of expression ( |5.66 ) 
as the total power emitted respectively by the SU(2) and U(l) cavities. That is, for 
temperature T and frequency v > v*(T\ the total power transfer from the U(l) cavity 
into the SU(2) cavity is equal to the total power transfer from the SU(2) cavity into the U(l) 
cavity. Above the cutoff frequency v*(T), the system is in equilibrium, and no temperature 
change takes place. 



From Eqs. (|5.5|) and (5.6 



Chapter 6 

Summary and Outlook 



The main subject of this thesis was the investigation of possible experimental methods to 
detect the presence of a nontrivial effective thermal ground state, under the hypothesis that 
the abelian U(l) gauge group responsible for photon dynamics in the Standard Model is 
not fundamental, but emerges from a pure SU(2) Yang-Mills theory as the unbroken gauge 
symmetry at finite temperature in the effective theory of the deconfining phase after coarse- 
graining. In this thesis, we refer to this identification as the SU(2) CMB hypothesis, following 
the literature. The CMB in the subscript refers to the fact that the single parameter of 
the nonperturbative approach to SU(2) thermodynamics can be fixed by the observation 
that the temperature of the CMB radiation gas filling the universe should be very close to 
the critical temperature for the deconfining/preconfining transition of the thermal SU(2) 
gas in its deconfining phase. 

Under the SU(2) CMB hypothesis, thermalized photons are understood as free propagat- 
ing thermal quasiparticles whose dispersion law in the effective theory is determined by 
interactions with massive vector modes. 

In this thesis, we have calculated these dispersion laws and discussed the results, in the 
context of possible experimental detection. It was found that the transverse polarization 
modes propagate via a modified dispersion relation. The modified dispersion relation 
reflects the fact that a temperature- and momentum-dependent dynamical photon mass 
a/ G(p, T) emerges in the effective theory. Further, this dispersion law shows a temperature 
dependent minimal propagation frequency, which we call the cutoff frequency v* (T) . It 
is shown that below this frequency, no propagation of transverse polarized photons is 
possible. This defines a spectral screening region for frequencies below v*{T). We have 
calculated numerically a functional relationship for u* (T) by fitting to a power law. A novel 
result of this thesis is the calculation of the dispersion law for the longitudinal mode. The 
dispersion relation shows that longitudinal mode propagation is possible only for very short 
momentum (or long wavelength) and therefore cannot be probed by neither radiometric 
nor bolometric means. 

The prediction of a temperature-dependent screening region has important experimen- 
tal consequences, which were investigated in this thesis. It was shown that an antenna 
immersed in an SU(2) photon gas of temperature T collecting radiation within a band- 
width Av = V2 — v\ > about a middle frequency v m = v\ + ^ would measure a power 
output P Q identical to that of the U(l) gas, as long as the entire bandwidth Av lies above 
the cutoff, v\ > v*(T). By decreasing the temperature to T' < T, however, an experi- 
mental situation can be achieved in which the entire bandwidth Av lies below the cutoff, 
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v<i < i>*(T'). In this case, photon propagation is forbidden and the antenna would deliver 
zero power output. This fundamental difference between U(l) and SU(2) behaviour be- 
comes the experimental signature for the existence of SU(2) photon dynamics. Namely, 
since this region can be probed by radiometrical or bolometrical methods, the SU(2) CMB 
hypothesis becomes experimentally testable, and thus falsifiable. This is the main result 
of our thesis. 

Two experimental techniques proposed in this thesis for probing the power output 
of a thermal photon gas at low temperatures and frequencies are total power and Dicke 
switching radiometry. It is shown that both can be made sensitive enough to detect SU(2) 
ground-state effects on photon propagation. Dicke switching, however, is preferable because 
by design it suppresses gain variations of the amplifier. 

A bolometric technique is also discussed in which the observable effect is cooling by 
power emission. This technique depends on the availability of a low-pass filter which should 
bidirectionally absorb radiation above the cutoff frequency v*(T). 

Beside our main experimental result, we offer a novel theoretical investigation on the 
possibility of studying thermally inhomogeneous ground-state dynamics. This subject is 
by itself interesting (thermomagnetic effect), but also necessary to elucidate questions of 
experimental parameter tolerance. We have found a solution to the pure SU(2) Yang- 
Mills equations of motion in Minkowski signature with an adjoint scalar field, which in 
the context of nonperturbative Yang-Mills thermodynamics and under the adiabatic ap- 
proximation allows for the generalization of the thermal ground state to space-dependent 
temperature profiles. It was shown that such temperature gradients are capable of gener- 
ating magnetic fields. This is the SU(2) thermomagnetic effect. 

Finally, an investigation on the consistency of 1-loop resummation of the energy density 
and pressure of the SU(2) gas was made. We found that thermodynamical consistency 
cannot be guaranteed if radiative corrections to the vector modes in the effective theory 
are ignored. 

It is remarkable that the experimental detection of SU(2) effects and thus confirmation 
of the SU(2) CMB hypothesis is achievable by a simple tabletop experiment consisting of a 
cryogenic chamber, black-body cavities and a well calibrated radiometer. The experimental 
realization of the black-body experiment described in this thesis is a worthwhile endeavour 
because of its considerable implications were the SU(2) CMB hypothesis be confirmed. 

We introduced this work by referring to the low-frequency power-like deviation of the 
CMB spectrum from perfect black-body physics, which can be explained by assuming the 
SU(2) CMB hypothesis. Another important consequence of SU(2) dynamics can be read in 
|58j . There it is shown that the realization in nature of the SU(2) CMB hypothesis would 
invalidate the Higgs mechanism for electroweak symmetry breaking in the Standard Model. 
SU(2) CMB yields six additional degrees of freedom to the total number of cosmologically 
active relativistic degrees of freedom g* at Big-Bang nucleosynthesis. This in turn would 
imply a value for the Fermi coupling at zero temperature Gf of about 12% larger than 
predicted by electroweak SM physics. 

It is our hope that this thesis aids further experimental investigations regarding the 
SU ( 2 )cmb hypothesis. 
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Appendix A 



SU(2) Equations of Motion with 
Scalar Field 



A.l Introdution 

Here we derive a special solution to the SU(2) Yang-Mills Equations of Motion with a scalar 
field <I> taken as perturation. The derivation is subject to some simplifying assumptions. 
Our only purpose here is to present the mathematical details, the physical interpretation 
of the result is discussed in chapter [4] 

We work with Minkowski metric, defined here as g^ u = diag( — 1, 1, 1, 1) to streamline 
calculations 



A. 2 A particular solution 



We now embark in our search for particular solutions of the SU(2) Yang-Mills equations 



of motion (3.9) with scalar field. These are given by 



DuC" = ie [$,£>"$] 



(A.l) 



We write the scalar field in capital greek <3? so as to differentiate it from the inert scalar 
<p described in chapter [3j since we will allow the scalar field to variate in space. Here we 
work in the physical unitary gauge (3.14), <I> a = 2<5 a3 |<I ) | 

Recall now that both $ and the trivial topology gauge field = a^t a live in the adjoint 



representation, so the covariant derivative acts on them as = 



-»e[a M , -J 



With this 



in mind we expand (A.l) 



-2ied^ [a M , a„] + ie [d^a^, a u ] - ie [<%a^, o M ] - e 2 [a M , [a M , a v }} 



[a 

ie +e 2 [$,[*,a v ]]=0. 



(A.2) 
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A. 2.1 The Yang-Mills equations for static fields in Minkowski space- 
time 

We are going to look for stationary solutions, that is, 8q = 0. The time-like v = and 



space-like v = k components of (A. 2) then become 

djOo — 2iedj [aj, ao] + ie [djaj, ao] 
-e 2 [a M , [a„, a ]] + e 2 a ]] = . 



(A.3) 



5|afc — dkdjdj — 2iedj [ctj, a&] + ie [djdj, afc] 
-ie [d fc a M , a M ] - e 2 [a M , [a M , a fc ]] + e 2 [$, [$, a fc ]] = . 

To simplify our notation, we define the complex fields 



Now write equations (A.3) and (A. 4) explicitely in color indices 

„2|^.|21 



(A.4) 



(A.5) 



[d 2 + 2iea 3 jdj + iedja* - e 2 a 3 l a 3 l + 4e 2 |$| 



= ie 



+ e 2 [ 



ctftaC — a n a,', — anor, I a 



oj o j 

(A.6) 



[d 2 - 2e 2 a+a^] al + e 2 [a^a^ + a 



+ c.c. , 



(A.7) 



-ie 



[d 2 + 2iea 3 j d j + iedja 3 - e 2 a 3 l a 3 l + 4e 2 |$| 2 ] a+ 
[a|9j + 29j-a|] a+ + ie [a^dk - d k o?J\ a+ - d k djaj 



e2 [ a t% ~ a k a t~ a l a V\ a t > 



(A.8) 



[d 2 - 2e 2 a+a M ] af + e 2 [a+a^ + a fc a+] a 3 - d^-af = -ie [a k dj + 2dja k ] a} - iea^ i 

+ c.c. (A.9) 

The symbol c.c. denotes the complex conjugate of the expression on the right hand side of 



equations A.7 and (A.9). 



A. 2. 2 Simple configurations with U(l) magnetic field 

In order to solve the above equation system, we make further simplifying assumptions. 
Namely, that the time-like-component of the tree-level massless gauge field (photon), as 
well as the space-like components of the complex tree-level heavy fields ( A,5[ ), both vanish. 



«o = 4 = • 



(A.10) 



Since we have already assumed stationarity, this restriction limits the set of possible so- 
lutions to pure magnetic field configurations. Under our assumptions, the electric field is 
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always zero E = -V.4 - d t A = 0, and only magnetic configurations B = V x A are 
possible. From (A. 10) it follows that the expression a^a^ = vanishes, equations (A. 7) 
and (A. 8) are thus trivially satisfied. Inserting (A. 10) in equation (A. 6) we obtain (recall 
our metric convention) 



[d 2 + 2iea 3 j d j + iefyaj - f?a)a) + 4e 2 |$| 2 ] = . 



(A.ll) 



Recall the identity fd 2 f = \d 2 f - (d/) 2 , which we use to write the above expression as 



V + e 2 (4|d>| 2 - a JaJ) (flo+) 2 = {QA y _ ied . \ {4) 



2 3 
a 3 



Reordering terms, (A. 12) becomes 



(A.12) 



e 2 (4|$| 2 - a)af) 



a, 



o 



+ ^9 2 Qq 2 /oq 2 + iedj 



+ 2 3* 







0. 



Expanding in a positive real magnitude ao(x) and a real phase (f(x) 
ag = a e^ , 



we can simplify the gradient, Laplacian and divergence terms in (A. 13) 



0, 



+ 2 3' 



/a f = dia z + 2a) d i4- 



djCiQ djao 



«o 



+ idj<p , 



O^af/af = 2 (&!> 
a 



+ 2 



d 2 a 
a 



4 (c^f + 2i d z <p + 4djp 



djCLQ 

a 



Separate (A. 13) into real and imaginary components 

d 2 a 



e 2 (4|$| 2 - a]a)) - (d^f - 2ea)d^ + 



+i 



a 



ao 



Of) 



0. 



Taking Re[( A. 15 )]=0 and Im[(A.15)]=0 we obtain two real equations 
d 2 a 



e 2 \a 3 \ 2 



a 



+ djcp [djip + 2eaf\ = 4e 2 |$| 2 , 



ec^-a? + 2eaj— h a V + 2dj(p— — 



0. 



«o 



«o 



Inserting flA.10) in (A. 9) gives us 



[d 2 + 2e 2 a f ] a\ - d k d_ 



, j0 | =ie [a c^ag - <9 fc a 



Writing Og in components as in (A. 14), we finally get 



d 2 + 2e 2 (a ) 2 a\ - d k d ja 3 



2e (a ) 2 d k tp. 



(A.13) 



(A.14) 



(A.15) 

(A.16) 
(A.17) 

(A.18) 

(A.19) 
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Reordering (A. 19) 



a k 



2e (a ) 2 [ea\ + d k ip] = d k djdj - d 
Now, from ( A.17| ) one may easily see that 



(A.20) 



0, 



(a ) 2 (ea 3 j + dj<p) 



0. 



That is, there is some vector field V(x) such that 

(ao) 2 (ea 3 k + d k ip) = rot k V , 

where we have defined rot k V = (V x y)fc. 
Moreover, from (A.20) we infer 

2e rot k V = 8*8 ja) - d 2 a\ . 
So V(x) must be of the form 



(A.17) 



2eV k = rotfc a 



(A.20 ) 



(A.21) 



Equation (A.21) can be used to rewrite (A. 16) as 

d 2 a 



4e 2 |$| 2 = e 2 |a 3 | 2 



o 



e(a ) 



(d j(p y 



(A.16) 



A. 2. 3 One-dimensional temperature distribution 



We now specialize to the case of spatial dependency in just one dimension, and assume that 
the scalar source $ and all $-coupled physical fields to be functions of coordinate z = X3 
only, so that |$(x)| = |$(z)|, and so on. Our double primed equations will hold for this 
special case only (together with stationarity and unitary gauge fixing (3.14), remember). 



Under this assumptions, we can write equation (A.21) in the form 

a\{z) , 



4 



•2<rot k V = rot fc B = (5 k3 - 1) d l z a% . (A.22) 

The magnetic field B is perpendicular to the z-axis. Further, together with equation 
(A.17 ), equation (A.22) directly implies 

eal(x) = -d z tp(x) . (A. 23) 

It follows that <p(x) functions as a U(l) gauge parameter, and can be set to zero (together 
with aj). The identity (A. 23) allows us to use equation (A.17 ) to solve for ao 



d 2 A 
2e 2 a\ 



&14 

'2e 2 a 3 2 ' 



(A.17') 



since under our assumptions d\tp = d%ip = 0. Equation (A.17'1 implies the identity 
a\d"la\ = a^Sfal, and can be used to determine ao from the fields a\ 2 - These fields also 
define the magnetic field B, B\ = — (? z a|, B2 = d z a\. 
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The expression (A. 16 ) for |$| 2 also simplifies in this special case. From (A. 22) we read 



rot3B = 0, which together with di<p = c^V 9 = and equation (A. 23), reduces (A. 16 

d 2 z a 



to 



4 4e2fl o 



(A.24) 



The temperature profile can be calculated from (3.2). Inserting (A.24) we obtain 

(r^/r 1 ) 



T(z) 



A 3 



2vr|$| 



A 3 /l 

2^ 



3 3 _ 
A ^ a ^ 4e 2 a 
i=l,2 



E 



(A.25) 



This means that if we find a reasonable Ansatz for a 3 (z), we automatically obtain the 
fields ao, B, |<3?| and the profile for the temperature T(z). Such solutions and their physical 
interpretations are discussed in chapter (fil. 



Appendix B 

Review of Classical Gauge Fields 



The basic result of the nonperturbative approach to SU(2) Yang-Mills thermodynamics 
in the deconfining phase reviewed in chapter [3] is the emergence of an inert scalar field (j), 
which together with a pure gauge configuration a^ s (see equation (3.10)) defines an effective 



thermal ground state acting as a background for the dynamics of weakly interacting thermal 
quasiparticles. The scalar field <p is also responsible for breaking the original SU(2) gauge 
symmetry down to its abelian subgroup U(l). Two of the original off-Cartan gauge modes 
acquire a temperature dependent mass, while a third one is left massless at tree-level. Under 
the SU(2) CMB hypothesis, this tree-level massless (TLM) gauge mode is identified with the 
SM photon, (see section 3.4). Radiative corrections, however, induce an effective screening 
mass m = 



G for the transverse polarizations of the TLM mode, which modifies their 



dispersion law (see section 3.5.3). Further, a possible longitudinal polarization is rendered 
nonpropagating (see section 3.5.5). 

Massive particles of spin 1 are described by the Proca equations. But the standard 
radiometric results used in this work where derived for massless spin 1 fields, described by 
the Maxwell equations. We must therefore explicitly check how these results are altered 
for spin 1 (vector boson) particles of nonzero dynamical mass. 

In this Appendix, we review the theory of classical spin-1 (vector boson) fields, following 
the presentation in [20J. First we introduce the Lagrangians for the massless and the 
massive case, i.e., the Maxwell and Proca fields. The respective equations of motion are 
derived, along with expressions for the energy and momentum density for each field. In 
section |B.2| polarization vectors for the massless and the massive case are introduced, we 



then expand the energy and momentum density in Fourier modes. Section B.3 recalls the 



continuity equation. Finally, since propagating excitations of the SU(2) TLM mode may 
only be transverse, transversality must therefore be imposed on these Proca solutions. This 
is done in section IB. 41 

The results of this appendix are mostly needed to rederive Planck's black-body formula 
(see section 5.1.4) and the Fresnel-Kirchhoff diffraction integral used in calculating the 



effective antenna area (see appendix |C| and section 5.2.2 ). It will be seen that the standard 



results remain unchanged, save for the modified dispersion relation 



(hw p ) 2 = (cp) 2 + G(^,p;T,A) . 



(B.l) 
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B.l Lagrangian formulation 
B.l.l The Maxwell field 

As we saw in chapter [2] the dynamics of the massless spin-1 field A^ subject to an external 
current density j M = (p, j) are described by the Lagrangian density 



M 



1 



(B.2) 



with the gauge field A^ = (Aq, A), the field strength F^ v = d^A u — d u A^ and the coupling 
e is the charge of the external field inducing the current density j^. The physical electric 
and magnetic fields are defined as in (2.4) and (2.5). 
By the principle of least action 



5S = 5 gTxC 



d„F^ 







6C 



5C 



SdpAv 



we obtain the Maxwell field equations of motion 

f, 

which for the massless gauge field A^ read 
UA V + d u d ll A^ = f . 



Taking the divergence of (B.3) leads to the continuity equation 
duf = 0, 



(B.3) 



(B.4) 



(B.5) 



since F^ u is antisymmetric. Thus, in the massless spin-1 field theory, current conservation 
is automatically guaranteed. This is a direct consequence of the conservation of U(l) gauge 
symmetry. 

Writing (B.4) in terms of the electromagnetic fields (2.4) and (2.5), we recover the 
familiar inhomogeneous Maxwell equations, which we write here in SI units 

V-E 



P_ 



V x B - eonodtE = p j ■ 



(B.6a) 
(B.6b) 



Recall also from chapter [2] that the Bianchi identity satisfied by the field strength tensor 

plCEJErj 



0, 



results in the homogeneous Maxwell equations, here given in SI units 

V • B = 0, 
V x E + d t B = . 



(B.7) 



(B.8a) 
(B.8b) 



In the following, we will be mostly interested in vacuum configurations, that is, j M = 0. 



B.l. LAGRANGIAN FORMULATION 
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Note first that the covariant formulation in (B.2) introduces redundant degrees of free- 



dom that allow for local gauge transformations without affecting the Lagrangian density 



C 



A' 



„ Afj, + d^A(x) 



(B.9) 
(B.10) 



for a Lorentz scalar function A(x). In order to interpret the gauge field A^ physically, we 
must fix the degrees of freedom so as to eliminate spurious polarization modes. This is 
called gauge fixing, and can be done in various ways. A manifestly covariant gauge fixing 
is given by the Lorentz gauge 



dpA" = 0, 



In this gauge, the field equations for the Maxwell Field (B.4) become 



OA* 



(B.11) 



(B.12) 



The Lorentz condition (B.ll) however still does not fix the gauge completely. Gauge 



equivalent field configurations may still be found that are related to each other by gauge 
transformation functions A(x) which solve the Laplace equation, OA(x) = 0. The Lorentz 
constraint d^A 11 = cannot fix A^ uniquely, as the Maxwell field possesses just two 
independent degrees of freedom, not four. 



Energy and momentum density 

The familiar expression for the energy-momentum tensor 
5C 



8^ 



Sd,A a 



8 v A c -g^C, 



(B.13) 



turns out to be not symmetric (also not gauge invariant) when the Lagrangian B.2 



is 



inserted directly in (B.13). To construct a physical energy-momentum tensor T^ v we add 
a surface term to (B.13) 



1 



g ^ F ^ F + F \F™ + g» v j c A° - fAv . 



(B.14) 



Writing (B.14) in terms of the physical fields (2.4) and (2.5), and taking = 0, we recover 
the familiar expression for the electromagnetic energy density in vacuum 



-i00 



1 



770-O 



i y 2F Foi + F lJ Fij) + F C F' 



-E 2 + 



SjJ cfc _ cfc « \ Tp T?rr 
°m u n °m u n I r jk r 



E 2 + 



: ijkp 



jk 



Fmr, 



(E 2 + B 2 ) , 



(B.15) 
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with e ij Ei mn = 5m^n — dm$h- The Poynting vector of the Maxwell field in vacuum, that 
is, the vacuum momentum density, can also be easily calculated 

S k _ T 0k = F 0* F(7 k = E ip.k =E il ^ § m § k _ 6 n 6 k^ 

= e ijk E* , 
S = ExB. (B.16) 
In SI units the fields are rescaled as E — >• -, S — » -, T^ v -> -T F and 

p = Ue E 2 + -B 2 ) , (B.17) 



2 V /V 

S = — E x B . (B.18) 
Ho 

This basic results can now be compared with the same expressions for the Proca field, 
which we review next. 



B.1.2 The Proca field 



The massive spin-1 field, or Proca field, is described by the Lagrangian density (B.2) plus 
a mass term quadratic in the gauge field 



C = -\f^ v - j^A" + l -m 2 A^ . 



(B.19) 



The electric and magnetic fields are again defined as in (2.4) and (2.5). 

By the principle of least action we obtain the Proca field equations of motion 

,2 A u 



V" + m A A v 



J 



In terms of the field A^ this reads 
(□ + m 2 ) A v = d v d^ + f 



(B.20) 



(B.21) 



The inhomogeneous Proca equations result from writing (B.20) in terms of the electric and 
magnetic fields, as in the Maxwell case. In SI units, they are given by 



V -E 





/mc\ 2 








/mc\ * 


Vo3 - 





A Q 



(B.22a) 
(B.22b) 



These are the inhomogeneous Maxwell equations (B.3) modified by a mass term ( ^) ■ 
Note the appearance of h along with c. As in he massless case, the Bianchi identity ( |B,7 ) 
results in the homogeneous Proca equations 



V-B = 0, 
V x E + d t B = . 



(B.23a) 
(B.23b) 



B.2. PLANE WAVE EXPANSION AND DISPERSION RELATIONS 
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which remain unchanged from the massless case 

The addition of a mass term to the spin-1 field Lagrangian (B.2) breaks the U(l) gauge 
invariance. As consequence of this, the current j v is no longer inevitably conserved. Taking 
the four-divergence from (B.20), we get 

m 2 d v A u = d v f . (B.24) 

This expression does not automatically vanish, in contrast to the massless case, and one 
needs to explicitly postulate d u j u = if current conservation is to be maintained. By 
demanding current conservation, the massive spin-1 field automatically satisfies the Lorentz 
condition 

d u A u = 0. (B.25) 
There is no gauge freedom left as in the Maxwell case, and the Proca field equations become 
(□ + m 2 ) A 1 * = f . (B.26) 

Energy and momentum density 



The energy momentum tensor is given by (B.14) with an added mass term 



T nv = \g^F a ^F a p + F^ a F uv + g^ u j a A u - fA h 



4 - - -c0 

+ to 2 (a^A" - ^g» u A a AA . (B.27) 
Setting j 11 = we obtain the energy density in vacuum 

p = T 00 = I (E 2 + B 2 ) + \m 2 (Al + A 2 ) . (B.28) 
Similarly, the Poynting vector of the Proca field in vacuum can be calculated to be 

S = ExB + m 2 A A. (B.29) 

By converting to SI units, the fields are rescaled as S - >• -, E — > A — > ^r, T^ v -> 
j-T^ u . The energy and momentum density of the Proca field are given in SI units by 

"4(^ 2+ ^ B2 ) + ^(f) 2 H + ^ A2 )' (Mo) 

S = ^ExB+(^) 2 i Aj. (B.31) 

B.2 Plane wave expansion and dispersion relations 
B.2.1 Polarization vectors 

To ease calculations, we will expand the vector field A 11 in its Fourier modes, 

A k ^(x) = N k e k ^e~ lkx , (B.32) 
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with ej-x^ a set of four 4-dimensional polarization vectors labeled by A. Demand that the 
polarization vectors e^x^ 1 build a 4-dimensional orthonormal system 



9w 



(B.33) 



Following |20| . we may define a set of polarization vectors obeying (B.33) for the Maxwell 
and Proca fields. 



Polarization vectors for the massless field 

In the Lorentz frame where the plane wave has 3-momentum k, we may pick two space-like 
transverse polarization vectors 



= (0, £fe,i) , 

Cfc,2 = (0, €fc,2) , 

such that for i,j G {1,2} it holds 
k • e k ,i = , 



(B.34a) 
(B.34b) 

(B.35a) 
(B.35b) 



For the longitudinal polarization e k fi a covariant expression can be written by first defining 
the time-like unit vector n with n 2 = 1 and n = (1,0) in the local frame. Then 



k — n(k ■ n) 



(k • n) 



k 2 



In the local frame therefore 
be taken to be n, 



(B.36) 



(0, ). The remaining time-like polarization degree can 



e kfi = n. (B.37) 
One can easily check that this set of four polarization vectors satisfy the normalization 



condition (B.33) 



Polarization vectors for the massive field 

Again, working in the Lorentz frame where the plane wave has 3-momentum k, we define 



transverse polarization vectors as in the massless case (B.34), (B.35). The third vector 



may be chosen so that its spatial components are parallel to k. Its zero-component can 
be further specified by the condition that 3 is orthogonal to the momentum four-vector 
&/i e fc,3 M = 0. This completely defines the longitudinal polarization vector 

/Ikl k k°\ , . 

ek >*= 'TUT- ' B - 38 
\m \k\ m J 

Note that this definition was impossible in the massless case, where 4-dimensional transver- 
sality would have implied infinite norm. 

The fourth time-like polarization vector can be defined as 



e k,0 



I*. 

m 



(B.39) 



This set of polarization vectors also satisfy condition (B.33). 
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B.2.2 Plane wave expansion of the Maxwell field 



In the Lorentz gauge d^A^ = the equations of motion (B.4) for the Maxwell field A^ 
become (B.12) and the general solution may be given as a sum over polarizations A, which 
labels the set of four polarization vectors for the massless case given in (B.34), (B.35), 
flB~37l ), flB~36l ). 



A k 



fc°=w k 



(B.40) 



The familiar dispersion relation for the Maxwell field in vacuum can easily be recovered in 
the Lorentz gauge d^Aji = by inserting (B.40) into (B.12) and taking = 0, 



(B.41) 



Consider now the Coulomb gauge V • A = 0. Then (B.4) becomes 

□A + d t VA° = , 
V 2 A° = . 

We can take ^4° = and expand A in its Fourier modes. Taking the limit 

d 3 k 



lim ~ -)• / 
V->oc V ^ J 



(2vrr 



(B.42a) 
(B.42b) 



(B.43) 



we can write the sum over all modes k as 
d 3 k 



A(x) 



A 



(2k)* 



AT I P „ — ikx | -* ikx 

Nk [ek,\e + efc,A e 



k°= 



(B.44) 



Since V • A = we sum only over the transverse (physical) modes A = 1,2. Choose linear 
polarization so we can take Ek t \ real. The fields are then given by 



A(x) = W 



A=l 
2 



d 3 k 
l3l 



AT 7? I „—ikx | ikx 

Nk(-k,\ e + e 



fc°= 



f d k / 

E ( x ) = E / 79^3 Nk iWk ik ' x \ e 



A=l ' 
2 



(27T) 
d 3 k 



ikx gifcrr 



k°= 



/ - — 73 N k ik x e k ,\ ( e 
, _ X J (2vr) 



ikx ^ikx 



fe°= 



(B.45a) 
(B.45b) 
(B.45c) 



We may now express the total energy Q of the Maxwell field as the integral of (B.15) over 
some finite volume V 



Q 



d 3 xp 



[ d 3 xJ(E 2 + B 2 ) 
Jv * 



(B.46) 
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Inserting the Fourier expansion of the fields and taking V — > oo we get 
1 , 2 ^ ^ r r d 3 k f d 3 q 



ix(k-q) _|_ e —ix(k—q) _ e ~ix(k+q) _ e ix(k+q) 



k°=U) k 
q°=Uq 



The resulting delta functions are integrated over d 3 q. Recalling the vector identity (a x 
b) • (c x d) = (a • c)(b • d) — (a • d)(b • c), the expression simplifies further to 

1 v--v v—-, f d 3 k 

Q = / \3 ( w k4,A • 4,A' + k 2 e fc ,A • 4,A' ~ (k ■ 4,a) ■ (k • 4,A')) 



A=l A'=l ' 



Remembering the dispersion relation (B.41) and the polarization tensor definition ( |B.35[ ) 



2 2 q 2 2 q 

hkJ ^ ' hkJ ( 2 -) 3 



^AA' 



E 

A=l 



d 3 k 

"(2^ 



(27V fe 2 u; k ) Wk 



A similar calculation for the total momentum vector results in 



E 

A=l 



d 3 k 



(2iV 2 ^ k ) k. 



(B.47) 



(B.48) 



B.2.3 Plane wave expansion of the Proca field 



The vacuum solutions of the equations of motion for the Proca field (B.26) may also be 
written as plane waves 

d 3 k 



N k U k ^e- ikx + el^e 



ikx 



fc° = 



(B.49) 



The dispersion relation for the Proca field of mass m in vacuum results from replacing 
(IB. 491) into (IB. 261) and taking jt 1 = 0, 



v 7 ^ 



+ m z 



(B.50) 



Recalling the consistency condition d^A^ = and again choosing e k x ^ real this can be 
written as 



\ i J 



(2vr) d Wk 



d " k N k —k-e, x (e- ikx + e ikx 



A=l ' 
3 



k°= 



-u> k 



A=l ' 



fc4,A ( e ikx + e ikx 



k°= 



(B.51) 
(B.52) 
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The electric and magnetic fields expand as 



E(x) = J2 



3 JVfciWk I efc,A n — k ||e — e 



A=l (^) 

= 2J / ^ k x e ~M (e 

A— 1 ^ * 



at 



fc°= 



fc° = 



(B.53) 
(B.54) 



Then the total energy may be calculated as before from (B.28) to give 
3 



A= 



d 3 k 



(2JV£w k ) a*. 



(B.55) 



This is the same expression as (B.47) but with an added (longitudinal) polarization. 



The total momentum also follows from (B.29) after a little calculation, given by 

3 



P = E 



A= 



d 3 k 



(B.56) 



summing over all 3 modes A. Note that the total momentum it is always longitudinal, even 
though the gauge field and the electromagnetic fields are not always transverse. 



B.3 Energy density and Flux 

Recall that the energy-momentum tensor derived from the invariance of the Lagrangian 
under space-time translations obeys the conservation equation 



duT^ = , 



which integrated over a volume V in 3-space becomes 



= / d 3 xd M T^ 



(B.57) 



v 



d 3 x (d T 00 + diT 0i ) 



= d t / d d xp+ / d J xV • S. 
Jv Jv 

Applying the known Gauss theorem we recover the familiar equation of continuity 



d t 



/ d 3 xp+ <p i 

JV JdV 



da • S = 0. 



(B.58) 



which is valid for both Maxwell and Proca fields. 
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B.4 Transversal solutions for the Proca field 

It will become necessary for us to ignore the longitudinal modes of the Proca Equation since 
those modes are unexcitable in the thermal Yang Mills theory in the deconfining phase we 



reviewed in chapter |3j From (B.26) it is clear that we can choose transverse-only solutions 
for the sourceless case so long as we respect the condition d^A^ = 0. Choosing ^4° = 
immediately imposes the Coulomb "gauge" V • A = 0, making the gauge field transverse. 
The electric field also recovers transversality since V • E = -m 2 A° = 0. The de erees 
of freedom are now restricted to just two, and we recover all the same results from the 
Maxwell theory in Coulomb gauge but with a modified dispersion relation w k = Vk 2 + m? 

Q = J2 [ ( 2N i^) (B.59) 



A=l 

2 



(27TY 



P = E/7^3 ( 2 ^-k)k. (B.60) 



A=l 



(2*)' 



Consider now a single mode A = A k e * fcx |fco =aJk . The electric and magnetic fields are 
given by 

E k = iu; k A k , (B.61) 

B k = ik x A k . (B.62) 

The time averaged Poynting vector of this mode is given by 

S k = ^Re (E k x B k ) = ^a; k (A£k - A k • kA k ) 

= ^u; k |A k | 2 k, (B.63) 

|S k | = ^|E k | 2 . (B.64) 
2wk 

We will use these expressions in our discussions of antennas in chapter [5] 



Appendix C 



The Fresnel-Kirchhoff Diffraction 
Integral from the Proca Equations 



Here, we sketch a derivation of the Fresnel-Kirchhoff Diffraction Integral for the case of 
massive spin-1 particles. This integral formula is needed to calculate the effective antenna 
area, see section 5.2.2 We follow the presentation given in [60, 67J, but accounting for 



nonzero mass. 



C.l Helmholtz equations 



Take the curl of (B.23b) and eliminate B with (B.22b) to get an expression containing 



only E. Do likewise for B with (B.22b), we arrive at 

V x V x E + (d? + m 2 ) E = -da + V (V • E) , 
VxVxB+(3(|m 2 )B = Vxj. 



(C.la) 
(C.lb) 



Using the identity Vx VxV = V(V-V) - V 2 V and recalling V • B = we get finally 



(V 2 - 8 2 — m 2 ) E = da , 
(V 2 - d 2 t — m 2 ) B = -V x j 



(C.2a) 
(C.2b) 



Consider now the Fourier expansion of the components Ei, B^. It is clear that, in the 
absence of charges, each such Fourier mode E^, obeys the homogeneous Helmholtz 
Equations 



(V 2 + k 2 )£ k = 0, 
(V 2 + k 2 ) J B k = 0, 

with dispersion relation 

Wk = \/k 2 + m 2 . 

We have omitted the component indices for clarity. 

Finally, recall that the Green's function Gk(x,Xo) to the Helmholtz Equation 



(C.3a) 
(C.3b) 



(B.50) 



(V 2 + k 2 ) G k (x,x Q ) 



-4ir5 (z ) 



(C.4) 
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is given by the spherical wave expanding about a point xo 
-t|k[[x-xo| 

G k (x,x ) = -, r . (C.5) 

\X — Xq\ 

C.2 Integral theorem of Helmholtz and Kirchhofff 

We now investigate the problem of determining the field Vectors E, B at a point xq, given 
the values of the electric and magnetic field vectors over an equiphase surface S. 

Recall Green's theorem: Let two scalar functions F, G be defined over some volume V, 
further let them have continuous first and second derivatives in V. Then F, G satisfy the 
following identity 



f dV (FV 2 G - GV 2 F) + / dA(FVG - GVF) = . 

Jv JdV 



(C.6) 



Let E, B be defined over a volume V which contains xo, and have u k (x) stand for any 
Fourier mode E^, B^. Inserting and Gk as scalar potentials in (C.6), we gelQ 



v 



dV (G k (x, x ) V 2 u k (x) - u k (x)V 2 G k (x, x )) 

= -& dA(G k (x,x )Vu k (x) -u k (x)VG k (x,x )) . (C.7) 

JdV 



Using the identities (C.3), (C.4) 



Uk(xo) = -^-f cL4 (G k (x,x ) Vuk(x) - u k (i)VG k (i,i )) • (C.8) 
47r JdV 

The fields at point x$ have thus been expressed as an integral over a boundary surface dV 
which must contain xq. 

In application to antenna radiation problems, the surface dV (called antenna aperture) 
for which we know the field distribution (called illumination) is not closed. To calculate 



anything useful, formula (C.8) needs to be generalized to open surfaces. We expect linear 
integrals over the open surface boundary to appear as extra terms. However, it turns out 
that under still reasonable conditions we may ignore the extra terms and integrate directly 



with (C.8) over the open surface [60J. Namely, in the Frauenhofer (far-field) approximation 



the diffraction field over the antenna aperture A is given by 

Uk(xa) = --r[ ^(Gk (s, x ) Vu k (z) - u k (i)VG k (x, s )) . (C.9) 

The Frauenhofer approximation is valid provided the point xo is far enough from the 
aperture, 

R > -d 2 \k\, (CIO) 

7T 



lr The function Gk(i,io) is not continuous at a;o- In this particular case, use of Green's theorem may 
be made rigorous by a simple limiting process |65| 
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where R is the distance from xo to the aperture surface A, d the maximum linear dimension 
of the aperture, and k the wave vector of the monochromatic wave illuminating the aperture 



67] , Expression (C.10) defines the Frauenhofer region. Integral expression (C.9) allows the 



calculation of -u k in the Frauenhofer region, given the values o/u k at the aperture. Thus, we 



may use (C.9) to solve antenna aperture problems by replacing -u k with the corresponding 



electromagnetic field component £? k , i? k . 

Let us now approximate the field illuminating the aperture by using simple monochro- 
matic wave 



u k (x) = A{x) e ikL{x \ 

with amplitude A(x) and phase L(x). Then 

du k . „ ., . „ r IdA 
— — = n • Vu k = iku^n • VL + u k — — 
on A on 



(C.ll) 



(C.12) 



where n is defined as the normal to the surface element in (C.9), dA = cL4n. It can 



further be shown that s = VL is the normal unit vector to the propagation wavefront [60 
It defines the direction of propagation of the wave illuminating the aperture A. 
The second expression we need is 



dn 

where r = x 
get 



n- VG k 



Akr 



ik I n • f - 



xq. Finally, we define cos# r 



n • r, cos ( 



(C.13) 

n • s. After joining terms we 



u k (x,y,z) 



4tt 



dx a dy a 



ik ) cos 9 T — ik cos ( 



Akr 



-u k (x a ,y a ) 



(C.14) 



C.3 Fresnel-Kirchhoff diffraction integral 



Since our derivation of (C.9) depended only on the field u k obeying the Helmholtz equation, 
we may replace it with any solution of ( C.3 ) , including the vector components of the electric 
and magnetic fields. We can then write the Fresnel-Kirchhoff Diffraction Integral 



E k (x,y,z) 



1 

47T 



dx a dy a 



ik I cos 9 r — ik cos 6. 



ikr 



%a i ya 

). (C.15) 



Here E a ^(x a ,y a ) is the electric field at the aperture. In the Frauenhofer region, this 
expression can be further simplified to give 



E u (r, 



1 e 



( k r 



. |k| MM) . (c.16) 

Z7TZ r 

with r,9,(p polar coordinates defined with the aperture at the origin and normal to the 
z-axis |67j . The polar function /i k (9, <p) is defined as surface integral over the illumination 
at the aperture A 



hu(0,d>)= [ dx a dy a e- lksi » d ^ cos * + y« sin ^E a ^(x a ,y a ) 
J A 



(C.17) 
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In this region, the propagating modes are simply spherical waves, with the Poynting vector 
S parallel to k and f, S || k || f. The radial component is then read from (B.64) as 



1 



Ik| 2 |k| 



8tt 2 r 2 



\K (0,. 



(C.18) 



Appendix D 

Fundamental Radiometric Quantities 



Quantity 


Symbol / Definition 


Units 


Radiant Energy 


Q 




J (Joules) 


Radiant Energy Density 


w 


= dQ/dV 


J m~ 3 


Radiant Power (Flux) 


$ = 


= dQ/dt 


W (Watt) 


Radiant Exitance 


M = 


= d$/dA 


W m" 2 


Irradiance 


E = 


= d$/dA 


W m" 2 


Radiant Intensity 


I = 


= d$/dn 


W sr" 1 


Radiance 


L = 


= d 2 : $/ \dVtd A cos 9) 


W m~ 2 sr -1 


Emmisivity 


e = 


= M/M bb 





Taken from 1 22 
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